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Chapter  I of  this  report  contains  a brief  introduction  to  the  general 
theory  of  asymptotic  expansions,  and  serves  as  the  theoretical  back- 
ground for  the  main  part  of  the  report. 

In  Chapter  U,  the  most  important  methods  for  the  asymptotic  expan- 
sion of  functions  defined  by  definite  integrals  are  developed.  Double  and 
multiple  integrals  are  not  included. 

The  remaining  two  chapters  are  devoted  to  the  asymptotic  expansion 
of  functions  defined  by  ordinary  linear  and  homogeneous  differential 
equations  of  the  second  order.  In  Chapter  lH,  the  **large”  quantity  is  the 
variable  in  the  differential  equation,  and  the  only  case  discussed  in 
detail  is  that  of  a differential  equation  for  which  infinity  is  an  irregular 
singular  point  of  rank  one  with  a characteristic  equation  which  has  two 
distinct  roots.  In  Chapter  TV,  the  “large”  quantity  is  a parameter  in 
the  differential  equation,  and  the  variable  is  real  and  bounded.  Both 
Liouville’s  approximation,  and  its  generalization  appropriate  in  a trans- 
ition region  containing  a single  simple  transition  point,  are  discussed. 

The  report  is  based  on  a course  of  lectures  delivered  in  the  autumn 
of  1954  at  the  California  Institute  of  Technology.  The  author  wishes  to 
express  his  thanks  to  Mr.  C.  A.  Swanson  for  his  very  able  and  valuable 
assistance. 
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INTRODUCTION 

It  happens  frequently  that  a divergent  infinite  series  may  be  used  for 
the  numerical  compulation  of  a quantity  which  in  some  sense  can  be 
regarded  as  the  “sum”  of  the  scries.  The  typical  situation  is  that  of  a 
series  of  variable  terms  whose  “sum”  is  a function,  and  the  approxima- 
tion afforded  by  the  first  few  terms  of  the  series  is  the  better  the  closer 
the  independent  variable  approaches  a limiting  value  (often  oc  ).  In  most 
cases  the  terms  of  the  series  at  first  decrease  rapidly  (the  more  rapidly 
the  closer  the  independent  variable  approaches  its  limiting  value)  but 
later  the  terms  start  increasing  again.  Sucli  series  u'^ed  to  be  called 
semUconvergent  (Stieltjes),  and  numerical  computers  often  talk  of  con- 
vergently  beginning  series  (Lmde);  but  in  the  mathematical  literature 
the  term  asymptotic  series  (PoincareO  is  now  generally  used.  I'.e  shall 
see  later  that  asymptotic  series  may  be  convergent  or  divergent. 

Let  us  consider  an  example  first  discussed  by  Euler  (1754).  The  series 

(1)  S(x)=  1-  l!x  + 2!  x*-3!*^+  =1  (-l)"n!x" 

0 

is  certainly  divergent  for  all  * / 0,  yet  for  small  x (say  10  *)  the  terms 
of  the  series  at  first  decrease  quite  rapidly,  and  an  approximate  numerical 
value  of  S (x)  may  be  computed,  ^hat  function  of  x does  this  numerical 
value  represent  approximately? 

Euler  considers  <^(x)=  xS(x).  Then 

<;i,'(x)=  11- 2!  X + 3!x*  - — = 


or 

X*  + 0(x)=  X, 
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and  <p(*)  may  be  obtained  as  that  solution  of  this  differential  e<7uation 
which  vanishes  as  x 0.  Alternatively,  we  use  Euler’s  integral  of  the 
second  kind, 

n ! = / dt 

0 

and  obtain 


S(x)=f  e~^dt-xf  e~*tdt  + x^j  e~*  dt  - 
0 0 "*0 


= 2 (-1)"  f"e-^(xt)''dt. 

If  we  formally  sum  under  the  integral  sign,  S (x)  becomes 


(2) 

Now, 
(3)  fix) 


OO  —( 

e 


dt. 


+ xt 


- jIL 

Jo  1 + 


dt 


is  a well-defined  function  of  x,  as  a matter  of  fact  an  analytic  function 
of  X in  the  complex  x-plane  cut  along  the  negative  real  axis,  and  it  is 
closely  related  to  the  so-called  exponential  integral.  The  question  then 
arises:  in  what  sense  does  the  divergent  series  (1)  represent  the  func- 
tion (3)?  To  answer  this  question,  we  note  that  for  m = 0,  1,  2,  ... 


1 

\ + xt 


i-Xt)'*' 

1 + xf 


and  hence 

(4)  fix)  = S ^ix)  + R ^ix) 


where 


(5)  S (x)=  i (-l)"n!x" 

n=  0 
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is  a partial  sum  of  (1),  and 


(6) 


R 

m 


(x)  . i-x) 


e~* 


1 + xt 


dt 


is  the  remainder > 

If  Re  * > 0,  we  have  |1  + *t|  < 1 and 

(7)  |f?.(x)|  <(m+  D!  Rex>0. 

On  the  other  hand,  if  Re  a:  < C,  4>  = arg  *.  and  7r/2  < ± < n,  then 

|1  + xt\  ’ < jcosec  <j)\. 


and 


(8)  I/?,(*)|<(m  + l)!|*|“-"'  |cosec0|.  He  a:  < 0. 

In  either  case,  the  remainder  is  of  the  order  of  the  first  “neglected” 
term  of  Six),  and  approaches  0 rapidly  as  x -»  0.  The  limit  is  uniform  in 
any  sector  |arg  x|  < rr  ~ f > 0.  If  Re  x > 0,  the  remainder  is  numeri- 
cally less  than  the  first  neglected  term,  and  ifx  > 0,  the  remainder  has 
also  the  sign  of  the  first  neglected  term.  Thus,  for  x > 0,  the  series  (1) 
behaves  very  much  like  a convergpnt  alternating  series,  except  that  the 
smallest  term  of  (1),  which  occurs  when  m is  approximately  equal  to  x, 
determines  a limit  to  the  accuracy  beyond  which  it  is  impossible  to 
penetrate. 

The  theory  of  asymptotic  series  was  initiated  by  Stieltjes  (1886)  and 
Poincare'  (1886).  I^e  may  distinguish  two  parts  of  the  theory.  One  part, 
which  we  may  call  the  theory  of  asymptotic  series,  treats  topics  such  as 
“sums”  of  asymptotic  series  (“asymptotic  limits”,  “asymptotic  con- 
vergence”), and  operations  with  asymptotic  series  (algebraic  operations, 
differentiation,  integration,  substitution  of  asymptotic  expansions  of  a 
variable  in  convergent  or  asymptotic  series  involving  this  variable,  and 
the  like).  The  most  comprehensive  presentation  of  this  part  of  the  theory 
is  to  be  found  in  van  der  Corput’s  Lectures  (1951,  1952)  and  current 
publications  by  the  same  author.  In  these  pages  we  shall  restrict  our- 
selves to  a brief  introduction  to  the  theory  of  asymptotic  series,  and  shall 
devote  most  of  our  attention  to  the  other  part  of  our  subject,  to  the  theory 
of  asymptotic  expansions.  Here  the  central  theme  is  the  construction  and 
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investigation  of  series  which  represent  given  functions  asymptotically. 
The  functions  arc  often  given  by  integral  representations,  or  by  power 
series,  or  else  appear  as  solutions  of  differential  equations;  and  in  the 
latter  case  the  “variable”  of  the  asymptotic  expansions  may  occur  either 
as  the  independent  variable,  or  else  as  a parameter,  in  the  differential 
equation. 
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CHAPTER 1 

ASYMPTOTIC  SERIES 


1.1.  O- symbols 

In  general,  the  “independent  variable”  will  be  a real  or  complex 
variable,  but  in  this  chapter  x stands  for  a variable  element  of  a topo- 
logical Tj-space  (liausdorff  space)  except  when  stated  otherwise.  The 
variable  x ranges  over  a set  R,  and  x^  is  a limit  point  of  R (which  may 
or  may  not  belong  to  R).  <f>(x\  xpix),  and  similar  symbols  denote  real-  or 
complex-valued  numerical  functions  of  x defined  when  x is  in  R. 

The  following  order  relations  involving  the  order  symbols  O,  o will  be 
used.  We  write  <f>  = in  R if  there  exists  a constant  (i.e.,  number 

independent  of  x)  A so  that  |<;6|  <A\i!f\  for  all  x in  /?;  0 = O (i/i)  as  x -*  x^ 
if  there  exists  a constant  A and  a neighborhood  U of  *jjSO  that  j<^|  < A 
for  all  X common  to  U and  R;  and  we  write  0 = o (0)  as  x -♦  if  for  any 
given  f > 0 there  exists  a neighborhood  1/ ^ of  x^  so  that  |0|  < f |0|forall 
X common  to  and  /?.  If  0 ^ 0 in  then  the  three  conditions  may  be 
formulated  more  simply  :0  = 0(0)  in  R [as  x in  /?]  if  0/ 0 is  bounded 
in  R '[as  x ->  x^  in  /?],  and  <j)  = o (0)  as  x -»  x^  if  0/  0-»O  asx-»x^. 

In  the  following  examples  x is  a complex  variable,  and  is  the 
sector  0 < |x|  < oo:,  |arg  x|  < n/2  — A.  The  reader  should  verify  each 
assertion,  (i)  = 0(x“),  e“*  = ofx*^)  as  x -»  oo  in  S^,  A > 0,  a arbitrary; 

and  neither  of  these  order  relatione  holds  (for  arbitrary  a)  when  A < 0. 
(ii)  e~*  - 0{x‘^)  as  X -*  oo  in  5^  provided  that  Re  a > 0;  this  order  relation 
fails  to  hold  when  Re  a < 0.  (iii)  e”*  = 0(x“)  in  provided  that  either 
A > 0 and  Re  a < 0 or  A = 0 and  Re  a = 0 . 

If  the  functions  involved  in  an  order  relation  depend  on  parameters,  in 
general  also  the  constant  A,  and  the  neighborhoods  V , V ^ involved  in 
the  definitions  will  depend  on  the  parameters.  If  A,  U,  U ^ may  be  chosen 
to  be  independent  of  the  parameters,  the  order  relation  is  said  to  hold 
uniformly  in  the  parameters. 
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1.1 


J 


Operations  with  order  relations  are  governed  by  a number  of  simple 
rules.  We  shall  set  out  the  more  frequently  used  rules  for  the  0-syrahol: 
the  corresponding  rules  hold  for  the  o-symbol.  In  the  following  rules  R 
snd  Xg  are  fixed,  and  the  qualifying  phrase  “as  x -*  x^”  is  omitted 
throughout. 

If  0 = 0{il/)  and  a > 0 then 

(1)  i<j!>r=o(i,An. 

If  0 ^ = O ( j),  i = 1 , ...  , A and  the  a ^ are  constants,  then 

(2)  2 a.<;6,=  0(2ia.l  I0J). 

i 

This  relation  holds  also  for  infinite  series  provided  that  ^ ^ 

uniformly  in  t.  In  the  case  of  infinite  series,  equation  (2)  and  similar 
statements  will  be  interpreted  in  the  following  manner.  If  con- 
verges then  so  does  2 and  (2)  is  true,  and  if  2 ja..  i/f  J diverges 

then  there  is  nothing  to  state. 

If  ^ = 0((/f  ^),  t = 1,  ...  , k,  the  are  constants,  and  |i/»  J < t/f  for 
i = 1,  ...  , k and  for  all  x common  to  R and  to  some  neighborhood  V ^ of 
Xq,  then 

(3)  2aj0j=O(^). 

This  relation  holds  for  infinite  series  provided  that  4*  i~  0{4>  uniformly 
in  i,  and  2 |<i  .|  < «>. 

If  <^  . = 0(0  p,  I = 1,  , k then 

(4)  n0.=  O(ndr,). 

i i 

The  proof  of  (1)  is  immediate.  To  prove  (2), we  remark  that  by  assump- 
tion there  are  numbers  A . and  neighborhoods  U ^ of  x^  associated  with 
the  ..  If  the  number  of  the  0^is  finite,  there  is  an  A larger  than  all  the 
A and  a neighborhood  U contained  in  all  the  U and 

|2  j < 21a -1  d .10  < A21a -1  10-1 

when  X is  common  to  R and  u,  and  this  proves  (2).  If  there  is  an  infinite 
number  of  0 then  the  existence  of  A and  U follows  from  the  uniformity, 
in  i,  of  the  order  relation.  (3)  can  be  deduced  from  (2)  since  under  the 
circumstances  envisaged  we  may  take  U above  to  be  contained  in  U ^ and 
then 
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where  <4,  = ^4  2 |a  J is  a finite  constant.  The  proof  of  (4)  is  similar  to  that 
of  (2). 

Order  relations  may  be  integrated  either  with  respect  to  the  inde- 
pendent variable  orwith  respect  to  parameters.  For  the  sake  of  simplicity 
we  shall  restrict  ourselves  to  integrals  with  respect  to  real  variables. 
Extensions  to  complex  and  abstract  variables  are  possible. 

Let  X be  a real  variable,  let  R be  the  interval  a < x < b,  and  let 
0 = Oid')  as  X -*  b.  If  0 and  0 are  measurable  in  R then 

(5)  / 0(0  dt  = <?(/'’  10(01  di)  as  X -*  b. 

Proof:  If  / * |0(t)l  dt  = oo , there  is  nothing  to  prove.  If  / |0|  dt  < oo 

X X 

for  some  x,  then  A and  X exist  so  that  / **  |0|  dt  < oo  and  |0(x)|  ^.4  l0(x)| 
for  X < X < b,  and  hence 

l/*"  0(t)  dt|  < |0(t)l  dl  <A  / **  |0(t)|  dt  for  X < x < b. 

Let  X be  a variable  element  of  the  set  /?  in  a liausdorff  space,  let  y 
be  a real  parameter,  a < y < /3,  and  let  0(x,  y)  = 0(0  (x,  y)),  uniformly 
in  y,  as  X -►  Xg.  If  for  each  fixed  x in  /?,  0 and  0 are  measurable  func- 
tions of  y in  a < y < ^ then 

(6)  0(x,  y)  dy  = Oi/J^  I0(x,  y)|  dy)  as  x -»  x^  • 

The  proof  is  similar  to  that  of  (5).  On  account  of  the  uniformity  of  the 
0-symbol,  A and  V are  independent  of  y,  |0|  < A j0|,  and  (6)  follows  by 
integration  of  this  inequality  with  respect  to  y. 

It  is  in  general  not  permissible  to  differentiate  order  relations  either 
with  respect  to  the  independent  variable  or  with  respect  to  parameters. 
However,  some  general  results  on  the  differentiation  of  order  relations 
exist  in  the  case  of  analytic  functions  of  a complex  variable  (see 
sec.  1.6). 
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1.1 


We  conclude  this  section  with  a few  formulas  concerning  combinations 
of  order  relations 

(7)  0(0(0))=  0(0) 

(8)  O(o(0))  = o(O(0))=o(o(0))  = o(0) 

(9)  0(0)  0(0)=  0(00) 

(10)  O (0)  O (0)  = O (0)  O (0)  = O (00) 

(11)  O(0)  + O(0)-=O(0)  + o(0)  = O(0) 

(12)  o (0)  + o (0)  = o (0) 

The  { roof  of  these  formulas  is  immediate,  and  they  can  be  extended  to 
combinations  of  any  finite  number  of  order  symbols. 

In  referring  to  the  above  rules  we  shall  quote  the  number  of  the  equa- 
tion which  expresses  th  final  conclusion,  and  we  shall  use  the  same 
number  to  indicate  the  corresponding  rule  for  the  o-symbol.  For  instance, 
(1)  will  indicate  either  the  rule  that  0 = 0(0)  and  a > 0 imply  |0|  “=0  (|01  ‘^) 
or  the  rule  that  0 = o(0)  and  a > 0 imply  101“^=  o(101“). 

1J2.  yVsymptoCic  sequences 

In  this  section  R,  x,  x^,  <f)  have  the  same  meaning  as  in  sec.  1.1.  A 
finite  or  infinite  sequence  of  functions,-  0j,  0^,,  ...  , will  be  abbreviated 
as  10„l. 

The  sequence  of  functions  I0„}  is  called  an  asymptotic  sequence  for 
X -»  Xq  in  R if  for  each  n,  0^  is  defined  in  R and  0„+,  = as  x •*  x^ 

in  R. 

If  the  sequence  is  infinite  and  0^  + | = uniformly  in  n,  then  10^1 

is  said  to  be  em  asymptotic  sequence  uniformly  in  n.  If  the  0^  depend  on 
parameters  and  0„  + | = o (0^)  uniformly  in  the  parameters,  then  10^^!  is 
said  to  be  an  asymptotic  sequence  uniformly  in  the  parameters. 

We  proceed  to  give  some  examples  of  asymptotic  sequences  in  which  * 
is  a complex  variable,  R is  the  complex  plane  except  when  otherwise 
specified,  and  is  the  sector  defined  in  sec,  1.1. 

(i)  {(x  - Xj,)"!,  X -*  Xg-, 

(ii)  |x”"l,  X -»  oo  ; 

_ -A 

(i-Ii)  lx  "j,  X -*  oo  in  S^, 
whore  Re  > Re  for  each  n; 

(iv)  lx  "I,  X -*  oo, 

A real  and  A +,  > A,  for  each  n; 

fl  ▼ 1 n 
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(v)  \e*  X "i 

and  X and  either  as  in  (iii)  or  as  in  (iv); 

(vi)  "},  at oo  in 

and  either  A > 0 and  the  A are  as  in  (iii)  or  A > 0 and  the  A are  arbitrary 

(vii)  ir(x )/!'(*+  n)|,  X -►  oc  in  5^,  A > -77/2. 

The  reader  should  verify  that  each  of  the  sequences  (i)  to  (vii)  is  an 
asymptotic  sequence,  and  should  justify  the  restrictions  imposed  upon  A 
and  A^  in  these  examples.  Why  is  (iii),  with  arbitrary  Im  A^,  not  an  asym- 
ptotic sequence  for  x -*  oo  in  the  complex  plane  (without  restriction  to 
sonje  S^)?  The  infinite  sequence  lF(x  — n)/F(x  + n)},  n = 1,  2,  ...  is 
not  an  asymptotic  sequence  for  x -►  oo  in  any  region  including  unbounded 
portions  of  the  real  axis,  but  it  is  an  asymptotic  sequence  for  x -►  oo  in 
any  region  whose  closure  lies  entirely  in  the  upper  or  the  lower  half- 
plane. The  finite  sequence  ir(x  — n)/r(x  + n)},  n = 1,  2,  ...  , A’,  is  an 
asy'mptotic  sequence  for  x -♦  oo  in  any  R. 

From  given  asymptotic  sequences  new  such  sequences  may  be  ob- 
tained by  processes  which  are  largely  based  on  the  operational  rules  of 
order  symbols  given  in  sec.  1.1.  In  describing  some  of  these  processes 
we  shall  restrict  ourselves  to  real  variables  although  extensions  to  more 
general  variables  are  possible.  In  most  cases  x^  and  R will  not  be 
mentioned:  in  such  cases  they  are  fixed. 

Any  subsequence  of  an  asymptotic  sequence  is  an  asymptotic  sequence. 
The  proof  follows  from  1.1  (8). 

is  an  asymptotic  sequence  and  a > 0,  then  1 is  an  asym- 

ptotic sequence.  The  proof  follows  from  1.1  (1). 

Two  sequences,  lo  t and  liA  i,  so  connected  that  ch  = Oiih  ) and 
t/<  = for  each  n,  are  said  to  be  equivalent.  and  are 

equivalent  sequences  and  10,^1  is  an  asymptotic  sequence,  then 
also  an  asymptotic  sequence.  To  prove  that  is  asymptotic  we  remark 
that 


+ i = C>(0„  + ,)=  0(o(9i„))=  O(o(O(0„)))=  o(0j 
by  1.1  (8). 

If  10  I and  10„1  o.re  asymptotic  sequences  containing  the  same  number 
of  functions,  then  10  0 1 is  an  asymptotic  sequence.  The  proof  follows 

from  1.1  (10). 
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n = 1,  ...  , N,  is  an  asymptotic  sequence,  n »=  1,  ...  , A’; 

* =0,1,  , k < N is  a set  of  positive  constants,  ^<  a.^  ^for  all 

n,  i and 


(1) 


ra  = 1, 


...  , 


N -k, 


then  10^1  is  an  asymptotic  sequence. 

In  this  statement  N may  be  finite  or  infinite;  h is  finite.  To  prove  th.e 
statement  we  remark  that  it  follows  from  k being  finite  that  for  any  n and 
any  f > 0 there  exists  a neighborhood  V ^ of  so  that  + < f|s6j  in 

the  common  part  ol  U ^ and  R for  r = n,  n + 1,  ...  , n +•  i.  We  then  have 


(2)  lA. 


+1 


k 

= 2 
i=  0 


a + t +1 1 — 


S a .1^  +J  = f 

i=0  n.a‘~n  + * n 


The  extension  to  infinite  sums  is  contained  in  the  following  theorem. 

Let  he  an  asymptotic  sequence  uniformly  in  n,let  1,  ...  , 

i >■  0,  1,  ...  be  a set  of  positive  constants  suck  that  + , ^for  all 

n,  i and  put 


(3) 


“n.  . + 


n « 1,  2,  ... 


If  the  infinite  series  for  i/f|  converges  in  some  neighborhood  of  x^,  then 
there  is  a subset  of  R so  that  is  a limit  point  of  R^,  all  infinite 
series  (3)  converge  in  R^,  and  is  an  asymptotic  sequence  for  x -*  x ^ 
in  R^,  uniformly  in  n. 

Proof:  From  the  uniform  asymptotic  property  of  10^},  it  follows  that  there 
exists  a subset  /?,  of  R so  that  is  a limit  point  of  R , and  + , 
for  all  X in  f?  , and  all  n . For  x in  /?  , 


^ + + t + “n.  < -•  < ^ i + , I 


SO  that  all  infinite  series  (3)  are  dominated  by  the  series  for  0 , . If  the 
series  for  converges  in  a subset  R^  o{  R,  and  x^  is  a limit  point  of 
we  take  R^  to  be  the  common  part  of  /?  , and  R All  functions  x!j ^ 
are  defined  in  R q,  x^  is  a limit  point  of  R and  on  account  of  the  uniforn 
asymptoticproperty  of!  <;&„!,  equation  (2)  with  k = oo  holds  uniformly  in  n. 
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New  asymptotic  sequences  may  be  formed  by  integration  in  two 
different  ways. 

y)\  is  an  asymptotic  sequence  uniformly  in  y,  a <v  < /',  for 
S'  -♦  in  R,  and  if  all  integrals 

(4)  «!>„(*)  = y)| 

exist,  then  14>  } is  an  asymptotic  sequence . The  proof  follows  froir 
1.1  (f*).  As  in  the  case  of  (3),  it  is  sufficient  to  assiinie  tlial  all  q>^{x,  y) 
are  measurable  functions  of  y,  and  that  is  integrable.  The  intepra- 
bility  of  all  (f> possibly  for  e more  restricted  set  R then  follows  by 
showing  that  /|0^|dy  is  dominated  by  dy. 

If  X is  areal  variable,  R is  the  interval  a < x < b,  is  an  asymptotic 

sequence  for  x -*  b,  and  if  atl  integrals 

(5)  «!)„(*)  = dt 

exist,  then  is  an  asymptotic  sequence  for  x -*  b.  The  proof  follows 
from  1.1  (5),  and  it  is  again  sufficient  to  assume  that  a\\  are  n;easur- 
able  and  is  integrable;  the  result  tlien  follows  at  least  for  some  inter- 
val a , <x  < b. 

Note  that  the  differentiation  of  an  asymptotic  sequence  does  not 
necessarily  yield  an  asjTnptotic  sequence.  For  instance,  tahe 

p = a:~"[a  T cos(*")]  n = 1,  2,  ...  , 

' n 

Then  1<^„I  is  an.  asymptotic  sequence  for  x -♦  oo  on  the  real  axis,  but 
is  not  an  asymptotic  sequence. 

1.3.  Asymptotic  expansions 

In  this  section  and  in  the  following  sections,*,  have  the  same 

meaning  as  in  sec.  1.1;  l^„l»  ^Xn^»  *'•  » always  asymptotic 

sequences  for  X -*  x^  in  R;  f{x),  g(x),  h (x),  ...  arc  numerical  functions 
of  X defined  in  R;  and  a,  b,  c,  ...  are  constants  (i.e.,  independent  of  x). 

The  (formal)  series  2 a^  said  to  be  an  asymptotic  expan- 

sion to  N terms  of  f{x)  as  x x^  if 

(1)  f{x)~  2 a„  p„(x)+ 0(9!)^)  as  x-*x^. 
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An  asymptotic  expansion  to  N terms  will  often  be  indicated  as 

(2)  fix)  ~ S ^ terms  as  * -*  in  R , 

and  the  qualifying  phrase  “in  /?’’  will  often  be  omitted.  An  asymptotic 
expansion  to  1 term  will  be  written  as 

(3)  fix)  ~ a,  <;6,(x)  * - *0 

and  will  be  called  an  asymptotic  representation-,  and  an  asymptotic 
expansion  to  any  number  of  terms  (i.e.,  witb  N — «>)  will  be  written 


(4)  /(x)  ~ 2 a 6 (,x)  ^ *n 

and  called  an  asymptotic  expansion.  An  asymptotic  expansion  may  be 
convergent  or  divergent.  In  most  textbooks  only  the  cases  N = \ and 
A'  = oo  are  discussed,  but  we  shall  let  N stand  for  any  positive  integer. 

If  an  asymptotic  expansion  to  N terms,  with  N finite,  involves  certain 
parameters,  we  shall  say  that  it  holds  uniformly  in  these  parameters  if 
the  remainder  in  (1)  in  oi<f>f^)  uniformly  in  the  parameters.  An  asymptotic 
expansion  ij^  = oo)  involving  certain  parameters  will  be  said  to  hold 

uniformly  in  these  parameters  it  f — ^ = o(0^)  uniformly  in  the 

#»«  I 

parameters  tor  each  sufficiently  large~^  (but  not  necessarily  uniformly 
in  M.). 

The  formal  (finite,  or  infinite)  series  So  0 will  be  called  an 

_A  n F» 

asymptotic  series.  If  0^  — x ",  we  shall  speak  of  an  asymptotic  series 
of  powers,  and  if  0^  = x*",  of  an  asymptotic  power  series.  For  instance, 
2 (a  — 1)!  (— x)"~'  is  the  asymptotic  power  series  expansion  for  x -*  0 in 
of  the  function  fix)  defined  by  equation  (3)  of  the  Introduction. 
Some  authors  speak  of  asymptotic  power  series  when  0„=^o(*.)  * 
but  it  is  more  appropriate  to  call  the  series  divided  by  ff^ix)  an  asym- 
ptotic power  series. 

From  (1)  it  follows  that  the  coefficients  in  an  asymptotic  expansion 
to  N terms  may  be  computed  by  means  of  the  recurrence  formula 

(5)  a = lim  l[/’(x)  - *2  0^(x)]/0_  (x)|  m = 1,  ...  , /V 

" *-*  *0  * 


✓ 
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Conversely*  suppose  that  we  have  + 1 functions, 
fix),  <^,(x),  ...  , <f>ffix) 

defined  in  R.  If  (5)  holds  and  ^ 0 for  m = 1,  , N then  10^1  is  an 

asymptotic  Sequence  for  x -*  x^,  and  S “„<?!>„  is  an  asymptotic  expansion 
to  N terms  of  fix)  as  x -*  x 

To  prove  that  is  an  asymptotic  sequence  we  have  to  show  that 
4‘ ^ +1  = oi4>J)  for  m = 1,  ...  , — 1.  Now  from  (5), 

f-  ^ 

ns  1 

and  if  we  replace  m by  m + 1 in  (5)  we  have 

f-  ^ <^-+1 

ns  1 

Comparing  the  last  two  equations  we  find 
+0(1)]?!',+,  = oi4>^). 

If  a ^ ^ 0 then  a,  + o (1)  0 for  x in  some  neighborhood  of  x^,  and 

we  may  divide  by  this  factor  to  see  that  0,,+,  “ oi<f>^).  Thus  is 

an  asymptotic  sequence.  Moreover,  (5)  with  m = N shows  that  (1)  holds, 
and  2)  is  an  asymptotic  expansion  to  A'  terms  of  f. 

If  S a^  (f>^ix)  is  an  asymptotic  expansion  to  N terms  of  fix),  then 
the  same  formal  series  will  also  provide  an  asymptotic  expansion  to 
any  lesser  number  of  terms  of  the  same  function.  l\e  also  have  the 
somewhat  sharper  result 

(6)  fix)=  2 ,V  = 1,  ...  , A^  - 1 

n = 1 

which  is  an  immediate  consequence  of  (1). 

\\ith  Xq  and  R fixed,  (5)  shows  that  the  asymptotic  expansion  to  a 
given  number  of  terms  of  a given  fimction  is  unique  if  the  asymptotic 
sequence  is  gi^'en.  On  the  other  band,  one  and  the  same  function  may 
have  asymptotic  expansions  involving  two  different  asymptotic  se- 
quences, and  the  two  sequences  need  not  be  equivalent  in  the  sense  of 

• TTjis  theorem  was  suggested  hy  Dr.  A.G.  Mackie. 
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sec.  1,2.  For  instance, 

— 2 (-!)"•'  *■" 

1 + X 

1 X 

-J_  (x*-r  > l)x-*" 

1 + X 


In  this  example  all  three  asymptotic  expansions  are  convergent  series 
when  |x|  > 1.  It  often  happens  that  some  asymptotic  expansions  of  a 
function  diverge  while  others  converge.  The  transformation  of  divergent 
asymptotic  expansions  into  convergent  ones  is  of  great  analytical,  although 
of  very  little  computational  interest.  Transformations  of  asymptotic  expan- 
sions into  convergent  expansions  or  else  into  expansions  more  suited  to 
numerical  com.putation  have  been  investigated  among  others  by  Airey 
(1937),  van  der  Corput  (1951),  ^■iiller  (1952),  van  Wijngaarden  (1953), 
Watson  (1912  b). 

An  asymptotic  expansion  does  not  determine  its  “sum”,  fix),  uniquely. 
For  instance,  the  functions  (1  + a:)”’,  (1  + e~*)/(l  + x),  (1  + + x)”' 

all  possess  the  asymptotic  expansion  2(— 1)"”'  x”"  as  x -►  » in  S^, 
A > 0.  A given  (finite  or  infinite)  asymptotic  sequence,  l^^t,  for  x -♦  x^ 
in  R establishes  an  equivalence  relation  among  functions  defined  in  /?: 
fix)  and  g (x)  are  asymptotically  equal  v/ith  respect  to  if 


fix)  - gix)  = oi<f) ^)  as  x-»Xg  in  R, 

for  all  n occurring  in  the  sequence.  An  asymptotic  series  represents  a 
class  of  zisymptotically  equal  functions  rather  than  a single  function. 


1.4.  Linear  operations  with  asymptotic  expansions 

If  f'^'S.a^rf^  and  g ~ 2 6^  both  to  N terms,  and  if  a,  j8  are 

constants,  then 

(1)  a/(x)  + /3g(x)  ~ 2 (aa^  + /3  6^)  0^(x)  to  N terms. 

The  proof  of  this  theorem  is  obvious,  as  is  its  extension  to  a linear 
combination  of  an  arbitrary  finite  number  of  asymptotic  expansions.  The 
extension  to  an  infinite  series  of  asymptotic  expansions  is  as  follows: 
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If  f iM  2 ^ terms,  uniformly  in  i,  t = 1,  2,  ...  , and  if 

the  a .are  constants  for  which  2 a .converges  absolutely  and 


(2) 


eo 

" i=  1 


a .a . 

n,  I t 


converges  for  each  n,  then  2 a ^f  ^x')  converges  in  some  neighborhood  of 
*o» 

(3)  f‘{x)=  2 a.f  d>  {*)  to  N terms. 

I't  n ’ n 


Proof:  We  have 

fi-  ^ 

fl  = 1 

uniformly  in  i,  and  £|a  J < oo.  By  1.1(3), 


2 


«=  1 


and  the  infinite  series  on  the  left  is  convergent  at  least  in  some  neigh- 
borhood of  x^.  Adding  ^ <f>^  to  both  sides  we  have  (3)  when  N < oo, 

n=  1 

If  /\'  = oo  then  both  the  assumptions  and  the  conclusion  hold  for  all 
sufficiently  large  N,  and  hence  2/4  <f>  is  an  asymptotic  expansion  of 

F (x). 

More  generally,  we  may  extend  (1)  to  finite  or  infinite  asymptotic 
expansions. 

Let  \<f>  I,  n = 1,  ...  , A < oo  and  m = , M < x be  asymptotic 

sequences  for  the  same  R,  x^  ; and  let  = 0{ib„)  for  each  m:  if 
f a d>  to  N terms  and  for  each  n,  (b  ~ 2 6 ih  to  hi  terms,  then 

j n ' n ^ ’ft  Bn  ' tt 


(4)  /(x)  ~ 2 (x)  to  , If  terms. 


where 
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Let  10^1,  /I  = 1,  2,  , and  \*p ^ },  'll  = 1,  < oo  be  asymptotic 

sequences;  suppose  that  for  each  n there  is  an  integer  fi{n)  < M so  that 
p(,n)  -*  M as  n -»  oo  and  if  f I.  a ^ 

to  M terms  uniformly  in  n,'S,  a ^ is  absolutely  convergent,  and  the  infinite 
series  in  (5)  is  convergent  for  each  m;  then  (4)  holds. 

The  proof  for  M,  N finite  is  immediate,  since  then 

/■(*)=  <f„ 

M 

= 2 c d/  + o (i/f  ) 

• « 1 " 

by  1.1(3).  If  A/  = oo,  then  the  same  reasoning  holds  for  any  M,  and  hence 
S c^  is  an  asymptotic  expansion  (with  M = oo)  of  fix).  Ir  the  exten- 
sion to  /V  = oo  we  use  the  extension  of  1.1(3)  for  infinite  series. 

We  now  turn  to  the  integration  of  asymptotic  expansions  either  with 
respect  to  a real  parameter  y,  or  with  respect  to  the  variable  x.  In  the 
latter  case  x will  be  assiuned  to  be  a real  variable. 

If  fix,  y)  ~ 2 o„(y)  ^„(*)  ^ terms,  uniformly  in  y,  a < y < /3,  if 

fix,  y),  for  each  fixed  x,  and  a„(y)»  for  each  fixed  n,  is  a measurable 
function  of  y,  and  if  hiy)  is  an  integrable  function  of  y for  which  each  of 
the  integrals 

(6)  A^  = ff  hiy)a^iy)  dy 
exists,  then  also  the  integral 

(7)  Fix)=  hiy)  fix,  y)dy 

a. 

exists  for  each  x in  some  neighborhood  of  x^,  and 

(8)  F(*)~2/4  4>  ix)  to  N terms. 

The  proof  is  very  similp.  to  that  given  above  for  infinite  series, 
except  that  1.1(6)  must  be  used  instead  of  1.1(3).  Some  generalizations 
of  this  theorem  are  obvious:  the  interval  (a,  ^)  may  be  replaced  by  any 
measurable  set,  of  finite  or  infinite  measure,  and  there  is  a similar  result 
for  multiple  integrals. 
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Let  X be  a real  variable,  let  R be  the  interval  a < x <.  b,  let  \cb  (*)l  be 

. ' n 

an  asymptotic  sequence  of  positive  functions  for  x -*  b,  and  assume  that 
each  of  the  integrals 

(9)  <D  (r)  = /'’ 

n X ” n 

exists.  If  fix)  ~ 2 <p  ^{x)  to  A'  terms  as  x -*  b,  and  fix)  is  a measurable 

function,  then 

(10)  F(*)=  J**  fit)dt 

X 

exists  in  some  interval  c < x < b,  and 

(11)  Fix)  ~ 2 a (■»:)  to  N terms  as  x b. 

n A 

The  proof  follows  from  1,1  (5). 

It  is,  in  general,  net  permissible  to  differentiate  asymptotic  expan- 
sions either  with  respect  to  the  variable  x,  or  with  respect  to  parameters. 
Some  general  results  on  the  differentiation  of  asymptotic  expansions  of 
analytic  functions  of  a complex  variable  exist  and  will  be  given  in  sec. 

1.6. 


1.5.  Other  operations  with  asymptotic  expan.sioiis 

Multiplication  of  asjTrptotic  series  does  not  in  general  lead  to  an 
asymptotic  series,  for  in  the  formal  product  of  2 <f>  ^ and  ^ b ^<f>  ^ all 

products  (f>  occur,  and  it  is  in  general  not  possible  to  arrange  the 
system  of  functions  1^^  m,  n = 1,  ...  , A so  as  to  obtain  an  asym- 

ptotic sequence.  There  are,  however,  important  special  asymptotic 
sequences  \ \ with  the  property  that  the  products  q> ^ rp ^ either  form  an 

asymptotic  sequence,  or  else  possess  asymptotic  expansions  in  terms  of 
an  asymptotic  sequence  (which  need  not  be  1<^  1).  First  we  shall  prove  a 
general  result  on  the  multiplication  of  two  asymptotic  expansions. 

Let  10^1,  n = 1,  ...  , A,  10^1,  m = 1,  •••  . hf,  and  k = , K be 

three  asymptotic  sequences  such  that  0,  = O(x^).  0i  = Oix^^), 

and 

(1)  c^j  X*  to /C  terms. 
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If  f 2 a ^ cf>^  to  N terms  and  g ^ b ^ to  M terms,  then  fg'^'S. 
to  K terms,  where 


/V  N m 

(2)  C.=  2 2 

n = 1 *=  1 


nmk  * 


Here  K may  be  finite  or  infinite;  N,  M are  finite.  The  result  remains  true 
if  M,  or  A,  or  both  M and  l\'  are  infinite,  provided  that  each  of  the  in- 
finite series  {or  double  series)  in  (2)  converges. 

The  coefficients  C are  those  obtained  upon  multiplication  of  2 a^  cf>  ^ 
and  2 6^  4i  ^ and  substitution  of  (1)  so  that  instead  of  (2)  we  may  say 
“where  the  coefficients  are  obtained  by  formal  substitution”,’  and 
this  description  will  be  used  in  similar  cases  throughout  this  section. 

We  first  prove  the  theorem  for  finite  N,  f.i,  K. 

/g=[  2 + )][  2 (A,  +o(«/r  )] 

n s I «=  1 


ns  1 «s  1 

by  (1)  and  (2).  This  proves  the  result.  If  (1)  holds  to  any  number  of  terms 
and  and  di  ^ are  0(Xj).  for  any  h,  then  the  above  computation 

holds  for  any  K,  and  the  extension  to  = <»  bolds.  The  extension  to 
infinite  M , N can  similarly  be  justified  provided  that  the  infinite  series 
defining  converge. 

A sequence  of  functions,  n = 1,  ....  A will  be  called  a multipli- 
cative asymptotic  sequence  is  an  asymptotic  sequence,  c6,  = 0(1) 

and  <f>^  <f>^  ^ ^ nmk  k terms,  m,  n = 1,  ...  , A . In  the  case  of  a 

multiplicative  asymptotic  sequence  the  former  result  on  the  multiplication 
of  asymptotic  expansions  can  be  extended  considerably. 

If\<j>^,  n - 1,  , N is  a multiplicative  asymptotic  sequence, 

to  N terms  i = \,  ...  , h, 

and  P {z  ...  , z if)is  a polynomial  in  the  h complex  variables  z ...  , z 
then  F{x)=  P {f^,  ...  , fif)  possesses  an  asymptotic  expansion  2 .4^  9^ 
to  N terms,  and  the  coefficients  A ^ may  be  computed  by  formal  substi- 
tution. 
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To  prove  this  theorem  we  remark  that  in  the  case  of  a multiplicative 
asymptotic  sequence  we  have  and  we  also  have  an 

asymptotic  expansion  of  (f> ^ to  N terms.  By  our  general  theorem,  it 
follows  from  / ~ (f}^  to  N terms  and  g ^ b ^ ^ to  A'  terms  that  fg 

possesses  an  asymptotic  expansion  S to  N terms,  and  the  coef- 

ficients may  be  computed  by  formal  substitution.  The  evaluation  of 
any  polynomial  P{f^,  ...  , can  be  reduced  to  a finite  number  of  oper- 
ations each  of  which  involves  either  a linear  combination,  or  the  multi- 
plication, of  two  asymptotic  expansions.  Each  of  these  operations  pre- 
serves the  asympi  .tic  character  of  the  expansion,  and  in  each  operation 
the  resulting  expansion  may  be  computed  by  formal  substitution.  Hence 
the  theorem. 

The  result  obtained  for  polynomials  can,  under  certain  circumstances, 
be  extended  to  (convergent),  power  series,  and  even  to  asymptotic  power 
series.  For  the  sake  of  simplicity,  we  shall  restrict  ourselves  to  the  case 
of  a single  variable  z;  there  is  a generalization  to  the  case  of  several 
variables. 

Let  n = 1,  ...»  A'  be  a multiplicative  asymptotic  sequence  such 

that  <^,  = o(l),  and  for  some  positive  integer  M.  If 

f(.z)  ~ i c^  z“  to  M terms  as  z ■*  0 in  the  complex  plane,  and 

z = z (x)  ~ 2 a 4>  to  N terms  as  x ■*  x.  in  R, 

then  F {x)  = fiz{x))  possesses  an  asymptotic  expansion  ^ cf>^  to  N 
terms  as  x -*  x^,  and  the  coefficients  may  be  computed  by  formal 
substitution. 

Proof:  From  the  assumptions  it  follows  that  z"  possesses  an  asymptotic 
expansion  2 to  A terms,  and  also  that  Hence  we  can 

apply  the  theorem  in  sec.  1.4  on  the  substituti6n  of  an  asymptotic  expan- 
sion into  an  asymptotic  expansion. 

An  important  particular  case  concerns  functions  f (x)  which  possess 
asymptotic  expansions  of  the  form 

yy 

(3)  /■(x)=c+  2 a^<p^+ o{<f>^.). 

n=  1 

The  theorem  shows  that  [/^(*)]  * also  possesses  an  asymptotic  expan- 
sion of  this  form  provided  that  c ^ 0 and  \ <f>„\  satisfies  the  assumptions 
of  the  theorem.  In  other  words,  asymptotic  expansions  of  the  form  (3)  rfiay 
be  divided.  This  enables  us  to  extend  the  last  theorem  but  one  from 
polynomials  to  rational  functions. 
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'*  = 1»  ••• » /V  is  a multiplicative  asymptotic  sequence,  ^|=o(l), 
[01 1 ^ for  some  M,  f.'^  2.  a .(j>  to  N terms,  i = 1 , , k,  and 

•••  f ^ rational  function  in  the  k complex  variables  Zp  ...  , 

such  that  the  denominator  is  different  from  zero  when  Z|  = ^ =»»*  = Zj^*=0; 
then  F{x)^P{f,  ...  ,4)  possesses  an  asymptotic  expansion  ^ A 
to  N terms,  and  the  coefficients  A ^ may  be  computed  by  formal  substi- 
tution. 

Under  the  same  conditions  we  also  have  an  asymptotic  expansion  for 
g{F  {x))  if  g (4)  is  a function  of  the  complex  variable  which  is  regular 
in  some  neighborhood  of  = P (0,  ...  , 0).  In  this  manner  we  may  justify 
the  asymptotic  expansions  of  expressions  such  as  ex  plPf/",,  ...  , f^)]. 

1.6  Asymptotic  power  series 

The  sequence  of  functions  lx~"|,  n « 0,  1,  2,  ...  or  n ■ 1,  2,  ...  is  a 
multiplicative  asymptotic  sequence  forx-»  00  in  any  region  of  the  complex 
plane  which  does  not  include  the  origin.  This  sequence  satisfies  all  the 
conditions  imposed  upon  asymptotic  sequences  in  the  two  preceding 
sections,  except  that  in  some  of  the  theorems  of  sec.  1.5,  n 0 must  be 
excluded.  Besides,  this  system  has  some  special  properties. 

The  asymptotic  expansion 

a a 

(1/  fM-^a  + — i+— !-+•••  to /V  terms  as  X -►  00 

is  an  asymptotic  power  series.  From  the  results  of  sections  1.5  and  1.6 
it  follows  tliat  an  asymptotic  power  series  expansion  may  be  multiplied 
by  a constant,  and  that  two  such  expansions  may  be  added  or  multiplied, 
and  also  divided  provided  that  ^ 0 inthe  expansion  in  the  denominator. 
Asymptotic  power  series  may  be  substituted  in  finite  linear  combinations, 
in  polynomials,  in  rational  functions  provided  that  the  denominator  does 
not  vanish  as  * -♦  00,  and  in  asymptotic  or  convergent  power  series 
S c z",  z -*  0,  provided  that  in  the  expansion  (1)  of  z = f(x)  we  have 
Oq  = 0,  Substitution  of  (l)inother  types  of  convergent  or  asymptotic  series 
is  valid  under  the  conditions  set  out  in  sec.  1.5.  In  all  these  cases  the 
coefficients  of  the  new  expansion  are  obtained  by  formal  substitution 
and  a rearrangement  of  terms.  An  asymptotic  power  series  expansion  (1) 
which  is  valid  uniformly  in  a parameter  may  be  integrated  with  respect 
to  this  parameter.  Lastly,  if  (1)  holds,  then  ff*)  ~ °o  ~ inte- 

grable,  and 


t 
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a,  a a 

~ — H = H •••  to  N — 2 terms  as  x 

X 2x^  3*® 


A simple  corollary  of  this  last  result  is  the  following  theorem  on  differ- 
entiation.//  fix)  in  (1)  is  differentiable  and  if  f'i^c)  possesses  an  asymptotic 
power  series  expansion,  then 

Q 2 o 3 o 

(3)  f'ix)  ~ i — ' ^ — to  A — 1 terms  as  ac  -►  <». 

X X x 


In  the  case  of  analytic  functions  a more  definite  statement  can  be 
made  in  that  it  is  not  necessary  to  assume  that  f'ix)  possesses  an 
asymptotic  power  series  expansion.  Let  R be  the  region 

1*1  > a,  a < arg  * < /3, 

let  a I > o,  a < a,  < /3 , < /3,  and  let  R , be  the  region 
1*1  > a,,  a,  < arg*  < /3,. 

If  fix)  is  regular  in  R and  (1)  holds  uniformly  in  arg  x as  x -*  oo  in  R,  then 
(3)  holds  uniformly  in  arg  x as  x oo  in  R ^ . Ibe  proof  of  this  theorem 
follows  from  Cauchy’s  integral  formula  for  the  derivative. 


(4)  f 


2^^*  Jc 


1 

(*  - z) 


For  given  R,  there  exists  an  f > 0 so  that  for  each  * in/?,,  the 
circle  with  center  * and  radius  f |*|  is  in  R,  and  we  may  take  this  circle 
as  the  contour  of  integration  in  (4).  Along  the  circle,  z = * + exe  **,  and 
0 < t < 2n-,  so  that  (4)  becomes 


(5)  f'(x)=- 


-f 

Jo 


e"*7l*(l  + €e‘*)]dt. 
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Now,  e “/[(*  4-  c e **)]  possesses  an  asymptotic  power  series  expansion 
uniformly  in  t,  and  this  may  be  integrated  with  respect  to  t,  showing  that 
f (*)  possesses  an  asymptotic  power  series  expansion  which  turns  out  to 
be  (3). 

Asymptotic  power  series  expansions  are  usually  valid  in  sectorial 
regions,  and  analytic  functions  possess  different  asymptotic  expansions 
in  different  sectors  {Stokes’  phenomenon).  That  something  like  this  must 
happen,  except  in  the  case  of  an  analytic  function  which  is  regular  at 
infinity,  follows  from  the  following  theorem. 

If  f{x)  is  single-valued  and  regular  when  |x|  > a,  and  (1)  holds  for  all 
values  of  arg  x,  then  the  power  series  in  (1)  converges  for  sufficiently 
large  values  of  |x|,  and  its  sum  if  f{x\  To  prove  this,  we  set  x = 1/^ 
and  g(0)  = a^,  g = f{\/^),  0 < |^|  < |a|“'.  Then  g(^)  is  a single- 
valued continuous  function  in  Kl  < |a|  *,  and  is  regular  except  possibly 
at  ^=0.  However,  at  0,  g has  certainly  no  pole,  nor  an  essential 
singularity,  since  it  is  bounded  in  any  neighborhood  of  ^ = 0.  Thus, 
g(^)  is  regular  at  ^ = 0 and  possesses  a Maclaurin  expansion.  From  the 
uniqueness  theorem  on  asymptotic  expansions  it  follows  that  (1),  with 
X - 1/^,  must  be  the  Maclaurin  expansion. 

1.7.  Summation  of  asymptotic  series 

It  has  been  pointed  out  in  sec.  1.3  that  an  asymptotic  sequence 
10^1  determines  an  equivalence  relation  between  functions  defined  in  R. 
Two  functions  defined  in  R are  asymptotically  equal  if  their  difference  is 
o(0^)  for  all  n.  Asymptotically  equal  functions  possess  identical  asym- 
ptotic expemsions,  and  given  an  asymptotic  expansion  /■  ~ 2 0^,  we 

may  define  the  class  of  all  functions  which  are  asymptotically  equal  to  f 
as  the  sum  of  the  asymptotic  series  2 0„. 

We  shall  conclude  this  chapter  by  proving  that  every  asymptotic 
series  possesses  a sum.  Results  of  this  nature  have  been  proved  for 
asymptotic  power  series  by  Borel  and  Carleman  (1926),  for  series  domin- 
ated by  an  asymptotic  series  of  powers  by  van  der  Corput  (1954  b),  and 
for  asymptotic  series  of  analytic  functions  by  Carleman  (1926).  The 
proof  given  below  is  adaptation  of  van  der  Corput’s  propf. 

An  asymptotic  series  is  a formal  finite  or  infinite  series  2 0^(x) 

where  10„i  is  an  asymptotic  sequence  and  the  a^  are  constants.  Since 
any  subsequence  of  an  asymptotic  sequence  is  also  such  a sequence, 
we  may  assume  that  a^/:  0 for  each  n.  The  asymptotic  sum  of  2 0„  is 

a class  of  asymptotically  equal  functions,  and  we  shall  demonstrate 
the  existence  of  the  asymptotic  sum  by  constructing  a member  of  this 
class.  If  2 0^  is  a finite  asymptotic  series,  then  the  sum 
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a, 

in  the  ordinary  sense  may  be  taken  as  a representative  of  the  asymptotic 
sum.  It  is  sufficient,  then,  to  give  the  proof  for  an  infinite  asymptotic 
series  2 a 0 in  which  a ^0  for  each  n. 

Let  be  a neighborhood  of  x^,  and  for  each  n 1,  2,  ...  let  U ^ be  a 
neighborhood  of*,  such  that  the  closure  of  V is  in  {/  and 

D FI  n ■ 

la„+, 

for  all  X common  U ^ and  R:  such  a neighborhood  exists  since 

For  each  n let  u (*)be  a continuous  function  of  * such  that  0 < u (*)  < 1 

* FI  * FI  “ 

in  R,  ti  (*)  = 0 when  * is  outside  U and  u (*)  = 1 when  * is  in  f/  . , : 

•FI  Fl'FI  Fl“1 

such  a function  exists  since  the  closure  o{  U a,  is  contained  in  f/  . Then 

when  * is  in  U for  this  inequality  holds  by  the  construction  of  the  U*s 
if  X is  in  {/  . , and  the  left-hand  side  vanishes  when  x is  outside  U . . 

n+p’  n+p 

Let 


(2)  /■(*)=  2 a^fi^ix)  <f)^(x). 

n=  1 

The  series  converges  for  all  * by  (1),  and  defines  a function  f (x)  in  R. 
(Actually,  the  series  terminates  except  for  those  * which  are  in  all  the 
U .)  To  show  that  f ~ 2 a di  as  * -*  x„,  fix  /V  and  let*  be  in  the 

JX  * ' FI  * FI  IF 

common  part  of  and  R.  Then  p„(*)  = 1 for  n = 1,  ...  , N,  and  by  (1) 


+«■ 


= 21«,v+,  9a.+,1  = 

Thus  it  is  seen  that  2 is  an  asymptotic  expansion  to  any  nunrber 

of  terms  of  f defined  by  (2).  The  asymptotic  sum  of  2 is  the  class 

of  all  functions  asymptotically  equal  to  f. 
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The  U ^ may  be  constructed  in  such  a manner  that  x ^ is  the  only  point 
common  to  all  the  U ^ in  which  case  the  series  in  (2)  terminates  for  all 
X / If  all  the  are  continuous  in  R,  also  f will  be  continuous  in  R. 
If  X is  a real  variable,  or  a point  in  n-dimensional  Euclidean  space,  the 
p^(*)  may  he  chosen  as  infinitely  differentiable  functions,  and  if  all  the 
<f) ^ are  k times  continuously  differentiable  (k  < oo),  then  fix)  will  also  be 
k times  continuously  differentiaule.  Carleman  has  proved  that  for  certain 
analytic  functions  <f>^  of  a complex  variable  *,  the  asymptotic  sum  con- 
tains a function  which  is  an  analytic  function  of  *. 

In  general  there  is  no  v'ay  of  ascribing  a unique  asymptotic  sirni  to 
an  asymptotic  series,  but  under  rather  special  circumstances  it  may 
happen  that  under  more  precise  assumptions  on  the  coefficients  of  the 
asymptotic  series,  and  under  certain  restrictions  on  the  functions  f (x), 
a unique  sum  may  be  obtained;  and  frequently  in  such  cases  the  asym- 
ptotic series,  though  divergent,  is  in  some  sense  summable  t.o  its  asym- 
ptotic sum.  Such  theorems  for  asymptotic  power  series  summed  by 
analytic  functions  regular  in  some  sectorial  region  were  obtained  by 
Watson  (1912  a)  and  Nevanlinna  (1916). 
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CHAPTER  n 

INTEGRALS 


There  are  several  methods  for  obtaining  asymptotic  expansions  of 
functions  defined  by  definite  integrals.  Copson  (1946)  gives  a survey  of 
these;  and  further  material  is  contained  in  van  der  Corput’s  l,ectures  and 
in  the  references  given  at  the  end  of  this  chapter. 

2.1.  Integration  by  parts 

Asymptotic  expansions  may  frequently  be  obtained  by  repeated  inte- 
grations by  parts.  As  an  example,  let  us  consider  the  function  f (x)  defined 
for  —n  < arg  * < tt  by  the  integral 


(1)  A*) 


dt. 


Integrating  by  peirts  repeatedly, 


(2)  fix)  = 1 


(1  + xt)^ 


dt 


l-x+2**  / ^ 

Jo  (1  + xt)^ 


dt 


=:  • • ^ 


t + (m+  1)!:.  , 

n~  0 ^ (1  + Xf) 

0 


+ 2 


dt 
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The  last  integral  may  be  proved  to  be  O (1)  as  ar  -*  0 in  A > —tt/2,  so 
that  we  have  obtained  a new  derivation  of  Euler’s  asymptotic  expansion 
discussed  in  the  Introduction. 

The  field  of  application  of  this  method  is  somewhat  limited,  and  it  is 
not  at  all  easy  to  formulate  precise  theorems  of  sufficient  generality.  In 
what  follows  we  shall  describe  some  results  which  seem  to  be  basic. 

For  any  function  f{t),  let  denote  the  m-th  derivative,  and  the 
m-th  repeated  integral,  so  that 

(3)  /o  - A /.  m = 1,  2,  ... 

(4)  m = 1,  2,  ...  . 
at 

Note  that  contains  rt  constants  (one  from  each  integration)  which  we 
suppose  to  have  been  chosen  in  some  suitable  manner.  The  formula 

(5)  f/  gU)h(t)dt^ 


where 

(6)  s„«(-l)"[g„(/3)A.„.,(/3)-  g„(a)A.„_,(a)] 

(7)  /?„  = (-l)"//g„(OA_„(Odt, 

is  obtained  by  repeated  integrations  by  parts.  lf'(a,  /S)  is  a finite  interval, 
(5)  is  valid  provided  that  g is  A times  continuously  differentiable  and  A 
is  integrable;  if  (a,  /3)  is  an  infinite  interval  then  all  the  integrals  in- 
volved. and  also  the  limits  of  g (t)  A _ (t)  as  t -*  a,  B,  must  be  assumed 

to  exist. 

If  g is  IV  + 1 times  continuously  differentiable,  a further  integration 
by  parts  shows  that 

(8)  R jif  ~ ^ N + > 

and  in  certain  cases  it  is  possible  to  use  this  relation  to  compare  the 
“remainder”,  with  the  first  “neglected  term”  s 
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If  g and  h are  real,  and  g.^  h_^  and  have  constant  and 

equal  signs  for  a ^ t £ /S,  then  R has  the  same  sign  as  and  |/i^| 
< The  proof  follows  from  (8)  on  noting  that  in  this  case  and 

have  opposite  signs,  and  hence  /?^  and  s^  must  have  the  same 
sign,  and 

If  g is  real,  is  an  increasing  function  of  t,  and  g^, 

have  constant  and  equal  signs  for  a < t <f>,  or  else  if  g is  real,  _,| 
is  a decreasing  function  ‘of  t,  and  g have  constant  and  opposite 
signs  for  a < t < then  lf?y|  < 2|s  „|. 

We  shall  prove  this  result  when  is  an  increasing  function  of  t 

and  > 0,  > 0.  Then 

(/3)|  (g^(/9)-g^(a)) 
<|A-^.,(/3)g//3)-A_^.,(a)g^(a)| 

and  hence  |^yy+,|  < From  (8)  we  then  have  the  desired  result.  If 

gy  <0  and  gyy+i  < 0,  replace  g by  — g.  The  result  for  decreasing 
follows  on  replacing  x by  — 

As  an  application  of  these  results,  let  us  consider /'(x)  as  defined  by 
(1).  If  X > 0,  put  git)  * (1  + xf)“',  h it)  - e~\  h_^  it)  » (-  1)*  e“*.  In  this 
case  g ^ > 0 for  all  f > 0,  and  hence  0 < (—1)"  R ^ < (~1)"  s,  • If* 

is  complex  our  results  do  not  apply.  However,  if  in  (1)  we  replace  t by 
l/x  and  accordingly  set  g =*  (1  + t)~’,  h - x“'  exp(-t/x),  and  then  letx 
become  complex,  with  Re  x > 0,  then  g^  and  g^  are  of  constant 
and  opposite  sign,  and  is  a decreasing  function  of  C,  for  t >0, 

and  hence  I R | < 2|5  |.  (Actually,  in  this  case  it  is  easy  to  prove  from 
(7)  that  |f?_|  < ls,|.) 

Let  us  suppose  now  that  the  integrand  in  (5),  and  hence  also  s^,  Ry, 
depend  on  a variable  x.  If  ts^i  is  an  asymptotic  sequence,  and  if  in 
addition  we  are  able  to  prove  R^  « Ois  ) by  one  of  the  above  results, 
or  in  some  other  manner,  then  (5)  provides  an  asymptotic  expansion  of 
the  integral  to  N terms.  For  instance,  in  the  case  of  (2),  tx"i  is  an 
asymptotic  sequence  for  x -*  0;  we  have  proved  \Ry\  ■^2|Sjy|  for  any  N 
and  Re  x >0,  and  hence  (2)  is  an  asymptotic  expansion  of  fix)  as  defined 
in  S^,  A > 0.  [Actually,  we  proved  in  the  Introduction,  and  could  prove 
fron<  the  last  integral  in  (2),  that  the  asymptotic  expansion  holds  in  the 
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more  extended  region  S^,  A > -17/2.] 

An  asymptotic  sequence  often  occurs  if  h (t)  = k{xt).  Denoting  by 
(u)  the  m-th  repeated  integral  of  k (u)  with  respect  to  u,  we  obtain 
from  (5), 

<9)  g(t)k{xt)dl=''±'  [g^(/3)A.^_,(/3x) 

fi  * 0 

If  the  k_^  (u)  are  bounded,  and  the  g^(/3)  i^x)  — g„(a) 

are  bounded  away  from  zero,  and  if  R ^ can  be  estimated  as  above,  then 
(9)  is  an  asymptotic  expansion  as  x ->  bo,  the  region  of  x being  determined 
by  the  estimate  of  R The  apparently  more  general  case  A U)  = A [x  9f)(t)] 
can  be  reduced  to  the  former  case  by  breaking  up  (a,  /S)  into  sub-inter- 
vals in  which  <p(t)  is  monotonic,  and  introducing  <^(t)  as  a new  variable 
in  each  of  these  sub-intervals.  In  applying  one  or  the  other  of  the  above 
criteria  for  the  estimate  of  R^,  we  then  need  information  about  derivatives 
of  the  form 

r git) ' 

d<t>‘  L _ * 

If  g and  have  derivatives  of  constant  or  alternating  signs,  this  infor- 
mation can  be  obtained  without  explicit  computation  from  results  on 
absolutely  and  completely  monotonic  functions,  see,  for  instance, 
Widder  (1941)  Chapter  IV.  General  theorems  of  this  nature  have  been 
obtained  by  van  der  Corput  and  Franklin  (1951).  The  most  important 
application  of  these  methods  is  to  integrals  of  the  form 

J-'’g(r)e'*<‘>dr.  f"  g(t)^f^lxk(t)]dt. 

a o 

More  general  results  involving  functions  of  the  form  g(t)=  (t  — a)~^g  ^(t) 
with  g,(t)  possessing  continuous  derivatives  were  obtained  by  van  der 
Corput  (1934). 

2J2.  Laplace  integrals 

Integrals  of  the  form 

(1)  /(x)  = dt  = 52  1^1 

0 
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are  called  Laplace  integrals , Such  integrals  occur  in  the  solution  of 
differential  equations  by  definite  integrals,  and  in  many  other  problems. 
The  infinite  integral  in  (1)  will  be  interpreted  as  the  limit  of  as 
T -*  oo,  and  it  will  always  be  assumed  that  <f>{t)  is  integrable  over  any 
interval  0 < t < T,  T <oo.  A function  (f>  will  be  said  to  belong  to  L (ac^) 
if  the  inte^al  in  (1)  exists,  in  the  sense  mentioned,  for  x = x^.  It  is 
known  (Widder,  1941,  Chapter  II)  that  for  a function  <f>  in  L (*„),  010} 
exists,  and  represents  €ui  analytic  function  of  x,  in  the  half-plane 
Re  (x  — Xq)  > 0.  In  particular,  if  0(t)  is  integrable  over  any  interval 
0 < f < r,  r < oo,  and  0(t)  = 0(e“*)  for  some  constant  a,  as  t -»  <»,  then 
O|0iexists  as  an  (absolutely  convergent)  infinite  integral,  and  repre- 
sents an  analytic  (unction  of  x,  in  the  half-plane  Re  (x  - a)  > 0. 

Under  certain  circumstances,  the  asymptotic  behavior  of  f (x)  as  x -♦  <» 
may  be  investigated  by  integrations  by  parts.  If  0(t)  is  N times  con- 
tinuously differentiable  for  0 < t < a and  belongs  to  L{x^)  for  some  x^ 
then 

(2)  /(x)~5:  0<"kO)x-"~' 

to  N terms,  uniformly  in  arg  x,  as  x -*  oo  in  S^,  A > 0.  To  prove  this,  let 
Re  (x  *“  *o)  > Of 

(3)  fix)  = f“e  ”**  <f)it)  dt  ■¥  e “**  0 it)  dt. 

The  second  integral  exists  and  may  be  integrated  by  parts  to  give 


1 


X - x„ 


“(x- 


0 it)  dt. 


rhere 


0(t)“/J  « ® 0(u)du 

is  a bounded  function,  say  |0|  < A tor  t > a,  so  that 


l*-*ol  Ip - Pol 


= 0(e"^“) 
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as  * - p + tCT  -♦  oo  in  Sa,  a >0.  In  the  first  integral  use  2.1(5)  with 
g = = (-*)""  e"« 

s ^ ’(0)  - ^\a)  e **«]  x‘^~' 

= 9S<"^(0)*"»“’  + Oie’P") 

and 

^AT  = *''^  4>^Kt)  e--*  dt  . 

Here  being  continuous,  is  bounded,  say  |0  < B for  0 ^ t < a, 

and 

|/?^|<fi|a:|  ^ p ’ < B |ie|  cosec  A = 

uniformly  in  arg  x,  as  x -*  oo  in  S^,  A > 0.  The  proof  is  completed  by 
noting  that  0(c~^“)  •=  o(*~^),  uniformly  in  arg  x,  as  x -*  oo  in  A > 0, 
A considerable  extension  of  the  last  theorem  may  be  based  on  the 
following  LEMMA.  Let  <^U)  and  be  in  L (x^)  for  some  x^,  xpit)  > 0, 

/■  = Q {0|,  g Q \ip\.  If  e"P  gip)  -*  oo  as  p ■*  BO  for  each  a > 0,  and  if 
(fit)  - o(i/f(f))  05  t -*  0,  then  f (x)  ■=  o(g(p)),  uniformly  in  arg  x,  as 
X - p + i a -*  oo  in  S^,  A > 0. 

Proof:  Given  f > 0,  there  exists  an  o > 0 so  that  |<^|  ^ exp  for  0 <t  ^o. 
^itb  this  o,  decompose  Q i^i  as  in  (3).  As  in  the  proof  of  the  previous 
theorem,  the  second  integral  is  0 (e~“0  as  p -♦  oo,  and 


|/°  e~**  (pit)  <ft|  <e  J°  e~P*xpit)  dt  = f gip). 
0 ® 


Thus 


\fM\ 

g(p) 


<e  + 0 


and  this  is  < 2f  for  sufficiently  large  p.  The  uniformity  in  arg  x follows 
from  the  remark  that  l*|/p  < cosec  A in  S^. 

From  this  lemma  the  followingthecrem  can  be  deduced.  For  ra  = 1, ... , A, 
let  tp  it)  be  in  L i?c^)  for  some  x^,,  xp ^it)  > 0 for  t > 0,  and  g^  = Q.\x//\. 
If  1«A„1  is  cm  asymptotic  sequence  for  f -♦  0,  and  c“^  OO  as  p oo 

for  each  o > 0 and  each  n,  then  lg^(p)l  is  an  asymptotic  sequence,  for 
p ^ +00  ; and  if  under  these  circumstances  (pit)  is  in  L (x^)  and 
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0 (f)  ~ 2 a ^ (t)  to  JV  terms  as  t 0, 

then 

/(p)  ~ 2 8 ^ terms  as  p -»  +<». 

If  in  addition  for  each  n = 1,  , /V,  g„(p)/g„(*)  is  bounded  in  for 

sufficiently  large  \x\,  then  also  |g^(a:)j  is  an  asymptotic  sequence,  and 
fix)  ~ 2 g^ix)  to  N terms,  uniformly  in  arg  x,  as  x -*  oo  in  S^,  A > 0. 
Proof:  From  the  lemma  we  have  g^  + ,{p)  = o(g^(p)),  and  hence  {g„(p)l 
is  an  asymptotic  sequence  for  p -»  oo  To  prove  the  asymptotic  expansion 
for  p -*  oo,  replace  <f)  by 

n*  I 

and  0 by  in  the  lemma.  When  x->oo  in  S^,  we  have  from  the  lemma  and 
the  additional  assumption  on  g„0c)  that 

g„+,  (*)  - o (g„(p))  oigjx))  - oig^ix)). 


and  hence  lg„(*)l  is  an  asymptotic  sequence  for  * -»  oo  in  S^.  The  proof 
of  the  asymptotic  expansion  is  the  same  as  in  the  previous  case. 

The  most  impo.'tant  particular  case  of  our  general  theorem  is 

, 0 <A, 


All  conditions  of  the  theorem  are  satisfied:  in  particular,  for 

g„(x)  = r(A„)x'^» 


we 


have 


g ip)  / 1*1  \ ^ 

< ( — ) < (cosec  A)  " 

\ P / 


for  X in  A > 0.  We  then  obtain  the  following  theorem  on  asymptotic 
series  of  powers. 
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Let  0 < A,  < Aj  < •••  . //  ^ (/)  is  in  L (x^^)  for  some  and 
, ^ ^ "1 

o~Sat"  to  N terms  as  t -♦  0, 

fl  " 

then 

-A  . , . 

/ ~ £ I (A^)  X " to  N terms,  uniformly  in  arg  x, 

as  X -*  oo  in  A > 0. 


Other  notable  examples  of  asymptotic  sequences  to  which  the  general 
tueorem  applies  are 

(4)  =(l-e-')"-\ 


x(x  + 1)  ...  (x  + n - 1) 


(5)  0 =(c*-l) 


n“1 


U-D!  ■ 

X (x  — 1)  ...  (x  — a + 1) 


(6)  «A 


/ t Y'" 

— I 2 sinh  — - ) 

V 2;  . 


(2a  - 2)! 


(x  - a + l)(x  — a + 2)  ...  (x  + a — 1) 


In  the  case  of  asymptotic  power  series  or  in  the  case  of  (4),  N may  be 
finite  or  infinite,  in  the  cases  of  (5)  and  (6)  A'  must  be  finite. 

The  result  obtained  by  integrations  by  parts  is  a particular  case  of 
of  the  above  theorem  for  asymptotic  series  of  powers.  If  f^{t)  is  A times 
continuously  differentiable  for  G < t < a,  then  it  can  be  proved  by  the 
mean  value  theorem  of  differential  calculus  that  9!>(t)  ~ S t"/a! 

to  A terms  as  t -»  0,  and  then  (2)  follows  from  the  theorem  on  asymptotic 
power  series. 

In  many  cases  it  is  possible  to  extend  the  region  of  validity  of  the 
asymptotic  expansion  to  an  5^  with  A<  0.  If  (^{t)  is  an  analytic  function 
of  t which  is  regular  in  and  is  O (e  “*)  for  some  constant  a,  as  t -♦  00 
in  Sg,  then  it  is  permissible  to  rotate  the  path  of  integration  in  (1)  to  any 
ray  in  S^,  and  by  this  means  fix)  may  be  continued  analytically  to  some 
region  which  contains  the  sector  — w + 0 < arg(x  — a)  < jt  — 0 [see,  for 
instance,  Doetsch  (1950,  p.  362 fif.)] . K 1 ’A „l  has  suitable  properties  along 
each  ray,  and  ~ ^ ® n ^ ^ in  S^,  then  ~ S a ^ g ^ as  x -♦  00  in 

S^where  A>0  - m/2.  Die  shall  formulate  a precise  theorem  for  tisymptotic 
series  of  powers. 
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If  <^(t)  is  a regular  function  of  t in  S g,  <f>  = uniformly  in  arg  t, 

for  some  a,  as  t -*  oo  in  S g, 

A.^-1 

<f>  2.  a ^ t to  N terms,  uniformly  in  arg  t,  aa  t -*  0 in  Sg, 

where 

0 < A,  < Ajj  < ...  < A^, 

then  f(x)  exists  at  least  in  the  sector 


-n  + 6 < arg(*  — a)  < tt  — 0, 

and 

-A 

f a ^ r(A^)  * " to  N terms,  uniformly  in  arg  x, 

as  X -»  oo  in  5^,  A > 0 — n/’l. 

The  particular  case  of  this  result  in  which  the  A^  form  an  arithmetic 
progression,  and  ^ is  represented  by  a convergent  infinite  series 


for  sufficiently  small  |t|  in  5^  is  known  as  Watson’s  lemma',  it  is  suffic- 
ient for  many  applications. 

Throughout  this  section  we  investigated  the  behavior  of  (1)  for  large  x. 
Similar  methods  may  be  used  for  the  investigation  of  fix)  as  x -*  x^.  For 
the  basic  lemma  see  Erdelyi  (1947),  and  for  some  of  the  most  important 
results  see  Doetsch  (1950)  Chapter  13  and  Widder  (1941)  Chapter  V. 


2.3..  Critical  points 

We  have  seen  in  the  last  two  sections  that  under  certain  circum- 
stances the  asymptotic  behavior  of  integrals  is  determined  by  the  beha- 
vior of  the  integrands  at  certain  distinguished  points,  the  end-points  of 
the  interval  in  the  cases  considered  in  the  preceding  sections.  Such 
distinguished  points  have  been  called  critical  points  by  van  der  Corput 
(1948).  There  is  no  general  theory  of  critical  points  but  a few  types  of 
such  points,  and  the  methods  adapted  to  deal  with  them  may  be  described 
€is  follows. 
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First  let  us  consider  an  integral  of  the  form 

(1) 

o. 

where  ac  is  a large  positive  parameter  and  h (t)  is  real.  If  h (t)  lias  a 
maximum  at  t = r,  and  h(t)  <k  (r)  when  t ^ r,  then  for  large  *,  the  modu- 
lus of  the  integrand  will  have  a sharp  maximum  at  a point  vcrv  near  r,  and 
most  of  the  contribution  to  the  integral  will  arise  from  f.n*  immediate 
vicinity  of  this  maximum.  The  integral  can  be  evaluated  approximately  by 
expanding  both  g and  h in  the  neighborhood  of  t = r.  This  is  the  central 
idea  of  Laplace’s  method  (sec.  2.4).  l^e  have  encountered  such  a case  in 
sec.  2.2  where  h{t)  = —t,  0 < t < oo,  and  h (t)  has  a maximum  at  t = 0. 
Accordingly,  we  evaluated  Laplace  integrals  asymptotically  by  expand- 
ing g (t)  for  small  values  oft. 

If  X and  h (t)  are  complex,  and  g (t)  and  h {t)  are  analytic  functions  of  t, 
then  it  is  often  possible  to  deform  the  path  of  integration  so  that  it 
passes  through  one  or  several  points  at  which  h ' {t)  = 0,  If  r is  such  a 
point,  it  is  a critical  point;  it  is  possible  to  determine  that  part  of  the 
path  of  integration  which  passes  through  r in  such  a manner  that 
X [A  (t)  — A(r)]  is  real  along  the  path,  and  the  integral  can  then  be  evalu- 
ated by  an  adaptation  of  Laplace’s  method.  This  is  Riemann’s  metAod  of 
steepest  descents  (sec.  2.5). 

IVext  let  us  turn  to  integrals  of  the  form 

(2)  g(t)e“*<‘>dt 

a 

where  we  again  assume  that  * is  a large  positive  parameter,  and  h it) 
is  a real  function.  In  general,  the  rapid  oscillations  of  exp[ixA(t)]  will 
tend  to  cancel  large  contributions  to  the  integral,  but  this  cancellation 
will  not  occur  at  the  end-points,  or  at  the  stationary  points  of  k it).  If 
hit)  has  no  stationary  points  in  the  interval  a < t < ^,  integration  by 
parts  (sec.  2.1)  will  in  general  give  a good  approximation.  Stokes’ 
method  of  stationary  phase  (sec.  2.9)  appraises  the  contribution  of  a 
stationary  point,  r,  to  the  integral  by  expanding  g and  h in  the  neighbor- 
hood of  this  point. 
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The  method  of  stationary  phase  has  been  extended  by  van  der  Corput 
(1936)  to  integrals  of  the  form  (1)  where  xh  it)  may  be  complex  (instead 
of  being  imaginary,  as  in  the  method  of  stationary  phase).  According 
to  van  der  Corput,  the  critical  points  in  this  case  are  those  points,  r, 
at  which  A ' (t)  [A"(t)]“^  is  real,  while  the  imaginary  part  of  this 
function  changes  its  sign  when  t passes  through  r.  (In  the  case  of  (2), 
h = ik,  k is  real,  and  the  only  points  at  which 

A'  (*")■’<  = (/x)^  k'  Vc”  )~^ 

is  real  are  the  stationary  points  of  A.) 

2.4.  Laplace's  method 

In  the  integral 

(1)  /(x)./^g(Oe**<‘>* 

a 

let  A(t)  be  a real  function  of  the  real  variable  t,  while  g (t)  may  be  real 
or  complex,  and  let  x be  a large  positive  variable.  According  to  Laplace, 
the  major  contribution  to  the  value  of  the  integral  arises  from  the  immed- 
iate vicinity  of  those  points  of  the  interval  a <,  t /S  at  which  A (t) 
assumes  its  Icmgest  value.  If  A if)  has  a finite  number  of  maxima,  we  may 
break  up  the  integral  in  a finite  number  of  integrals  so  that  in  each 
integral  A (f)  reaches  its  maximiun  at  one  of  the  end-points  and  at  no 
other  point.  Accordingly,  we  shall  assume  that  A (t)  in  (1)  reaches  its 
maximum  at  t = a,  and  that  h{f)  < fi  (a)  for  a < t < /3. 

Assuming  g continuous  and  A twice  continuously  differentiable, 
A * (a)  = 0,  A"  (a)  < 0,  Laplace  introduced  a new  variable  u by  the  sub- 
stitution A (a)  - A(t)  = u*.  A ' it)  will  be  negative  ina<£  < a + rj  for 
some  sufficiently  small  77.  As  x -►  00, 

h'(t)  lexpx[A(a)-u*]|  du 

where  (/  = [A  (a)  — A (a  + >0.  Since  only  the  neighborhood  of  u = 0 

matters,  we  may  replace  g(£)  approximately  by  g{a\  and  u/h'{t)  by 
— [— 2A"(a)]~^,  which  is  the  limit  of  u/h'{t)  as  t -*  a,  and  obtain 


f(x)'^  g(t)  e‘''^‘’d£-  - r 

0 
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By  the  same  argument,  we  may  extend  the  integration  to  u 
obtain  Laplace’s  result 


(2) 


/(*)  ~ g(a) 


K 


oo  and  finally 


X -»  oo. 


Later,  Burkhardt  (1914)  and  Perron  (1917)  showed  that  the  same  result 
can  be  proved  by  expanding  g and  h in  the  neighborhood  of  a.  Copson 
(1946)  reproduces  a simple  proof  of  Polya  and  Szego,  and  \(idder  (1941, 
Chapter  VII)  gives  a more  sophisticated  proof  under  more  general  con- 
ditions. Further  extensions  of  Laplace’s  formula  were  obtained  by  Bsu 
(1949  a,  b;  1951  a,  b),  Levi  (1946)  and  Rooney  (1953).  Laplace’s  method 
has  been  applied  to  integrals  depending  on  two  large  variables  by  Fulks 
(1951)  and  Thomsen  (1954),  and  to  double  and  multiple  integrals  by  Hsu 
(1948  a,  b;  1951  c)  and  Rooney  (1953). 

The  following  extension  of  Laplace’s  result  will  be  derived  from  our 
discussion  of  Laplace  integrals.  Let  g and  h be  functions  on  the  interval 
{a.  /3)  for  which  the  integral  (1)  exists  for  each  sufficiently  large  positive 
X,  let  h be  real,  continuous  at  t = cl,  continuously  differentiable  for 
a<t<a  + Tj,  T}>0,  and  such  that  h'  <0  for  a < t < a + rf,  h it)  < h (a)  - t, 
f > 0,  for  a + 77  < t < /3;  suppose  that  h ' (t)  — a(t  - a)^~'  and  git)  ~ 

b it  — a)^”'  as  t -»  a,  A > 0,  1/  > 0:  then 


(3)  fix)  = j\it) 

OL 


y V OX  y 


X 00. 

V(e  first  note  that 


<4)  iJ^  g (t)  < exp{*[A  (a)  - f]|/^  |g(t)|  dt 

a + 77  a*  r, 


X ^ 


00 , 


In  the  interval  (a,  a + t/)  we  introduce  a new’  variable  u = hia)  - hit), 
set  V = h (a)  - A (o  + 77)  > 0,  A (u)  = - g it)/h ' (t)  and  obtain 

(5)  g(t)  A(i/)e‘”du. 
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u = h{a)-h{t)  = - X,*  h '(r)  aV 


and  hence 


Also, 


so  that 


(6)  A:(u) 


~ — (t  - af 


i-T" 


as  t -»  a 


as  u -►  0. 


as  ? -♦  a 


as  u -»  0. 


By  the  results  of  the  preceding  section  on  the  asymptotic  behavior  of 
Laplace  integrals  it  follows  from  (5)  and  (6)  that 


as  * -►  oo. 


and  (4)  and  (7)  prove  (3).  Moreover,  both  (4)  and  (7),  and  hence  also  (3), 
remain  true  if  x is  a complex  variable  and  x -►  oo  in  S^,  A > 0. 

A further  extension  of  (3)  leads  to  an  asymptotic  expansion  of  fix). 
In  the  following  formulas  n - 0,  1,  ...  , A'  — 1.  If 


(8)  -h'it)'^  g(t)~  2 6„(£-ay 

to  N terms  as  £ -►  a,  then  there  is  an  expansion 


\+n““! 


(9) 

h 'it) 


Sc  (£  — to  A'  terms  as  £ -*  a 


and  the  c may  be  computed  by  formal  division.  Also 


u = — X * h '(r)  dr  ~ / — — (£  — to  N tenr.s  as  £ -»  a. 

a 1/  -h  n 
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From  this  last  expansion  it  may  be  shown  that  t — a possesses  anasymptotic 
power  series  expansion  in  powers  of  this  asymptotic  power  series 

may  be  substituted  in  (9),  and  leads  to  an  asymptotic  expansion  of  the 
form 

(11)  k (u)  1,  y u to  N terms  as  u ->  eo  . 

Kith  (11)  instead  of  (6),  an  asymptotic  expansion  to  N terms  is 
obtained  for  (5);  (4)  may  be  strengthened  to  o and 

we  have 

(12)  /(x)~c*'‘<“>  S r ^ 

to  N terms  as  * -♦  oo  in  S^,  A > 0.  The  coefficients  y^  ma^  be  computed  by 
formal  substitution  according  to  tL  ' scheme  described  above. 

There  is  an  alternative  procedure  for  the  computation  of  the  y^  which 
avoids  the  necessity  for  inverting  the  asymptotic  series  (10)  to  obtain 
the  expansion  of  t •-  a in  powers  of  u From  (10), 

A (t)  = A (a)  - it-ay  + h.  it) 

V 

where 

TV- 1 - 

hAt)=-  y — + 

^ v + n 

n = 1 

Ke  now  write 

jJ'%(Oe‘''^‘>dt  = e*'*<‘^)  //'''^/(Oexp  * («  - a)*' J dt, 

expand 

lU)  = git)  exp[*A  ,.(t)] 

formally  in  powers  of  it  — a),  and  integrate  term-by-term  to  obtain  (12). 
It  is  also  possible  to  construct  a proof  of  (12)  along  these  lines. 

2 J . The  method  of  steepest  descents 

Ke  again  consider  the  integral 

(1)  /■(x)= g(t)  e dt. 
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in  which  we  now  assume  x to  be  a large  complex  variable,  g and  h to  be 
analytic  functions  of  the  complex  variable  t,  and  the  integral  to  be  taken 
along  some  path  in  the  complex  t plane.  This  integral  may  be  evaluated 
asymptotically  by  the  method  of  steepest  descents,  which  was  originated 
by  Riemann  and  developed  by  Debye.  Copson  (1946)  gives  a detailed 
description  of  this  method  together  with  references  and  several  examples. 

Those  points  of  the  t plane  at  which  h'{t)=  0 will  be  called  saddle 
points  or  cols.  The  surface  representing  |exp[*A(t)]|  as  a function  of 
Re  t and  Im  t will  be  called  the  relief  of  on  this  surface  cols  will 
be  “saddles”,  and  the  most  convenient  trail  from  one  “valley”  to  the 
other  will  lead  over  one  or  several  saddles.  Let  r be  a col:  if  h'ir)  = 
h"  (r)  = ••»=  k \r)  = G atad  k \r)  0,  we  call  r a coZ (or saddle  point) 

of  order  m.  In  the  t-plane,  curves  along  which  Re  xh  (t)  is  constant  are 
called  level  curves:  along  such  curves  has  a constant  modulus  (they 
are  contour  lines  of  the  relief),  and  the  phase  of  e'^  changes  as  rapidly 
as  possible.  Those  curves  along  which  Im  xk  {t)  is  constant  are  called 
steepest  paths : along  such  curves  e’^  has  a constant  phase,  and  the 
modulus  of  e’^  changes  as  rapidly  as  possible  (they  are  gradient  lines 
of  the  relief).  At  a col,  r,  of  order  m,  m + \ level  curves  intersect  at 
equal  angles,  and  their  angles  are  bisected  by  m + 1 steepest  paths: 
along  each  of  the  latter  curves  |c**^*^  has  a stationary  point  at  r. 

The  method  of  steepest  descents  consists  in  deforming  the  path  of 
integration  so  as  to  make  it  coincide  as  far  as  possible  with  arcs  of 
steepest  paths.  If  a and  /3  lie  on  steepest  arcs  through  cols,  for  instance 
if  a and  j8  are  singularities  of  h (t),  then  the  path  of  integration  n^ay  be 
deformed  so  as  to  consist  entirely  of  steepest  paths  through  cols;  other- 
wise two  steepest  arcs  may  occur  which  do  not  pass  through  cols.  This 
latter  case  may  be  described  by  reference  to  the  relief  by  saying  that 
we  first  descend  along  a gradient  line  to  a singularity  and  them,  climb  the 
saddle  along  another  gradient  line.  In  any  event.  Re  xh  {t)  is  monotonic 
along  any  steepest  path  (except  at  saddles),  and  Laplace’s  method  may 
be  used  to  evaluate  the  integral  asymptotically.  The  asymptotic  expan- 
sions of  g and  h needed  for  the  application  of  the  theoremi  in  the  pre- 
ceding section  are  the  Taylor  expansions  of  g and  h around  that  point  of 
the  steepest  path  at  which  Re  xh  (t)  is  a maximum  (this  is  often  the  col). 
The  inversion  of  the  series  2.4(10)  may  be  effected  by  Lagrange’s 
expansion  (see,  for  instance,  Copson  (1935),  p-  123”125)» 
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Meijer  (1933  Stb)  has  shown  that  nnmerica]  hounds  for  the  error  term 
may  be  obtained  by  using  Lagrange’s  expansion  with  a remainder,  and 
he  has  also  shown  that  in  some  cases  recurrence  relations  for  the  coef- 
ficients may  also  be  obtained. 

We  shall  consider  several  examples  of  the  application  of  the  method 
of  steepest  descents:  these  are  taken  from  Copson  (1946). 


2.6.  Airy’s  integral 

We  shall  investigate  the  asymptotic  behavior  of 


(1)  Ai(z) 


^ r ^ ^ 

■■ — / cos  I S + ZS  I 

^ Jo  \3  / 


ds 


for  large  positive  values  of  z.  With 
(2)  s = z t,  x~z^^ 


we  obtain 


(3) 


57/ 


and  the  method  of  steepest  descents  can  be  applied  to  the  integral  in 
(3).  In  (3),  t may  be  envisaged  as  a complex  variable  of  integration.  The 
path  of  integration  (for  x > 0)  is  the  real  t axis,  but  it  can  be  deformed 
into  any  curve  which  begins  at  infinity  in  the  sector  2n/3  < arg  t < n 
and  ends  at  infinity  in  the  sectw  0 < arg  t < Tt/Z.  Here 


h (t)  = i 


1), 


and  the  cols  are  the  zeros  of  h\t),  i.e.,  the  points  t=  ±i.  The  steepest 
paths  are  determined  bylm  A(t)=const.  We  set  t = ^ + ir)  and  obtain 

Im  hit)  •=  — i)  = 0, 

3 

so  that  the*^ equation  of  the  steepest  paths  is 

(4)  (ie  - 3t7*  + 3)  - 0. 
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This  equation  represents  a degenerate  cubic  consisting  of  the  imaginary 
axis  and  of  the  two  branches  of  a hyperbola.  In  the  figure,  arrows  indicate 
the  direction  in  which  Re  h {t)  decreases.  The  asymptotes  of  the  hyper- 
bola are  the  lines  ^ ± rj  \/3  = 0,  and  clearly,  the  path  of  integration  in 
(3)  can  be  deformed  into  the  upper  branch  of  the  hyperbola,  and  runs  from 
oo  • exp(5;7t'/6)  to  oo  • exp(t;r/6).  With  this  path,  the  integral  in  (3)  can  be 
seen  to  be  convergent  whenever  Re  x > 0. 

We  now  write 


(5) 


Ai{x^)=  f 


OO*  exp  (i‘Tr/6  ) 


^OO.  exp  (5  i7r/6) 


£Uid  evaluate  ^ by  Laplace’s  method.  In  both  integrals,  h{t)  — h(i)  is 
real  and  reaches  its  maximum  at  t » i;  also  A (t)  — h{i)  is  a decreasing 
function.  We  introduce  a new  variable  u by 


(6)  u=A(»)-AW 


2 

Z 


= u -o* 


0^ 
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From  (6) 
(7)  ±u^ 


where  is  the  positive  square  root,  [ is  that  value  which  reduces 

to  1 at  t = i,  and  the  upper  sign  in  (7)  holds  for  J,,  the  lower  sign  for 
It  follows  from  Lagrange’s  theorem  that  sufficiently  near  to  the  col, 
t — i possesses  an  expansion  of  the  form  t — i = 2 h (±  u^)”  where  nb  is 
the  coefficient  of  (t  — i)"~’  in  the  expansion  of  [1  — i{t  — i)/3]  in 
powers  of  t — i.  In  this  manner  the  expansions 


(8)  t-i=  y — - 

^ n ! 


~ r(3n/2  - 1) 


% \n 


n = I 


!r(n/2)  3 


(±a«) 


are  obtained  where  the  upper  or  lower  sign  holds  in  J,  and  respec- 
tively. Now 

J-~g-xu  ^ ^ 

1 . 2 Jo  I..  » 


and  accord ing  to  sec.  2.4  the  asymptotic  expansions  of  J,  ^ are  obtained 
by  substituting  (8)  in  dt/du  and  then  integrating  term-by-term.  Thus 


r(3n/2-l) 
2(n-  l)!r(a  /2)3""’ 


du 


y (±l)"t"~'  r(3«/2-  1) 
2(n- 1)!3""' 


Substituting  this  in  (5)  and  expressing  the  result  in  terms  of  z,  we  obtain 
after  some  simplification 


(9) 


r(3m  + »/2) 

(2m)! 


(-9  2"/^)"-, 


and  this  asymptotic  representation  holds,  uniformly  in  arg  z,  as  z -*  oo  in 
1 arg  2 1 < T7/3  — A,  A > 0. 
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2.7.  Further  exan vies 


We  shall  now  consider  two  examples  where  the  limits  of  integration 
are  not  singularities,  and  accordingly,  the  asymptotic  expansions  are  not 
obtained  by  expansions  around  the  col.  Also,  in  the  second  example,  the 
col  is  of  order  two. 

First  let  * > 0 and 


(1)  /"(*)=  Jg  exp 


The  function  h{t)  occurring  here  is  the  same  as  in  sec.  2.6,  the  steepest 
paths  are  those  shown  in  the  figure  in  sec.  2.6,  and  it  is  easily  seen 
that  the  appi'opriate  path  of  integration  consists  of  that  portion  of  the 
imaginary  axis  from  0 to  i £uid  then  one-half  of  the  upper  branch  of  the 
hyperbola.  Thus, 


(2) 


dt. 


The  asymptotic  expansion  of  the  second  integral  has  already  been 
obtained.  In  the  first  integral  h (t)  is  real  and  decreasing  as  t runs  from 
0 to  i,  and  we  may  again  use  Laplace’s  method.  Accordingly,  we  set 


u = h (0)  — h{t) 


and  infer  from  Lagrange’s  theorem  that  —it  = S 6^  c"  where  rib  ^ is  the 
coefficient  of  (— in  the  expansion  of  (1  + fV3)~"  in  powers  of  —it. 
Clearly  6^  = 0 if  n is  even  and 


if=  ^ 


(3m)!.u^ 


«=  0 


m 


\i-2rr.  + D!  3“ 


Substituting  this  in  the  first  integral  in  (2)  and  integrating  tenn-by-term, 

(3m)! 


* ^ zA  it) , 


= e-’  ^ 

0 -^0  .^o'«H2m)l3" 


f « 


du 


• I 


(3m)! 


-2. -I 


"*•  3" 
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It  is  seen  from  sec.  2.6  that  the  second  integral  in  (2)  is  exponentially 
small  in  comparison  with  the  first  one,  and  hence  the  result 


(3)  f ' 


exp 


f M 

txl  t + t 1 

L \3  yj 


dt  ~ 


y (3m)! 


-2.  -1 


as  X oo  in  S^,  A > 0. 

Our  last  example  is  the  integral 


(4)  /(x)=  /J  exp  dt 


where  we  take  x > 0.  Here  /i(i}  = and  { = 0 is  a col  of  order  two.  The 
steepest  paths  through  the  col  arc  the  lines  Im  (it^)  = 0,  that  is  the  lines 
arg  f = ± 7t/6,  ± n-/2,  ±5n^/6.  In  the  figure,  arrows  indicate  the  direction 
of  decreasing  jexpCixf^)]  for  x > 0.  None  of  these  steepest  paths  passes 
through  t = 1.  With  t = ^ + i t],  the  equation  of  the  steepest  path  through 
t = 1 is  Im(it^)  — 1,  or  — 3 = 1.  This  is  a cubic,  and  the  branch 

of  this  cubic  passing  through  t = 1 is  also  indicated  in  the  figtire.  In 
order  to  get  from  0 to  1 along  st»*epest  paths,  we  first  integrate  from  0 
to  oo  along  the  line  arg  t = 77/6,  and  then  from  00  to  1 along  the  upper  half 
of  the  branch  of  our  cubic.  Accordingly,  we  set 
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(5)  axp(, •«=)<;, 

-I,-h  . 

In  J,  we  set  u = or  t = where  u >0.  and  obtain 

(6)  J,  =—  X“e"“  rfn  = r(4/3)  e *""^®  x■'/^ 

3 


In  Jj,  we  set  u = — i(t^  — l)ort  = (l  + iuY^^  and  obtain 


1 oo 

I = _ I f e (1  + du  . 

z 3 0 

Expanding  (1  + in  the  binomial  expansion, 

^7)  Lx  s r(/n- 2/3)(txr"-’, 

^ r(-l/3)  r.=  0 

and  substituting  (6)  and  (7)  in  (5)  we  finally  have 
(8)  jJ  exp  (txt  df  ~ r (4/3)  e ® x' 

-TTTTTrr  2 r(n  + y3)(tx)-‘ 
r (—  1/3)  n = O 


as  X ->  oo  in  S^,  A > 0. 

The  last  equation  describes  the  asymptotic  behavior  of  f (x)  as  x -»  so 
in  the  right  half-plane.  If  x -►  ± loo,  the  integrand  in  (4)  is  real,  and 
Laplace’s  method  may  be  applied,  and  if  x is  in  the  left  half-plane,  we 
may  use  the  relation 

/(x)  = /vi)" 

which  follows  from  (4)  and  in  which  bars  denote  complex  conjugation. 

2.8.  Fourier  integrals 

Integrals  of  the  form 

(1) 

a 
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are  called  Fourier  integrals . T^e  shall  always  assume  that  (a,  /3)  is  a real 
interval,  and  mostly  (a,  /3)  will  be  a finite  real  interval;  and  will  be 
anintegrable  function  so  that  (l)exists  for  all  real  x.  We  shall  investigate 
the  asymptotic  behavior  of  (1)  as  x -*  +oo;  to  obtain  the  asymptotic  be- 
havior as  X -*  — oo;  replace  t by—t.  Unlike  in  the  case  of  Laplace  integrals 
(sec.  2.2)  it  seems  that  repeated  integrations  by  parts  is  the  only  effec- 
tive method  for  obtaining  asymptotic  expansions  of  (1),  except  in  the 
case  of  analytic  i^(t)  when  the  method  of  steepest  descents  may  be  used. 

First  we  shall  prove:  If  is  N times  continuously  differentiable 
for  a < t < /3  then. 


(2)  j/  e ^ 0(t)  dt^Bf,fx)-A„(x)^o  (x'^) 


where 


(3)  AAx)^  I,  r"-'  e 

^ fi=  0 


and  ^ ^ = d"<f>/dt The  result  remains  true  when  a=— «>(orj8  = cx))  pro- 
vided that  <f}  -*  0 as  t ■*  — «»  (or  t -»  oo)  for  each  n •=  0,  i,  , A'  — 1, 

and  provided  further  that  <p‘'^  is  integrable  over  (a,  /3).  To  prove  (2), 

we  apply  2.1  (5)  with  g = 0,  i A = e “*  h_^  = (tx)“"  e For  the 

remainder  we  obtain 


and  this  is  o{x~^)  since  the  integral  approaches  zero,  as  x -♦  oo,  by 
Riemann’s  lemma. 

We  note  that  A ^fx^  = 0 if  0 and  its  first  N — \ derivatives  vanish  at  a 
(for  instance,  if  </>  vanishes  identically  in  some  neighborhood  of  a),  and 
also  that  B ^{x)  = 0 if  and  its  first  IV  — 1 derivatives  vanish  at  /3  (for 
instance,  if  ^ vanishes  in  some  neighborhood  of  /3). 

We  now  turn  to  Fourier  integrals  whose  integrands  have  singularities 
of  a simple  type  at  one  end  point  of  the  interval. 
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,If  <f}U).is  N times  continuously  differentiable  for  a < t < ^;  = 0 

for  n = 0,  1,  ...  , - 1;  and  0 < A < 1;  then 

Q 

(4)  9!.(t)  d«  = -/i^(*)  + OU"^)  asx-*=c, 

where 

14 

(5)  Af^(.x)=  2 

n = 0 ^ * 

If  <^{t)is  N times  continuously  differentiable  for  a < t < ft;  <p^''\a.)-=0 
for  n = 0,  1,  ...  , /V  — 1;  and  0 < //  < 1;  then 

(6)  //  0(t)dt  = S^(*)+OU"^)  asx-oo, 

where 

(7)  fi„(x)=  f g7ri(„-M)/2 

^ a ! 

n = 0 

With  A = 1,  (5)  becomes  the  first  equation  (3),  and  with  (j.  = 1,  (7) 
becomes  the  second  equation  (3),  but  the  O terms  in  (4)  and  (6)  give  less 
information  than  the  o term  in  (2).  Instead  of  we  could  write 

o {x~^ in  (4),  and  o (x  ^ jj,  ajjj  these  latter  forms  remain 

valid,  and  pertinent,  when  A=  1 or  p = 1 respectively. 

We  shall  prove  l4):  the  proof  of  (6)  is  similar.  In  (4)  we  apply  2.1(5) 
withg(O=0(«),  g„(0=  <?i  <"’(«),  A(t)  = A„(f)=e“‘(t-a)^-\  and 

(8)  (0  = //”  (a  - tr  {u  - a)^-'  e *“  du 

n =0,  1,  ...  , A'  — 1. 

In  (8)  we  assume  t > a,  and  take  a path  of  integration  which  lies  entirely 
in  the  quadrant  0 < arg(u  — a)  < jt/2.  The  integral  converges  absolutely, 
and 


dt 


= A_  (0 

ft  1 n 


n •*  0,  1,  ...»  /V  - 1. 
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If  we  take  the  ray  u = t + * a,  a > 0,  as  the  path  of  integration,  we  have 
|n  — a|  > t — a,  |w  — a|^~*  < ~ a)^~'  for  0 < A < 1,  and  hence 


n ! 


/ 


t + ioo 
t 


u 


-tr 


le““dn| 


Substituting  u = t + i a we  have 

(9)  (01  <(«  - a)^"'  l>a,x>0. 

Also,  from  (8), 

n 1 n]  °- 

and,  with  u = a + r cr,  <7  > 0, 

(10) /,  ,(a)=  (-1)"^'  ^7r.<n4A)/r^-n-A^i,a  „ . Q,  1 

n! 


Ke  can  now  apply  2.1(5).  The  contributions  of  /3  to  the  vanish  since 
<jb^"^(/3)  = 0.  From  (10)  it  follows  that  2 s ^ ^ where  is  given  by 

(5).  Moreover,  from  2.1  (7), 

and  by  (9), 

l«^l  if  1«;6^’(01  (t-o)^*’  dt  = o(x“^). 

Ibis  proves  (4). 

Lastly  we  turn  to  Fourier  integrals  whose  integrands  have  singular- 
ities at  both  ends  of  the  interval. 

If  <^(t)  is  N times  continuously  differentiable  for  a < t < onrf 

0<A<1,0<I‘<1» 

(11)  e“*(t-a)^“’  (/3-0^”  0(t)dt  = 6^(x)-A^(x)  + 0(*-'^) 

a 


I 

1 


! 
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j 

I 


I 

I 


i 


I 


where 


r(n  + A) 

(12)  2 


g rri  (n +A“2  )/2  ^~n— A g 


1=  0 
N-y 


n ! 


(f" 


B 


r(n  + u)  _ . , 


)/ 2 e • 


n=  0 


n ! 


t//5" 


[(/S-a)^-'<5i.(.Q)], 


and  O {x  in  (11)  may  be  replaced  by  o{x  if  K = p = 1, 

This  theorem  contains  the  three  previous  results  as  special  cases.  To 
prove  (11),  we  shall  use  a device  which  is  frequently  employed  tor  such 
purposes,  and  is  called  a neutralizer  by  van  der  Corput.  Let  i/(t)  be  an 
infinitely  differentiable  function  for  a < t < fi,  v{a)  = 1,  ^^"ka)  = 0, 
n = 1,  2,  ...  ; = 0,  n = 0,  1,  2,  ...  . An  example  of  such  a function 

is 


With  such  a neutralizer  v{t)  we  write 

(13)  //  e “'(t  - a)^~'  {ft  - O'"*’  <?(0  dt 

= e (t  - a)^"’  [v(0  {ft  - O'"”  SS(0]  dt 
+ e “'(/3  - O'"”  1[1  - v(0]  {t  - ai^-^  0(t)i  dt. 

The  first  integral  on  the  right  heind  side  is  of  the  form  (4),  with  0(f) 
replaced  by  [ •••  1:  since  all  derivatives  of  this  function  vanish  at  t = ft, 
and  are  equal  to  the  corresponding  derivatives  of  {ft  — f)^~’  0(f)  at  t = a, 
we  obtain  the  expression  (12)  for  /4^(x).  Similarly,  the  second  integral 
on  the  right-hand  side  of  (13)  is  of  the  form  (6),  with  0(f),  replaced  by 
I I;  all  derivatives  of  1 •••!,(  vanish  at  f = oc,  and  are  equal  to  the 
corresponding  derivatives  of  (f  — a)^”’  0(t)  at  t = ft;  and  the  expression 
(12)  for  B jy{x)  follows  from  (7).  This  proves  (11).  If  A = p = 1,  then  (2) 
shows  that  0{x~^)  may  be  replaced  by  o (x”^). 
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All  our  results  remain  true  if  i is  replaced  by  —i  throughout  the  form- 
ulas, and  in  (8)  the  integration  from  t to  — t oo  is  taken  over  a path  which 
lies  in  the  quadrant  — jr/2  < arg(u  — a)  < 0.  Thus  we  are  able  to  describe 
the  behavior  of  our  integrals  as  x — oo,  and  also  the  asymptotic  behavior 
of  Fourier  integrals  with  trigonometric  kernels  cos  xt  and  sin  xt. 

2.9.  The' method  of  stationaiy  phase 

Ve  now  consider  the  integral 

(1)  f(x)= 


in  which  x is  a large  positive  variable  and  h (t)  is  a real  function  of  the 
real  variable  t.  According  to  Stokes  and  Kelvin,  the  major  contribution 
to  the  value  of  the  integral  arises  from  the  immediate  vicinity  of  the  eiid 
points  of  the  interval  and  from  the  vicinity  of  those  points  at  which  h (t) 
is  stationcay,  i.e.,  h'it)=  0;  and  in  the  first  approximation  the  contri- 
bution of  stationary  points,  if  there  are  any,  is  more  important  then  the 
contribution  of  the  end  points. 

Suppose  that  g is  continuous  and  h is  twice  continuously  differentiable, 
let  r be  the  only  stationary  point  of  h,  a < r < /S,  A '(r)  = Oand  h"(r)  > 0. 
In  the  assumption  that  the  neighborhood  of  r will  give  the  principal  con- 
tribution to  the  integral,  we  introduce  a new  variable  of  integration  u by 
the  substitution  h{t)  ~ h (r)  = u * and  obtain 


k '(0 


exp|rx[h(r)  + u*]}  du 


where  u,  = [h(r  - c)  - A,(r)]^,  = [A(r  + f)  - h(r)]*^.  Since  only  the 

nei^borhood  of  u = 0 matters,  we  may  replace  g (t)  by  g(r)  and  2u/h ' {t) 
by  [2/ h '' which  is  the  limit  of  2u/h ' (t)  as  t -*  r,  so  that 


exp  [ixu  ^ 


ixh  (r)]  du  • 


By  the  same  argument  we  may  extend  the  integration  from  — oo  to  oo  and 
finally  obtain 


(2) 


g (r)  exp[txh  (r)  + i n/ 4] 


as  X -*  BO, 
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which  is  virtually  Kelvin’s  result.  The  contribution  of  the  point  of  sta- 
tionary phase,  r,  to  the  integral  is  more  important  than  the  contributions 
of  the  end  points  because  the  latter  can  be  shown,  by  integration  by 
parts,  to  be  O (*”’)  if  A '(a)  0,  A ' (ft)  / 0. 

The  principle  of  stationary  phase  has  been  applied  to  numerous 
mathematical  and  physical  problems  but  it  appears  to  be  difficult  to 
formulate  it  in  a precise  manner.  Perhaps  the  best  available  theorem  is 
one  given  by  Watson  (1920).  Poincare'  discussed  the  principle  of  station- 
ary Dhase  applied  to  integrals  involving  analytic  functions,  and  the 
conn»-i  ion  of  his  work  with  the  method  of  steepest  descents  is  indicated 
in  Copson  (1946).  The  qiethod  of  stationary  phase  has  also  been  dis- 
cussed by  Bijl  (1937),  and  in  a much  more  general  setting  by  van  der 
Corput  (1934,  1936)* 

We  shall  use  the  discussion  of  Fourier  integrals  given  in  the  pre- 
ceding section  to  derive  a theorem  which  may  be  regarded  as  a precise 
version,  and  at  the  same  time  generalization,  of  (2).  A point  r at  which 
A '(r)  = k"(r  = ••».!=  A \t)  = 0 and  A ^ r)  0 will  be  called  a station- 
ary point  of  order  m,  m = 1,  2,  ...  . In  the  neighborhood  of  such  a point 
h’{t)  = it  - r)*  A,(t),  where  A,(r)  ^ 0.  The  notion  of  a stationary  point 
may  be  generalized  to  fractional  order.  A point  r will  be  called  a station- 
ary point  of  (.fractional)  order  p.  if  in  some  neighborhood  of  that  point 
A '(t)  is  either  of  the  form  |{  - r|^  A ,(t)  or  of  the  fonnsgn(f  - r)  |t-r|*^A  , (t), 
where  A j(r)  ^ 0.  Assuming  that  A (t)  has  at  most  a finite  number  of  sta- 
tionary points  (of  positive  order)  in  the  interval  under  consideration,  we 
may  break  up  the  integral  in  a finite  number  of  integrals  in  each  of  which 
A (t)  is  monotonic;  €uid  we  may  assume  A (f)  to  be  increasing.  Thus,  we 
shall  consider  integrals  of  the  form  (1)  in  which  A (t)  is  strictly  increas- 
ing when  a < t < /3,  and  a and  /3  are  either  ordinary  points  (i.e.,  station- 
ary points  of  order  zero),  or  stationary  points  (of  positive  order). 

//  0 < A,  p < 1;  g (t)  is  N times  continuously  differentiable  for  a<t<p-, 
A (t)  is  differentiable  and 

(3)  A ' (t)  = (t  - a)^‘*  (/3  - t)°^-'  A , (t), 

where  p,  o > 1,  and  A ^ (t)  is  positive  and  N times  continuously  differen- 
tiable for  a < t < /3;  then 

(4)  g(t)e^^*\t-a)^-'  (/3-  t)^-'  dt  = B(x)-A(.x) 

CL 
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where 

(5)  and  B(x)~fl^(ae)  to  terms  as  x -»  <», 

and  /4jy(x)  and  B are  given  by  (17)  and  (20)  below. 

In  the  proof  of  this  theorem  we  shall  use  the  abbreviation 

S 

(6)  g,(t)  = g(0(t-a)^-'  (/3-t)^-’. 

and  shall  employ  a neutralizer  v{t\  which  is  infinitely  differentiable  for 
° ^ ^ P such  that  for  some  77,  0 < 17  < (/S  — a)/2,  v(t)  = 1 when 
a < t ^ a + 77,  and  u{t)  = 0 when  /3  — 77<.t</3.  Ke  then  set 

(7)  - /I  (x)  = u(t)  g ,(t)  e ^ d( 

(8)  B(x)  = [1- v(t)]g,(t)e^<*’dt. 

To  obtain  an  asymptotic  expansion  of  A (x),  we  introduce  a new  vari- 
able of  integration,  u,  in  (7)  by 

(9)  u^^h(t)-h(a),  uf=AO-77)-A(a). 

From  (3)  we  have 

u^=  A(t)- A(a)  = fjh'(s)ds 

= (t  - a)^  t/3  - a — (t  - a)y]'^~’  A ,[a  + (t  - a)  y]  dy 

where  s = a + (t  — a)y.  The  last  integral  is  an  A'  + 1 times  continuously 
differentiable,  positive,  emd  increasing  function  of  t,  so  that  (9)  repre- 
sents an  yV  + 1 times  continuously  differentiable  mapping  of  the  interval 
^ ^ P ~ V onto  the  interval  0 < u < u , , and  the  inverse  mapping  is 

also  A’  + 1 times  continuously  differentiable. 

We  now  put  (u)  = vU)  and 

(10)  A (u)  = g (t)  u ’ — — , 

* rfti 
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where  ^,(f)  is  given  by  (6),  and  A (u)  is  N times  continuously  differen- 
tiable for  0 < u < u j . Then 

/I  (x)  = -c  ^ ' V-,  (u)  A (u)  exp  (ixuO  du 


can  be  integrated  by  parts  N times,  differentiating  and  integrating  the 
remaining  factor  of  the  integrand.  With 

(—  1 ^ ^ ^ 

(11)  0_„_,(“)=  (z  — u)"  exp(ix2^)d2 

the  result  of  the  integrations  by  parts  is 

A(x)  = A ^(x)  + R^(x), 

where 

(12)  ^ff(*)=^J^  (-l)"A<"’(0)^.^.,(0)e“'‘<‘"> 

and 

(13)  «^(x)  = (-l)^^'  ‘ <^.^(u)  ■■  du. 


In  (11),  the  path  of  integration  is  the  ray  argfx  — u)  = n/{2p)  in  the 
complex  plane.  Clearly, 


(-I)"-"'  A + A r^-ifrt+X)! 

(14)  . 


~ (n  +A.)/ p 

\-i  ^ .A-1 


To  estimate  (u)  for  u > 0,  we  note  that  jz]  ” < and  also  that 


+ X \z 


-u\P=ipx.J^  ^ 


\z  — u|  exp 


ni  Y*’ 

t;) 


d^. 


Since  the  real  part  of  the  last  expression  is  certainly  negative,  we  have 
I exp(ix2^|  < exp(— X \z  — u]^) 
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and  hence 


A-i 


(15)  (u)|  < ; — f \2  - s«|  " exp(-*  \z  - u\P)  d\z  - u\ 


n ! 


< _L  r 


)/p 


Alternatively,  the  method  of  steepest  descents  may  be  applied  to  (11)  to 
show  that 

(16)  <;6.„_,(a)  = u"^^0[(*a^)-»-’] 
for  large  *u'°. 

Substituting  (14)  and  (15)  in  (12)  and  (13)  we  obtain 
A‘"\0)  fn+x\  r ni(n  + A) 

V p y ’ L 2;rj 


n *=  0 


-(n*\)/Pp  ixhia) 
X ^ 


1 

.\x-N/p  r' u^-' 

d^iu^k) 

(N  - D!  \p 

) 1 

duf^ 

du. 


This  proves  that  A A ^ to  N terms  when  \ < 1.  When  A = 1 and  p = 1 
the  same  result  follows  from  sec.  2.8.  Let  A = 1 and  p > 1 and  choose  S 
so  that 

S 


(/V 


^<7)1. 


dHv^k) 

duf^ 


du  < 'A  { . 


Since  (16)  gives  0_jy(u)  = 0{x~^)  uniformly  in  u when  u > d,  we  have  for 
sufficiently  large  x 


d^(v,  k) 


du’ 


du  < Vt( 


so  that  = o(x”^^^)  also  in  this  case.  This  proves  the  result  for  A. 
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A similar  result  holds  for  B (x).  In  (8)  we  introduce  a new  variable  of 
integration  by 

(18)  t;‘^=A(/S)-A(0 
and  put 

(19)  /(t;)  = g.,(t)u'-'"^. 

dv 

where  g ^ (c)  is  given  by  (6).  In  the  repeated  integrals  of  v^~’  exp(—ixv‘^) 
we  integrate  along  the  ray  arg(z  — w)  = —Tri/(2a),  and  obtain  by  a process 
very  similar  to  that  used  in  the  case  of  A that  B B ^ to  N terms  where 


We  conclude  this  section  by  applying  the  general  result  to 
/■<*)*  exp  {ixt  ’)  dt. 

Here  l,p“3,  o-ljU**,  A (u)  = 1 and 

A^{x)^:  1/3  r(l/3)e^‘/*  **’'=* 

by  (17).  Also  V = 1-t^,  f = (l-  w)''^’, 
dt 


/(t,)=51  =_i 
dv  3 


and 


B 


V r(n 

2 777- 

n=  0 ' 


r(n  + ^3) 

V3) 


(ix) 


-n-\  g ix 


SO  that 

/■(x)~r(4/3)e^‘/‘x"’/^-  y Ei!LL?^(,x)-"-’  e“  as  x - c 

„i^o  r(-V3) 

which  is  the  expansion  obtained  in  sec.  2.7  by  the  method  of  steepest 
descents.  Note  that  x -»  oo  through  positive  values  with  our  present  n’.ethod 
while  X could  be  complex  in  sec.  2.7. 
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CHAPTER  m 

SIMGULARITIES  OF  DIFFERENTIAL  EQUATIONS 

In  this  chapter  we  give  a brief  introduction  to  the  asymptotic  theory  of 
ordinary  homogeneous  linear  differential  equations  of  the  second  order. 
Analogous  theories  exist  for  equations  of  arbitrary  (finite)  order  and  for 
systems  of  differential  equations  of  the  first  order.  For  these  more 
general  theories  see  Ince  (1927,  especially  p.  169 ff.,  p.  428  ff.,  p.  444  ff., 
p.  484ff.),Kamke  (1944,especially,  p.  17  ff.,  p.  60  ff.,  p.  100  ff.,  p.  133  ff.), 
Wasow  (1953),  the  references  given  in  these  works,  and  the  references 
given  at  the  end  of  the  present  chapter.  Asymptotic  expansions  occur 
also  in  connection  with  non-linear  differential  equations,  and  partial 
differential  equations. 

We  shall  investigate  the  asymptotic  behavior  of  solutions  of 
y"  + p(x)  y ' + q(x)  y * 0 

as  X -»  x^.  Here  * is  either  a real  variable  ranging  over  an  interval  (of 
which  is  usually  an  end-point),  or  else  a complex  variable  ranging 
over  a region  (of  which  x^  is  often  a boundary  point).  Without  loss  of 
generality,  we  take  x^  = oo  throughout  this  chapter. 

The  reader  is  expected  to  know  the  basic  existence  theorems  regard- 
ing the  above  differential  equation  both  for  real  and  complex  variables, 
and  he  is  also  expected  to  be  familiar  with  the  fundamental  properties  of 
the  solutions. 

3.1.  Classification  of  singularities 

In  the  present  section  we  shall  discuss  the  differential  equation 
(1)  y " + p (x)  y ' + q(*)  y = 0 

when  X is  a complex  variable  ranging  over  an  Euinular  region,  R,  given  by 
r < |x|  < oo,  and  p (x)  and  q (x)  are  single-valued  analytic  functions  in  R 
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(which  may  or  may  not  have  singularities  at  <»).  We  shall  briefly  review 
the  well-known  classification  (see,  for  instance,  Poole,  1936,  Chapter 
IV)  of  isolated  singularities  of  (1),  the  singularity  in  question  being  at  eo. 

If  y I (pc)  and  y ^ (pc)  are  two  linearly  independent  solutions  of  (1)  form- 
ing a fundamental  system,  and  if  we  continue  these  functions  analytically 
along  some  curve  in  R which  begins  and  ends  at  x and  encircles  cm  in  the 
positive  sense,  we  obtain  two  new  functions  which  may  be  denoted  by 
y j(xe~^^'),  / = 1,  2.  7’hese  need  not  be  identical  with  the  y .(pc),  but  at 
any  rate,  they  will  be  solutions  of  (1),  so  that  relations  of  the  form 

y,(*e"*^*)=  a„  y,(*)  + a,^  y ^(x) 

(2) 

y^(xe~^'”  *)=a2,  y,M+  a^y^(pc) 
will  hold.  In  (2), 


12 


a 


22 


is  a constant  non-singular  matrix. 

If  instead  of  y,,  y^  we  take  another  fundamental  system,  we  obtain  a 
matrix  B,  and  an  easy  computation  shows  that  B = MAM"',  where  ^ is  a 
non-singular  constant  matrix.  Thus,  all  matrices  obtained  in  this  manner 
have  the  same  affine  invariants , in  particular,  the  same  latent  roots,  and 
the  same  canonical  forms.  These,  then,  are  independent  of  the  fundamen- 
tal system  chosen,  and  are  characteristic  of  the  singularity  at  oo  (if  there 
is  a singuleu'ity  there). 

Let  us  assume  that  the  latent  roots  of  A are  distinct  so  thatwe  have  a 
diagonal  canonical  form 


If  y,,  yj  is  the  fundamental  system  corresponding  to  the  canonical  form 
of  A (canonical  fundamental  system),  then  (2)  assumes  the  form 

(3)  y^.(xe"*^*)  = A^.y^.(*) 


7=1,?. 
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We  now  set  = exp(2n’t  p^.),  and  call  p,,  p^  the  exponents  belonging  to 
oo:  these  are  determined  up  to  an  integer  in  each.  From  (3)  we  see  that 
the  canonical  fiuidamental  system  is  of  the  form 

(4)  y/*)=x"^>0.(x)  / = 1,2 

where  tjj^  are  single-valued  analytic  functions  of  x in  /?  with,  possibly, 
singularities  -at  cx>. 

If  the  latent  roots  are  equal,  = A = exp  (2  rri  p),  then  the  canon- 

ical fundamental  system  can  be  shown  to  be  of  the  form 

y , (x)  = x~^  0,  (x) 

(5) 

yj(x)=  cy,  (x)  logx  + x'^i/f^fx), 

where  \jj ^ and  xfj ^ have  the  same  properties  as  in  (4),  and  c is  a constant: 
c = 0 or  c ^ 0 according  as  the  canonical  form  of  A in  this  case  is  dia- 
gonal or  not. 

X = oo  is  called  an  ordinary  point  of  (1)  if  all  solutions  are  regular  at 

oo,  i.e.,  can  be  represented  by  convergent  power  series  in  x”';  x = <»  is 

called  a regular  singularity  of  (1)  if  it  is  not  an  ordinary  point  and  if 
0,  and  0 2 liave  at  most  poles  at  oo  so  that  by  a suitable  choice  of  p, 

and  P2,  0,  and  02  can  be  made  regular  at  00;  and  x = 00  is  called  an 

irregular  sirigularity  of  (1)  if  at  least  one  of  the  two  functions  0|,  02  has 
an  essentied  singularity  at  00. 

It  can  be  shown  (see,  for  instance,  Poole,  1936,  section  20)  that  a 
sufficient  condition  for  x = 00  to  be  br  ordinary  point  is 

(6)  p (x)  = 2x“’  + (?(x“*),  q{x)=0{x~*)  as  x -►  00, 

and  that  a sufficient  condition  for  x = 00  to  be  a regular  s ingularity  is 

(7)  p(x)=0(x”'),  q(x)  = 0(x”*)  as  X -►  00. 

In  the  case  of  an  irregular  singularity , p and  q may  have  essential  singu- 
Icirities  at  00:  if  p and  q have  at  most  poles  at  00,  we  speak  of  an  irregular 
singularity  of  finite  rank,  and  the  least  integer  k for  which 

(8)  p(x)  = 0(x*”’),  q(x)  = 0(x“~*)  asx-»eo 

is  called  the  rank  of  the  irregular  singularity.  Sometimes,  a regular  singu- 
larity is  regarded  as  a singularity  of  rank  zero. 
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3.2.  Nonnal  BolotioDS 

If  X = oo  is  an  ordinary  point  of  3.1  (l)f  we  may  expand  y in  a series 
of  powers  of  x ' . The  coefficients  of  that  series  may  be  determined  from 
recurrence  relations,  and  the  series  converges  in  f?.  If  x = oo  is  a regular 
singularity  of'3.1(l),  we  may  set 

y=  S 

n=  0 

We  then  obtain  a quadratic  equation  for  p,  recurrence  relations  for  the  c 
(similar  to  (7)  below),  and  a series  for  y which  converges  in  f?.  In  either 
of  these  two  cases  the  coefficients  can  be  computed  easily,  and  the  con- 
vergent series  can  be  used  with  ctdvantage  to  compute  the  solutions  for 
large  x. 

The  situation  is  entirely  different  if  x = eo  is  an  irregular  singular 
point.  Since  the  0^.  have  essential  singularities  at  oo,  we  must  set  in  this 
case 


y = 


C X 
n 


-p-n 


For  the  we  obtain  an  infinite  system  of  linear  equations  which  cannot 
be  solved  recurrently,  and  for  p,  a transcendental  equation  which  involves 
an  infinite  determinant  (the  determinant  of  the  system).  In  this  case  the 
coefficients  cannot  be  computed  easily,  nor  is  the  series  rapidly  con- 
vergent for  large  x. 

It  was  discovered  by  Thomd  that  in  the  case  of  an  irregular  singularity, 
of  finite  rank  certain  formal  solutions  exist  which  do  not  suffer  from  the 
disadvantages  mentioned  above;  the  coefficients  occurring  in  these 
solutions  can  be  computed  recurrently,  and  the  series  appear  suitable  for 
numerical  computations  for  large  x.  Thome’s  solutions  are  of  the  form 


y = exp  [P  (x)]  2 x“^",  ^ 0, 

ns=  0 

where  P (x)  is  a polynomial:  they  are  known  as  normal  solutions. 

We  shall  explain  the  construction  of  normal  solutions  in  the  case  of  an 
irregular  singularity  of  rank  one.  We  first  note  that  setting 


y = z exp  (-  Vi  f p dx) 
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in  3.1  (l)f  we  obtain  for  z a differential  equation  of  the  form  (1)  in  which 
z 'does  not  occur.  Hence  it  will  be  sufficient  to  discuss  the  differential 
equation 

(1)  y"  + q{.x)  y = 0 
in  which 

(2)  q(*)=  2 

0 

the  series  being  convergent  in  /?.  Wfe  shall  attempt  to  find  formal  solutions 
of  the  form 


(3)  y = e‘ 


n*  0 


<^0  ^ 


where  ej  and  p are  constants,  ft  is  to  be  noted  that  p in  (3)  is  not  nec- 
essarily one  of  the  exponents  belonging  (in  the  sense  of  the  analytical 
theory  of  sec.  3.1)  to  the  irregular  singularity  x = oo. 

In  manipulating  the  formal  series,  we  shall  adopt  the  convention 
q_^  = 0»  = 0,  m = 1,  2,  ...  , so  that  all  summations  may  be  extended 

from  — oo  to  +<M  except  those  explicitly  stated  otherwise. 

Substituting  (2)  and  (3)  in  (1),  we  have 

2 c — 2<u  2 (p  + n)  c 

+ 2(p  + a)(p  + n+  l)c  + .V  q x-"2  c x“^~"  = 0. 

Comparing  coefficients  here,  we  obtain 

(4)  c — 2tu(p  + a - 1)  c + (p  + n - 2)  (p  + R — 1)  c , 

A ' fl  I * * A Z 


+ 


c 


n—v 


= 0 


for  all  integers  n.  The  first  non-vacuous  condition  arises  when  .1  = 0 in 

(4) .  Since  0 we  have 

(5)  cj*  + = 0. 
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If  n «:  1 in  (4),  and  a satisfies  (5),  we  obtain 

(6)  — 2(iip  + 9 , »=  0. 

These  two  equations  determine  a>  and  p.  The  recuirence  relation  for  the 
coefficients  may  also  be  obtained  from  (4)>  Ve  replace  n by  n + 1,  and 
use  (5)  and  (6)  to  obtain 

n + 1 

(7)  2cc)  nc^  = (p  + n)(p  + 71  - 1)  c^_,  + t*  = 1,  2, 

l\e  now  see  that  normal  solutions  exist  if  either  or  = 0. 

In  the  former  case,  (5)  determines  oi,  (6)  determines  p,  and  with  = 1, 

(7)  determines  the  coefficients.  Moreover,  m,  p,  c ,,  ...  , are  com- 
pletely determined  by  9,,  ...  , 9,+,,  and  vice  versa.  There  are  in 
this  case  always  two  normal  solutions  corresponding  to  the  two  possible 
values  of  oi.  In  the  latter  case, which  is  the  case  of  a regular  singularity 
at  oo,  we  have  co  = 0 from  (5),  equation  (6)  is  vacuous,  (7)  with  n = 1 
determines  p as  one  of  the  roots  of  the  quadratic  equation  p(pt- 1)  + 

and  (7)  with  n = 2,  3,  ...  determines  the  coefficients. 

If  9o  = ^ and  q,  7^  0,  then  (5)  and  (6)  cannot  be  satisfied,  and  there 
exists  no  normal  solution.  However,  subnormal  solutions  may  be  obtained 
by  tr2uisforming  (1)  into 

(8)  q"  7=0 

by  the  change  of  variables 

q(^)=  y(*X 

If  9o  ~ 0 and  q ^ ^ 0,  then  (8)  has  an  irregular  singularity  of  rank  one  at 
eo,  and  possesses  normal  solutions.  These  give  rise  to  subnormal  solu- 
tions of  (1)  having  the  form 

(9)  y = exp(<u*^)  2) 


in  this  case. 
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For  the  construction  of  normal  and  subnormal  solutions  for  singularities 
ofhigher  ranks,  for  differential  equations  ofhigher  orders,  and  for  systems 
of  differential  equations  see  Ince  (1927,  p.  423  If.,  p.  427flF«  p.  46911., 
p.  478  ff.). 

Normal  and  subnormal  solutions  are  formal  solutions;  i.e.,  if  they  are 
substituted  in  the  differential  equation  as  if  the  infinite  series  were 
convergent,  the  differential  equation  is  satisfied.  However,  the  infinite 
series  involved  in  the  formal  solutions  are  in  general  divergent  .Never- 
theless, they  £ire  far  from  being  useless,  for  they  represent  asymptotic 
expansions  of  solutions  of  (1).  This  situation  has  been  investigated  by 
many  authors,  beginning  with  Poiacar6:  some  of  the  papers  are  listed  at 
the  end  of  this  chapter,  and  it  may  be  noted  that  the  must  general  results 
were  obtained  by  Sternberg  (1920)  for  the  differential  equation  of  order  n 
and  arbitrary  (finite)  rank,  and  by  Trjitzinsky  (1933)  for  a system  of  first 
order  differential  equations. 

By  and  large  there  seem  to  be  two  methods  for  proving  that  the  formal 
solutions  are  eisymptotic  expansions  of  solutions  of  the  differential 
equation.  One  of  these  was  originated  by  Poincar6  and  was  developed  by 
Horn  in  a large  number  of  papers  of  which  we  refer  but  to  a few  at  the  end 
of  this  chapter.  This  method  consists  in  finding  lateral  representations 
of  Laplace’s  type  for  the  solutions,  and  then  basing  the  asymptotic 
expansions  on  the  work  of  sec.  2.2.  The  other  method,which  was  developed 
by  C.D.  Birkhoff  and  his  pupils,  uses  the  leading  terms,  or  partial  sums, 
of  the  formal  solutions  to  construct  a differential  equation  which  in  a 
certain  sense  is  close  to  the  given  equation  when  x is  large,  and  then 
compares  the  two  equations.  As  it  happens,  singular  Volterra  integral 
equations  or  integro-differential  equations  play  an  important  part  in  both 
methods. 

We  shall  use  a variant  of  the  second  method  to  discuss  the  differential 
equation  (1)  with  The  proof  to  be  given  below  is  based  principally 

on  the  work  of  Hoheisel  (1924)  and  Tricomi  (1953,  sections  47  to  50). 
Since  the  analytic  character  of  q and  y does  not  enter  in  the  formal 
solution,  the  investigation  may  be  carried  out  either  for  real  or  for  com.- 
plex  independent  variables. 

3.3.  The  integral  eqiutioa  and  its  solution 

We  first  consider  the  case  of  a real  variable,  and  defer  a brief  dis- 
cussion of  the  case  of  a complex  variable  to  sec.  3.5.  Let  us  then  con- 
sider 


(1)  y"  + ?(*)  y = 0 


X >'a  > 0, 
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assuming  that  q (x)  is  continuous  forx  > a and 

(2)  9(*)~  X asx-»oo,  7 ^0. 

n=  0 

\^e  then  obtain  two  forma]  solutions 

(3)  X c^x-P-\  c.^0 

n=  0 

where  <u,  p,  c,,  c^,  ...  satisfy  3.2(5),  (6),  (7).  We  shall  show  that  these 
formal  solutions  are  asymptotic  expansions  of  certain  solutions  of  (1). 

Let  (o  = tj,  + (Wj,  p = p^  + ip^,  and  determine  a = (— so  as  to 
make  x ^ bounded  as  x -»  oo.  If  is  not  positive  real,  we  take  that 
value  of  the  square  root  which  makes  co  ^ < 0;  if  > 0 and  Im  q , 0, 

we  take  that  value  of  the  square  root  which  makes  p^  = Re[q,/(2£L))]  > 0; 
and  if^Q  > 0 and  q ^ is  real,  we  take  either  value  of  the  square  root.  Thus, 
we  always  have  either  < 0 or  tu,  = 0 and  p,  > 0,  These  conventions 
will  be  retained  throughout  the  discussion. 

It  will  be  convenient  to  transform  (1)  by  setting 

(4)  y (x)  = e"*  x"^z  (x) 

so  that  z satisfies  the  differential  equation 

,5)  j ..0. 

Here  <j  and  p satisfy  3.2  (5),  (6).  We  put 

(6)  x^[o(x)  -q^-q  y x“’]  + p(p  + 1)=  F{.x\ 
and  see  from  (2)  that  F (x)  is  bounded,  say 

(7)  |F(x)1<A, 

We  now  rewrite  (5)  as 

+ x-^P-^  F(x)  2(x)  = 0, 


d f , on  dz 

( e*"*x  

dx  \ dx  . 


X > a. 
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integrate  to  obtain 

dx  ° 

where  c ^ and  b > a are  arbitrary  constants,  and  integrate  once  more  to 
obtain 

(8)  zCx)  + K(x,  t)  F(t)  z(t)  dt  = 

where 

/ 5 

(9)  K(Xft)=—f  exp[26»(t  — s)]  f — J ds. 

Equation  (8)  is  an  integral  equation  of  Volterra’s  type.  Any  solution 
of  (5)  satisfies  (8)  for  some  b,  c,,  Cjj,  and  conversely,  any  twice  con- 
tinuously differentiable  solution  of  (8),  for  any  b,  c,,  c satisfies  (5). 
The  existence  of  such  a solution  follows  from  the  general  theory  of  integral 
equations  when  b < oo.  When  6 = oo,  the  integral  equation  (8)  in  a singular 
integral  equation,  and  the  existence  and  differentiability  of  the  solution 
must  be  demonstrated. 

In  order  to  prove  that  (5)  possesses  a solution  which  can  be  repre- 
sented asjmiptotically  by  2 we  take  6 = oo,  c , = 1,  c ^ = 0 in  (8) 

so  that  the  integral  equation  (8)  becomes 

(10)  z(*)=l  + Kix,  t)  Fit)  zit)  dt. 

This  integral  equation  will  be  solved  by  the  method  of  successive  approxi- 
mations. 

For  any  function,  ^ix),  we  set 

(11)  TCix)  = Kix,  t)  Fit)  Cit)  dt, 
and  then  define 

(12)  Zj,(x)=l,  z^,  (x)  = Tz^fx)  n = 0,  1,  2,  ... 

(13)  z(x)-  2 

ns  0 


z^(x). 
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It  will  now  be  proved  that  z (x)  exists,  satisfies  (10)«  is  differentiable, 
and  satisfies  (5).  Tbe  proof  will  be  conducted  in  several  steps. 

The  kernel,  K (x,  t),  is  bounded  for  t > x > where  x^>  a and  x^  is 
sufficiently  large. 

Proof:  Since  either  o),  <0  or  <Uj  = 0 and  p,  > 0 we  have 


d 

ds 


— 2«i),s  2p,. 

log(c  s ’)  = — 2<i)|  + 


>0 


for  sufficiently  large  s;  and  hence 
-2«o,s  2p, 

€ S 


is  an  increasing  function  of  s . We  now  write 


2oit  — s)  / s\lp  2«i>,(t-s)  / s 

; \T/  vT"/ 


where 


2tp, 


, . / X 2»  a>,(t  s)  / S \ 

^,(s,  t)+ t)-=  e vT  y ’ 

and  apply  the  second  mean  value  theorem  to  (9),  obtaining 

. 2p, 

r t 2&>, (t  ~ s)  / s \ ^ 

-K{x,t)=J^e  (s.  t)  + 

\ f / 


= e 


2p, 

2a)  {t  - x)  / X \ - ^ T) 

f— J [/^  <f^(s,t)ds+if^  <^^(s,t)ds] 


+ 0 ds  + i J As,  t)ds, 

€ V 

where  x < rf  < t.  The  integrals  on  the  right-hand  side  are  bounded 
functions  of  x and  t,  and 
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when  *>*>*(,  and  is  sufficiently  large,  so  that 

(14)  \K{.x,  t)]  < B,  t>x>x^ 
for  some  and  fi. 

If  K(t)|  < Ct  for  t > *g,  where  A > - 1,  then 
ARC 

(15)  ircU)!^- — *>x„. 

A + 1 

Ue  have  by  (11),  (7),  and- (14) 

\TW\  = \r  K{x,t)F{t)at)t-^dt\  < ABC  dt. 


and  this  proves  (15). 

For  the  functions  defined  by  (12)  we  have 


(16) 


k „(*)!< 


X > Xq- 


Proof  by  induction.  (16)  is  true  for  n = 0,  and  if  it  is  true  for  any  n,  the 
definition  of  + , combined  with  (15)  shows  that  it  is  also  true  for  n + 1, 
The  series  (13)  converges  uniformly  for  x > x^,  and  the  function  z{x) 
satisfies  (10).  Moreover,  z(x)  is  twice  continuously  differentiable  and 
satisfies  (5).  The  uniform  convergence  follows  from  (16).  If  we  substitute 
2 = S 2^  in  the  integral  in  (10),  term-by-term  integration  is  justified  by 
uniform  convergence,  and  shows  that  (10)  is  satisfied.  Furthermore,  the 
integral  in  (10)  is  a differentiable  function  of  x,  and  so  is  z (x).  Since 


K(x,  x)  = 0. 


dK{x,  t) 

dx 


we  obtain  from  (10), 

(17)  2' (x)=.j“e^"‘*-*>^-^)  ^ Fit)  z it)  r^dt. 


The  last  integral  is  again  a differentiable  fiuiction  of  x,  and  substitution 
sliows  that  2 (x)  satisfies  (5).  With  z given  by  (13), 
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satisfies  (1).  If  is  positive  and  9,  is  real,  we  may  take  either  of  the 
two  values  of  (— for  and  we  thus  obtain  two  linearly  independent 
solutions  of  the  form  (18).  In  every  other  case,  there  is  only  one  solution 
of  this  form,  and  a second  solution  may  be  written  down  in  the  form 

(19)  y2(*)  = y|(*)  /b*  [y I 

where  b is  any  number  large  enough  to  ensure  that  y , (*)  ^ 0 for  x > b. 
Since  2 (x)  = 1 + 0(x  ')  as  x-»  00,  such  a b certainly  exists. 

Thus  in  every  case  we  have  two  linearly  independent  solutions  of 
(1)  in  the  interval  x > and,  if  a < x^,  these  solutions  can  be  extended 
to  the  interval  x > a.  Itremainsto  show  that  the  formal  solutions  obtained 
in  sec.  3.2  are  asymptotic  expansions  of  the  solutions  obtained  in  this 
section. 

3.4.  Asymptotic  expansions  of  the  solutions 

fte  first  remark  that 


(1)  2 (-l)-(v)  X- 


as  X 00 


(2)  f 


as  X -*  00, 


where 


(3)  (»^)o=  1,  (1^)^=  1^(1/ + 1)  (v  + r - 1), 


r — 1,  2,  ...  . 


Both  results  can  be  proved  by  successive  integrations  by  parts,  with 
g = = (+  !)■  exp(+  t)  in  2,1(5).  In  particular, 

(4)  J dt  = 0{e~*  x~^  as  X 00 


(5)  e * dt  ==  Oie*  X as  x 

6 
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Next  we  prove  by  induction  that  the  functions  defined  by  3.3(12) 
possess  asymptotic  power  series  expansions  of  the  form 

OO 

(6)  z (x)  c.x*  asx-»oo. 

a 

n 

This  is  certainly  true  for  n = 0.  If  it  is  true  for  any  n,  then 
F {t)  z (t) 'S.  a . t~  as  t -►  OO. 

n , j| 

k — n 

Also 

If  <ii  and  p are  real,  the  last  theorem  in  sec.  1.4  justifies  substitution 
of  the  asymptotic  expansion  of  F (t)  z ^(t)  in  the  integral,  so  that 

(7)  z^^,(x)~^S  ^ t~'‘~^  dt  asx-*oo. 

If  <t)  or  p is  complex,  (7)  can  be  proved  by  taking  the  asymptotic  expansion 
of  F (t)  2„(<)  to  a finite  number  of  terms,  with  a remainder,  substituting, 
and  estimating  the  remainder  in  (7)  by  (4).  In  any  event,  each  integral  in 

(7)  possesses  an  asymptotic  power  series  expansion  which  can  be  obtained 
from  (1),  and  starts  with  x~  . By  the  third  theorem  in  sec.  1.4,  it  is 
permissible  to  substitute  this  expansion  in  (7):  p(n)  = n,  the  uniformity  of 
the  asymptotic  expansion  is  trivial,  and  the  series  14(5)  tenninate  so 
that  the  question  of  convergence  does  not  arise.  Thus  we  obtain  by 
rearrangement 

(8)  <.,(*)-! 

n 

and  by  integration  of  (8)  we  obtain  the  asymptotic  expansion  of 

By  the  second  theorem  of  sec.  1.4,  it  is  permissible  to  substitute  (6) 
in  33(13):  since  z^  = 0(*~"),  the  question  of  uniformity  ’is  trivial,  and 
the  series  1.4(2)  terminate  so  that  the  question  of  their  convergence  does 
not  arise.  Thus  we  see  that  z (x)  possesses  an  asymptotic  power  series 
expansion  of  the  form 

(9)  z (x)  ~ S 


as  X -»  OO, 
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It  remains  to  prove  that  the  coefficients  c^  occurring  here  satisfy  3.2(7). 
It  follows  from  3 .3  (17)  and  the  corresponding  relation  for  z " , that  z '■  and 
z " also  possess  asymptotic  power  series  expansions.  By  a result  in 
sec.  1.6  it  follows  that  (9)  may  be  differentiated  twice.  The  resulting 
asymptotic  series  must  satisfy  3.3(5)  formally,  and  this  leads  to  3.2(7) 
for  the  coefficients.  Also,  c^  = 1. 

We  have  thus  proved  that  3.3(18)  is  represented  asymptotically  by 
one  of  the  formal  solutions,  and  we  conclude  by  showing  that  3.3(19)  is 
represented  by  ^e  other.  To  do  this  we  put 

(10)  y2(*)=c~"*x^Z2(x), 
and  have  from  3.3(18)  and  (19) 

(11)  z^(x)  = z(x)J‘  ^ [z(0r"dt. 


By  the  choice  of  6,  z (t)  is  bounded  away  from  zero,  and  we  have  also 
seen  that  z (t)  possesses  an  asymptotic  expansion,  (9)  with  w 1,  as 
t -*  oe.  Then 


U(t)]' 


a t-"+0(r^)  fort  >6 


for  some  a , and 

. "i'a  r 

z (x)  n=  6 " * 


The  integral  under  the  summation  sign  can  be  expanded  asymptotically  by 
(2),  and  the  last  integral  is  0 (x~^)  by  (5).  Since  /V  is  an  arbitrary  positive 
integer,  z ^{^)/z  (x)  possesses  an  asymptotic  power  series.  Bence  z^Oic) 
possesses  an  asymptotic  power  series  expansion 


(12)  2.(*)~  S c^x"” 

n=  0 


as  X eo  . 
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By>  a similar  consideration  as  in  the  case  of  z ,,  it  can  be  proved  that 
the  coefficients  occurring  here  satisfy  a recurrence  relation  which  differs 
from  3.2(7)  only  in  that  oj,  p are  replaced  by  —co,  —p.  Thus,  3.3(19)  is 
represented  asymptotically  by  one  of  the  formal  solutions. 

If  &>  and  p are  both  imaginary,  the  two  fundamental  solutions,  both  of 
the  form  3.3(18),  are  defined  uniquely  up  to  a constant  factor:  both  are 
bounded  and  neither  approaches  zero  as  a:  ->  oo.  In  all  other  cases,  one  of 
the  fundamental  solutions,  3.3(18),  approaches  zero  as  x ->  oo,  and  is 
defined  uniquely  up  to  a constant  factor:  the  other,  3.3(19),  is  unbounded 
as  X -►  oo,  and  is  not  unique  (since  it  depends  on  b).  In  fact, 

X]  y , W + Xz  ^2^*^  ^2  5^2  X -»  oo 

3dS.  Complex  variable.  Stokes’ phenomenon 

The  results  of  the  preceding  sections  may  be  extended  to  the  case  of 
a complex  variable  x ranging  over  a sectorial  region  5, 

(1)  |x|  > a,  a < arg  x < /3. 

It  will  be  assumed  that  q (x)  is  analytic  in  5,  and  that  3.3(2)  holds, 
uniformly  in  arg  x,  as  x « in  S.  We  then  have  two  formal  solutions 
3.3(3),  where  <u  satisfies  3.2(5). 

The  line  Re[x(-q^)*^]  = 0 is  called  the  critical  line  or  the  Stokes  line. 
If  X -*  oa  along  one  of  the  rays  of  the  critical  line,  the  exponential  factors 
in  both  formal  solutions  remain  bounded,  and  bounded  away  from  zero. 
If  x-*oo  along  any  otherray,  the  leadingterm  in  one  of  the  formal  solutions 
increases  exponentially. 

First  we  assume  that  the  critical  line  does  not  intersect  5.  Cleemly,  in 
this  case  /3  — a < tt,  and  we  may  take  a>  as  that  solution  of  3 .2  (.5)  for  which 
Re  cox  < 0 for  all  x in  S.  If  x varies  along  any  ray  arg  x = const,  in  S, 
the  results  of  sections  3.3  and  3.4  held,  and  these  results  can  be  extended 
to  the  sector  S as  follows.  In  the  integral  equation  3.3(10)  we  always 
integrate  along  a ray,  so  that  arg  x = arg  t.  The  boundedness  of  the 
kernel  then  follows  for  each  x in  S,  and  uniformly  in  x when  a < arg  x < /3. 
The  integral  equation  can  be  solved  as  before,  each  z can  be  shown 
to  be  analytic  in  5,  and  z (x)  is  also  analytic  in  5,  since  it  is  the  uniform 
limit  of  analytic  fimctions.  In  3.3(8),  b is  chosen  so  that  y , (x)  0 when 

l*i  > ^ < /3.  The  result  is  the  existence  of  two  solutions,  y, 

and  y^,  in  S which  are  represented  asymptotically  by  multioles  of  the 
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formal  solutions  3.3  (3).  The  asymptotic  expansions  hold  uniformly  in 
arg  X as  X -*  oo,  a < arg  x < /9.  Any  solution  of  the  differential  equation 
is  a linear  combination  of  y ^ and  y^;  and  its  asymptotic  expansion  follows 
from  the  asymptotic  expansions  of  y,  and  y^. 

Next  we  assume  that  the  critical  line  intersects  5 and  decomposes  it 
into  a finite  number  of  sectors,  A = 1,  , K,  and  certain  rays  of  the 

line  itself.  In  each  of  the  sectors  we  have  a value  of  such  that 
Re  X < 0 for  all  x in  and  in  each  Sj  we  have  a fundamental  system 
y I*  ’ y 2k  asymptotically  represented  by  the  formal  solutions.  For 

a fuller  discussion  of  these  fundamental  systems  the  reader  is  referred 
to  lloheisel  (1924).  It  turns  out  that  the  fundamental  system  belonging  to 
a ray  of  the  critical  line  may  be  taken  also  as  a fundamental  system  for 
the  two  sectors  separated  by  that  ray.  Each  of  the  two  solutions  is 
dominant  (exponentially  increasing)  in  one  of  the  two  sectors,  and  re- 
cessive (exponentially  decreasing)  in  the  other. 

Let  us  consider  a solution,  y (x),  of  the  differential  equation  in  S.  In 
each  of  the  sectors  y is  a linear  combination  of  the  two  fundamental 
solutions  for  that  sector;  in  each  of  the  sectors  y will  b e represented 
asymptotically  by  a linear  combination  of  the  two  formal  solutions;  but 
the  coefficients  may  vary  from  sector  to  sector.  This  circumstance  was 
discovered  by  Stokes,  and  it  is  called  Stokes’  phenomenon.  The  sectors 
Sj  are  sometimes  called  Stokes  sectors,  and  the  critical  rays,  Stokes 
rays. 

For  the  determination  of  the  coefficients  involved  in  the  expression 
of  y (x)  as  a linear  combination  of  the  formal  solutions  see  'lurrittin 
(1950). 

3.6.  Bessel  functions  of  order  zero 

I^e  shall  illustrate  the  results  of  the  last  few  sections  by  a brief 
discussion  of  the  differential  equation 

(1)  2 " + x“‘  Z ' + 2 = 0 

satisfied  by  Bessel  functions  of  order  zero.  The  change  of  variable 

(2)  2 = x"’^y, 
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transforms  (1)  to  tlie  standard  form 
(3)  r-o- 

This  equation  is  of  the  form  3.3(1),  and  in  3.5(1)  we  may  take  a = 0 and 
a,  /3  arbitrary. 

We  obtain  formal  solutions  as  in  sec.  3.2;  equations  (5),  (6),  (7)  of  that 
section  become 

m*+l  = 0,  p = 0,  = (n  — 1/2)^  c^_,; 


and  with  the  £ibbreviation 


n 

^4)  “„=  n 


(1--1/2)*  [r(n  + l/2)]' 


v=  1 


2"n!  n 


and  appropriate  choices  of  c.,  we  obtain  two  formal  solutions 


2 a (-»•)"  X 
n=  0 


n „~n~X 


(5)  S,(x)-(2/;r)>^ 

(6)  S,(x)  = (2/n)^ 
of  (1).  Between  the  formal  scries  (5)  and  (6)  the  identities 


eo 

2 a i"  X 
n=  0 " 


n „~n-X 


(7)  S,(xe^^)«-S,(xX 

S^(xe^‘)  = S,(x), 


S,(xe~'^^)  = S^M 
S2(*e“^*)  = - S,  (x) 


hold. 

Since  <u  = ± i,  the  critical  line  is  the  real  x axis,  and  according  to 
the  theory  outlined  in  the  preceding  section,  every  solution  of  (1)  is  re- 
presented asymptotically  by  a linear  combination  of  S,  and  in  any 
sector  which  is  entirely  in  the  upper  (or  lower)  half-plane.  As  the  real 
axis  is  crossed,  the  coefficients  may  change.  We  shall  See  that  such 
changes  actually  occur. 

It  can  be  verified  by  substitution  that  (1)  is  satisfied  by  the  Bessel 
function  of  the  first  kind  of  order  zero 
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(8) 


« SS  0 


(m!)’ 


which  is  an  even  entire  function  of  x.  Poisson’s  integral  representation 

(9)  = du 

7T  “1 

may  be  verified  by  expanding  the  exponential  function  and  integrating 
term-by-tcrm . 

Let  us  assume  for  the  moment  that  Re  x > 0.  Then  e““  vanishes 
exponentially  as  Im  u ->  + oo,  and  we  may  break  up  the  integral  in  (9) 
according  to 

+J-,  ' 

In  the  first  integral  we  put 

u •=  1 + it,  1 - u = 1 + li  = 2 + it, 

and  in  the  second, 

— 1 + it,  1 - u = 2 — it,  1 + u = te 
thus  obtaining  two  functions  which  are  constant  multiples  of 

(10)  //'"(*)  = - e-**t"^  (2+it)-^  dt 

7T  0 


// (x)  = — e 

TT 


-ix+i7r/4  J-«>  {2-ity'^  dt. 

0 


The  functions  defined  by  (10)  are  known  as  Bessel  functions  of  the  third 
kind,  or  Hankel  functions , of  order  zero.  These  functions  are  defined  by 

(10)  for  Re  x > 0,  but  their  domains  of  definition  can  be  extended  to 

— n<  arg  X < 277  in  the  case  of  ffg'*,  and  to  —277  < arg  x < 77  in  the  case 
of  by  rotating  the  path  of  integration  as  in  sec.  2.2. 

It  can  be  shown  that  f/g”  and  are  also  solutions  of  (1).  Clearly 

(11)  + ^ 
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and  a closer  investigation  reveals  tliat  both  and  have  logarithmic 
singularities  at  the  origin.  The  knowledge  of  these  singularities  leads  to 
a definition  of  the  Hankel  functions  for  all  values  of  arg  *. 

The  integrals  representing  flankel  functions  are  Laplace  integrals, 
and  their  asymptotic  expansions  for  large  x may  be  obtained  by  means  of 
the  last  theorem  in  sec.  2.2,  the  result  being 

(12)  (*),  uniformly  in  arg  x,  as  a:  -*  oo, 

— w + f < arg  X <2n  — e,  f>0 

(13)  /fg**  (.x)  ~ Sj  (x),  uniformly  in  arg  X,  as  X -»  oo, 

— 2 n + f < arg  x <n  — e,  e > 0. 

From  (11)  we  then  have 

(14)  2£7g  (x)  ~ S , (x)  + Sj  (x),  uniformly  in  arg  X,  as  x -»  oo, 

— tT+e<  arg  x <n  — e,  e > 0. 

In  the  last  equation  we  have  used  the -symbol  ~ in  a somewhat  unusual 
manner,  in  that  the  right-hand  side  is  not  one  asymptotic  expansion  but 
the  sum  of  two.  A justificiation  of  this  use  is  based  on  the  circumstance 
that  on  the  real  axis  the  two  expansions  are  of  the  same  o''der  and  may 
be  com.bined  into  a single  expansion  (which  then  is  not  an  asymptotic 
power  series),  while  in  the  upper  [lower]  half-plane  S,(x)  [S2(x)]  is  re- 
cessive and  may  be  omitted. 

Vie  have  thus  obtained  the  asymptotic  expansion  of  <7^  (x)  in  the  whole 
plane  with  the  exception  of  a neirrow  sector  around  the  negative  real 
axis.  To  obtain  asymptotic  expansions  valid  in  sectors  including  the 
negative  real  axis,  w.^  remark  that  it  follows  from  (8)  and  (14)  that 

2c7g(x)  = 2Jg  (xe^‘)  ~ S,  (xe'^‘)  S^(xe^') 

as  X -*  oo  and  — n’  + e < arg(xe^*)  < w — so  that 

(15)  2c7g(x)  ~ S , (x)  - Sj  (x),  uniformly  in  arg  X,  as  x -»  oo, 

— 2 rr  + e < arg  x < - t 

by  (7);  and  similarly  J^{x)=  J^{xe~^^)  and 

(16)  2 Jg(x)  ~ - S , (x) Sj  (x),  uniformly  in  arg  x,  as  x oo, 

€ < arg  X < 2 77-  — f . 
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A comparison  of  (14),  (15),  and  (16)  shows  the  Stokes’  phenomenon. 
The  rays  excluded  by  narrow  sect<m  are  the  Stokes  rays.  At  first  it  may 
seem  strange  that  the  sectors  of  validity  of  these  asymptotic  expansions 
overlap,  but  there  is  no  contradiction  involved  in  this.  The  regions  of 
validity -of  (14)  and  (15)  have  the  common  part  —jT  + f<  arg  x <—(;  in 
this  common  part  5^  (x)  is  recessive  so  that  the  right-hand  sides  of  (14) 
and  (15)  are  asymptotically  equal.  Thus,  the  coefficients  of  the  formal 
series  jump  in  sectors  where  these  series  are  dominated  by  the  other 
series . 
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CHAPTER  IV 

DIFFERENTIAL  EQUATIONS  \SITH  A LARGE  PARAMETER 

In  this  chapter  we  describe  briefly  the  asymptotic  theory  of  ordinary 
horrogeneous  linear  differential  equations  of  the  second  order  containing 
a large  paran'eter.  For  this  theory  and  its  extensions  to  equations  of 
higher  orders,  and  to  systems  of  differential  equations  of  the  first  order, 
see  Ince  (1927,  p.  270(1.),  Kan  ke  (1944,  p.  62  IT.,  102(1.,  137  (f.,  213(1.), 
V^asow  (1953),  the  references  given  in  these  works,  and  the  references 
given  at  the  end  of  this  chapter. 

As  in  sec.  3.2,  we  may  transform  the  differential  equation  to  standard 
form 

y " + q{x,  A)y  = 0, 

where  x is  a real  or  complex  variable,  and  A is  a real  or  complex  pa- 
rameter. We  shall  investigate  the  behavior  of  the  solutions  of  this  differ- 
ential equation  as  A ->  A^,  and  without  loss  of  generality  we  take  A^  = 
The  reader  will  be  expected  to  be  familiar  with  the  basic  theorems 
regarding  the  dependence  of  solutions  of  a differential  equation  on 
parameters  occurring  in  the  equation  (see,  for  instance,  Kamke,  1930, 
sec.  17,  sec.  8). 

4.1.  Lioaville's  problem 

In  the  course  of  his  classical  investigations  of  the  Sturm-Liouville 
problem,  Liouville  discussed  the  behavior  of  solutions  of  the  differential 
equation 

(1)  y"  +[A*  p(x)+ r(*)]y  =0 

as  A -»  oo.  Here  * is  a real  variable,  o < x < 6,  p (*)  is  positive  and  twice 
continuously  differentiable,  and  r(x)  is  continuous,  for  c < x < 6.  Liou- 
ville’s  procedure  may  be  summarized  as  follows. 
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New  variables,  ^ and  tj,  are  introduced  by  the  substitution 

(2)  ^-f[p(x)]^dx,  »7-[pU)]V, 

which  carries  the  interval  a < x < b into  a < ^ < /3,  and  the  differential 
equation  (1)  into 

d^T] 


(3) 


d^ 


+ X‘r)  «=  p(^)il. 


where 


Ip;'  5 p ' " 
4 p- 


(4)  p(a=T^  -- 


16  p ■ 


is  a continuous  function  of  a < ^<  /S.  (p ' = dp/dx,  etc.) 

By  a method  similar  ^o  that  employed  in  sec.  3.3  it  can  now  be  shown 
that  solutions  of  (3)  satisfy  the  Volteira  integral  equation 

(5)  77(^)  = c,  cos  A^+ Cj  sin  A^+ A”'  sin  A{^- t)  p (t)  t? (t)  dt, 

where  a<  y<  /3  and  Cj,  Cj  are  arbitrary.  >7(^)  and  Cj  cos  A^+  sin  A^ 
have  the  same  value,  and  the  same  derivative,  at  ^ - y. 

The  solution  of  (5)  can  be  obtained  by  successive  approximations  in 
the  form 

(6)  77(^,  A)=  2 77  (^,  A), 

n=  0 

where 


77p(^,  A)=  c,  cosA^+  Cj  sinA^ 

^„+,  (^.  X)  = A“*  /^^.sinA(^-  t)  pit)  rj^it.  A)  dt 

If  lp(^)l  < >^»  it  is  easy  to  prove  by  induction  that 
lc,l  + lcj  ^'•l^-y!" 


71  — Oj  ly  •••  • 


A)|  < 


n ! 


A" 


n = 1,  2,  ...  , 


80 


ASYMPTOTIC  EXPANSIONS 


4.1 


.ind  in  the  case  of  a finite  interval  (a,  /S)  it  follows  that  (6)  is  uniformly 
convergent  for  a < ^ < /3,  A > A,>  0,  and  is  also  an  asymptotic  expansion 
of  Tj{^,  A)  as  A “♦  oa. 

The  T] ^ are  very  difficult  to  compute.  Other  approximations  for  large  A 
may  be  obtained  from  fonnal  solutions  (which  are  divergent  in  general). 
There  are  two  methods.  One  uses  formal  expansions 


S a (x)  A~"  cosA^+  2 /3  (x)  A~"  sinAf 

n=o"  n=o” 

for  y(x.  A),  while  the  other  uses  formal  expansions 

tA^+  2 y (x)A~" 

n=  0 " 

for  log  y (x,  A).  In  the  second  method,  y is  a solution  without  zeros.  In 
either  method  an  approximation  to  y.  is  constructed  by  taking  a partial 
sum  of  the  formal  expansion;  and  this  approxin-ation  is  compared  with  y 
by  means  of  an  integral  equation.  Either  of  these  two  methods  reduces  to 
Liouville’s  process  if  the  partial  sum  in  question  consists  of  one  term. 

4.2.  Formal  solutions 

Instead  of  4.1  (1),  we  shall  discuss  the  slightly  more  general  differ- 
ential equation 

(1)  y " + qix,  A)y  = 0. 

If  q{x.  A)  is  a formal  power  series  in  A ' tvith  coefficients  which  depend 
on  X,  then  two  linearly  independent  solutions  of  (1)  may  also  be  repre- 
sented by  formal  power  series  in  A ' . Cn  the  other  band,  if  the  formal 
expansion  of  q in  powers  of  A contains  positive  powers  of  A,  then  the 
formal  expansion  of  y will  be  a Laurent  series.  Nevertheless,  we  shall 
see  that  in  the  case  that  ^(x.  A),  as  a function  of  A,  has  a pole  at  A = oo, 
we  can  construct  formal  solutions  which  are  analogous  to  the  normal  and 
subnormal  solutions  of  sec.  3.2. 

In  (1),  we  shall  assume  that  q (x.  A)  is  of  the  form 

(2)  2 g„(*)A“-", 

0 
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where  the  q ^ (*)  are  independent  of  A,  and  A is  a positive  integer.  shall 
also  assume  that  q ^ (x)  does  not  vanish  in  the  interval  (or  simply- 
connected  region  in  the  complex  plane)  over  which  x varies. 

Corresponding  to  the  two  methods  mentioned  at  the  end  of  sec.  4.1, 
we  shall  obtain  two  kinds  of  formal  solutions  of  (1).  The  first  one  of 
these  is  of  the  form 

(3)  I a (x)A-"exp[‘s  /3^(x)  A*~^]. 

n=  O " I/=  0 

In  substituting  (3)  in  (1),  we  use  the  convention  q ^ = Q,  = 0 for 
n = -1,  -2,  -3,  ...  , and  /S^  = 0 for  v = -1,  -2,  ...  and  also  for  i/= 

+ 1,  ...  . 

All  summations  may  then  he  extended  over  all  integers,  and  we  obtain 

E K * C fil  S a,  A-"  2 X A*-  X A-" 

+ X a"  A""  + X <7  A**'"  X a A'"  = 0. 

n ^ n n 

Picking  out  the  coefficient  of  A^~" 

m V m 

+ 22  a'  /3'_.  +a''  =0 

n*“«  « n^2k 

m 

for  all  integer  values  of  ra. 

The  first  non-vacuous  condition  arises  when  n = 0,  If  we  set  n = 0,  1, 
...  , A:  — 1 in  (4)  we  obtain 

+2 /3;,/3;_^  = 0 m = 0,  1,  ...  , A:- 1 

V 

or 

(5)  /3;^+?,  = o 

(6)  2/S;  /s:  +g.  + “2*  m = 1,  ...  , i - 1. 

In  consequence  of  these  relations,  we  may  restrict  summation  to  m > A:  in 
the  first  sum  in  (4).  For  n = h in  (4)  we  have 

(7)  2a;  /s;  + a ,(/3;  + 9,  + /s;. j = o. 
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and  when  we  replace  ra  by  * + a in  (4), 


(8)  + 


v=  1 


«=  1 */=  m +1 


+ 2 S a'  /S'  4-  a"  , = 0 

. n“«i  ~ a n“ii 

«ss  1 


a = 1,  2, 


We  thus  see  that  (3)  satisfies  (1)  formally,  provided  that  the  and 
the  /3j^  satisfy  (5)  to  (8).  In  these  equations,  empty  sums  (i,e.,  sums 
whose  upper  limit  of  summation  is  less  than  the  lower  limit)  are  inter- 
preted as  zero.  Since  ^ 0,  we  may  choose  a branch  of  [— and 
then  (5)  determines  up  to  an  additive  constant.  Moreover,  /3J  0, 

and  hence  (6)  determines  /S^,  ...  , recurrently,  up  to  an  additive 

constant  in  each.  (7)  determines  up  to  a constant  factor,  and  (8) 
determines  a^,  ...  recurrently,  up  to  an  additive  constant  multiple  of 

Cg  in  each.  Corresponding  to  the  two  branches  of  i—q^J^r  we  obtain  two 
formal  solutions  of  the  form  (3)« 

A second  type  of  formal  solution  is 

(9)  exp[  1 /3„U)  A*'**]. 

n=  0 

Substituting  (9)  in  (1)  we  obtain 

2 /8"  A*”'*  + (2  /3 ' A‘~")*  + 2 q„  A**’"  = 0, 

and  comparing  coefficients  of  A^”", 

(10)  /3'*  + ?„  = 0 

(11)  2/3^/3>q„+"2  /s:^:_,  =0  n = l,...,A:-l 

1 

(12)  2/3;/3>q„+"2  /3:  /3;_,  +/s;;_^  = 0 n = k,  k + I,  ...  . 

»*  1 


There  are  two  linearly  independent  formal  solutions  of  this  type. 
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The  connection  between  these  two  types  of  formal  solutions  is  fairly 
obvious.  Equations  (10)  and  (11)  are  identical  with  (5)  and  (6),  and 
2 X~"  is  the  formal  expansion  of 

exp(  2 /S„A*-"). 

ns:  k 

Throughout  this  discussion  we  have  assumed  that  q(x.  A),  as  a func- 
tion of  A,  has  a pole  of  even  order  at  A = oo.  If  the  pole  is  of  odd  order, 
then  no  solution  of  the  form  (3)  or  (9)  exists,  and  instead  of  powers  of 
A we  must  expand  in  powers  of  A'^. 

4.3.  Assrinptotic  solutions 

We  shall  now  show  that  under  certain  assumptions,  the  differential 
equation  4.2(1)  possesses  a fundamental  system  of  solutions  which  are 
represented  asymptotically  by  the  formal  solutions  obtained  in  the  pre- 
ceding section.  It  does  not  matter  whether  we  compare  solutions  of 
4.2(1)  with 

^2  a (x)A'"exp[  2 

n=o  " v=o  ^ 

where  the  and  /3^  satisfy  4.2(5)  to  (8),  or  with 
exp[  5k  /8  (*)A*~"], 

n = 0 

where  the  /3^  satisfy  4.2(10)  to  (12),  for  the  a’s  and  /8’s  can  be  so  chosen 
that  the  ratio  of  these  two  expressions  is  1 + 0(A~^). 

We  fix  a positive  integer  N,  and  set 

ii+W-1  . 

(1)  y.(*)=exp[  2 /3„.(*)  A*""]  7=1,2, 

■'  n = 0 

where  /S^,  = — , and  for  each  /,  the  /3^^.  satisfy  4.2(10)  to  (12).  These 

coefficients  are  completely  determined  by  ...  , 92*+jy-i,  certain 
derivatives  of  these  functions,  and  we  shall  say  that  the  are  suffi- 
ciently often  differentiable  if  all  the  derivatives  entering  the  determin- 
ation of  the  /3^  . , n = 0,  ...  , 24  + /V  — 1,  ex'st  and  are  continuous 
functions  of  x.  We  let  x vary  over  a bounded  and  closed  interval  £ 
a < X < b,  and  A,  over  a sectorial  domain  S:  jAj  > A,,  < arg  A<  <f>\. 

The  theorem  to  be  proved  is  as  follows. 
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If  for  each  fixed  X in  S,  ?(*,  A)  is  a continuous  function  of  x overl-,  if 

(2)  q{x.  A)  = qAx)  A"*-"  + 0(A-^  ), 

n = 0 

uniformly  in  x and  arg  A,  as  A -►  <»  *n  S,  where  the  q are  sufficiently 
often  differentiable  in  I,  and 

(3)  Re  lA*  0 

when  A is  in  S and  x in  I,  then  the  differential-  equation 

(4)  y"  + q(.x.  A)  y = 0 

possesses  a fundamental  system  of  solutions,  y ,(*)  and  y ^{.x\  so  that 

(5)  y.(*)=  y.U)[l  + 0(A"^)], 

y '.(*)=  y;(*)  [l  + OCA-'^)], 

uniformly  in  x and  arg  X,  as  X-*  oo  in  S. 

We  shall  prove  this  theorem  by  a method  analogous  to  that  used  in 
sec.  3.3.  By  (3)  and  4.2(10)  we  may  choose  /S^,  and  so  that  for 
each  A in  S,  Re  [A*  (a;)]  is  an  increasing,  and  Re  [A*  ® 

decreasing,  function  of  x.  It  then  follows  from  (1)  that  for  each  sufficiently 
large  A in  S,  |y,(*)|  is  an  increasing,  and  jy^C*)!  a decreasing,  function 
of  X. 

To  establish  the  existence  and  the  asymptotic  property  of  y,(x),  we 
substitute 

(6)  y , (x)  = y , (*)  z (*) 
in  (4)  and  obtain 

y ' 

(7)  z"  + 2 — z'  +F(*,  A)z  = 0, 

^ I 

where 

Y'l  2*+w.-t  . 2fc+w-t  . , 

(8)  F(*,  A)  = -^ +<7=  % /3>*-"  + ( i /3„',  A*-")*  + 9 = 0(A-'") 

1-  /i=0  r»=0 


4.4  DIFFERENTIAL  EQUATIONS  WITH  A LARGE  PARAMETER  85 

uniformly  in  x and  arg  X,  as  A-*  «c  in  5,  by  (2)  and  4.2(10)  to  (12)>  Equa- 
tion (7)  may  be  written  as 

— \y]{.x)  + yf(x) F(*.  a)2=o, 

dx\_  dx  J 

and  by  two  successive  integrations,  and  a suitable  choice  of  tbe  con- 
stants of  integration,  we  obtain 

(9)  z{.x)=\- j‘K(x,t)F{t,k)z(t)dty 

<1 

where 

K{x,  t)^f‘yf(t)y;Hs)ds. 

Since  ]y,(a:)|  is  an  increasing  function,  we  have  |y,(t)|  < ly,(s)I,  and 
1^  (*,  t)l  < (5  — a) 

The  existence  of  z (x)  now  follows  from  the  general  theory  of  Volterra 
integral  equations,  or  can  be  established  by  successive  approximations. 
From  (8)  and  (9),  2(x)  « 1 + 0(A~^),  uniformly  in  x and  arg  A,  as  A-»  t» 
in  5.  Moreover,  2 (x)  is  differentiable, 

r'(x)  = -/*  yf(o  y;*(x)F(t,  A)2(0dt  = o(a-^), 

* I 

and 

y (x)  1 

2(x)+— j— 2'(x)  = y;(x)[i  + o(A"^)]. 

y,(*)  j 

This  proves  (5)  for  / = 1.  The  proof  for  / = 2 is  similar,  except  that  6, 
rather  than  a,  must  be  chosen  as  the  fixed  limit  of  integration  in  the 
integral  equation. 

4.4.  Application  to  Bessel  functions 

We  shall  now  apply  the  methods  of  the  last  two  sections  to  prove  the 
asymptotic  formulas* 

(1)  c^^(A  sech /3)  ~ (2jt  A tanh  ^ exp(A  tanh /3  — A/B)  as  A-*  i» 

\-X 

rrAtan/S  ) exp[±»  (Atan /3  — A/3— n/4)]  as  A-»  oo. 

Hie  symbol  y3a»ed  in  tbeoe  foimulss  must  not  be  confused  with  the  function  to  be 

used  below. 


(2)  sec  /3) 
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Equation  (1)  holds  for  /3  > 0,  uniformly  in  (3  ifO</S,  </3</32<oo. 
Equation  (2)  holds  for  0 < /3  < 7r/2,  uniformly  in/8  if0<f</3</r/2— e; 
and  in  this  equation  (and  similarly  later  in  (18))  the  upper  sign  holds  for 
the  lower  sign  for  Both  results  may  be  derived  from  integral 

representations  of  Bessel  functions  by  means  of  the  method  of  steepest 
descents  (Watson,  1922,  sec.  8.4,  8.41). 

The  functions 

J^(Ax),  y;^(Ax),  ^J^'*(AxX 

/ are  solutions  of  the  differential  equation 

(3)  y"  +[A"-(A*- ;4)x-*]y-0. 

This  equation  is  of  the  form  4.3  (4)  with 

4 = 1,  q^{x)=\-x~^,  q ^{x)  = (2x)~^ , 

all  the  other  q^  (x)  vanishing  identically.  The  points  x = 0,  <»  are  singular 
points  of  (3),  and  x = 1 is  a so-called  transition  point  at  v/hich  the  con- 
dition 4.3(3)  is  violated  for  any  value  of  A.  On  any  interval  a £ x 
which  does  not  include  any  singularity  or  transition  point,  the  theorem  of 
sec,  4.3  will  yield  the  general  form  of  the  asymptotic  solutions  but  it  will 
fail  to  indicate  the  expression  of  y,  and  yj  in  terms  of  the  standard 
Bessel  functions.  In  order  to  identify  our  solutions  in  terras  of  Bessel 
functions,  we  shall  extend  the  interval  to  one  of  the  singular  points  of 

(3) .  Since  this  case  is  not  covered  by  the  theorem  of  sec.  4.3,  we  shall 
use  the  m.ethods  rather  than  the  results  of  the  preceding  sections. 

Let  us  first  discuss  (3)  on  the  interval  0 < x < b <1.  From  4.2(5)  we 
have 

/3'(*)=  ±(x-"-l)«, 

and  hence 

(4)  ±/3o(x)  = /3(x)  = J (x"^  - l)’^  dx  = (1  - x*)^  + log  . 

1 + (1  - X I'’ 

From  4.2(7), 

2a'  /3'  + «o/3-  =0, 
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so  that 

(5)  (x).  [/S' (*)]-»«.  a (x)  -x^a-x*)-’^. 

With  a and  /3  so  defined  we  form  the  functions 

(6)  y,(x)- a(x)c^H  y2(x)=  a(x)e-^‘*\ 

which  correspond  to  the  leading  terms  of  the  formal  solutions. 
The  integral  equation  for  z = y,/y , is 

(7)  z (x)  = 1 - K{.x,t)F  it.  A)  z (t)  dt. 


By  a straightforward  computation  we  find 

y"  / 1 \ 1 

(8)  /r(x,A)— ^ + A*  ( r-—  j + — 
y,  \ X / 4x 


4 + x‘ 


4(1 -X*)*  ’ 

so  that  F is  bounded  on  the  interval.  Moreover,  a~*  *=  jS'  from  (5),  and 


(9)  K{x, 


r * 

r ait)  1 

f 

t 

L J 

expl2A[iS(t)-  ^(s)]}  ds 


aHt) 

2A 


(1  - expl2A[/8(t)  - /8(x)]l ). 


Now,  /3  (x)  is  an  increasing  function  of  x,  /8(t)— j8(x)<0  for  0<r<x<h<l, 
and  the  exponential  function  will  be  bounded  if  Re  A > 0.  Also  a^(t)  is 
bounded.  We  thus  arrive  at  the  estimate 

(10)  |F(t,  A)A(x,  t)|  < 


for  ReA>0,  0<t<x<6<l.  Here  C is  independent  of  A,  xj  t. 

We  are  now  ready  to  solve  the  integral  equation  (7)  by  successive 
approximations  in  the  form 

2 „ + , (*)  = / * ^ ^ ^ 

2 (x)  = S 2 (*). 

n 


71  — Oj 
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From  (10)  it  is  easy  to  prove  by  induction  that 


so  that  the  series  defining  z converges  uniformly  in  * and  A if  A is 
bounded  away  from  zero;  z (x)  satisfies  the  integral  equation,  is  twice 
continuously  differentiable,  satisfies  the  appropriate  differential  equation, 
and 

2 (x)  = 1 + 0 

This  establishes  the  existence  of  the  solution  y,  of  (3)  which  has  the 
property 

(11)  y,(x)=  y,(x) 

It  remains  to  show  that  this  solution  is  a multiple  of  J;,(AX). 

Since  y ^ is  a solution  of  (3), 

(12)  x~'^  y , (*)  = c , (A)  J^(Ax)  + C2(A)  y;^(A*). 


0 < X < 6 < 1,  Re  A > 0. 


Now,  fix  A,  and  let  x -♦  0.  It  is  well  known  (Watson,  1922,  p.  40,  p.  64) 
that 


J^\x)  ~ 
y;,(Ax)  ~ 


(Ax/2)^ 

r(A+  1) 

(Ax/2)^ 
r(A+  1) 


as  X -*  0 


^ (Ax/2)-^ 

cot  A n’ cosec  A n 

r(i- A) 


as  X -»  0, 


and  it  follows  from  (11),  (6),  (5),  and  (4)  that  * 

x“^  y , (x)  ~ x~^‘  Y j (x)  ~ ~ (x/2)^  as  x -*  0. 

Making  x ->  0 in  (12),  we  see  that 

c,(A)  = e^A“^r(A+ 1).  c^(A)=0, 


Note  that  is  the  approximate  solution,  * Bessel  functioa  of  the  second  kind. 
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and  by  Stirling’s  formula. 


,(A)=  (2  7rA)’<  + 


ViC  thus  find  that 


= (2  jta*)  y , (*) 


when  0 < a:  < h < 1 and  Re  A > 0.  If  we  put  * = sech  /S  and  take  A positive 
in  the  latter  form,  we  obtain  (1). 

Let  us  now  turn  to  the  discussion  of  (3)  on  the  interval  1 < a < x < oo. 
In  this  case 


(14)  aix)  = {x^  - 1)-^  = [-t/3'  (x)]-^ 

/3(x)=t/  (1-x**)^  dx  = i{x^-  1)'^  -i  cos"’  r"’. 

where  cos"’  denotes  the  principal  value  of  the  inverse  cosine, in  particular 
cos"’  x"’  -♦  rr/2  as  x ->  w.  The  comparison  functions  are  again  of  the 
form  (6),  with  a and  /3  defined  in(14)«  The  integral  equation  for  z = y ^/Y  ^ 
is 

(15)  2 (x)  = 1 + K (x,  i)  F (t.  A)  2 (t)  dt. 


Equation  (8)  holds,  and  shows  that  F {t.  A)  = 0{t~^)  for  I > h > 1 and  all 
A.  Since  a * = — /S  * in  this  case,  the  evaluation  of  the  integral  in  (9) 
leads  to 


K{x,  f)  = - 


2A 


(1  - expl2A[/3(t)  - /S  (x)]  } ). 


l\'ow,  —i  /3  is  an  increasing  function  of  x,  and  t > x,  so  that  the  exponential 
function  will  be  bounded  if  Im  A > 0.  Also  a^{t)  is  bounded,  and  we  have 
the  estimate 
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(16)  |F{t,  X)K(x,  0|  <— ^ 


for  ImA>0,  l<a<*<«,  C being  independent  of  A,  x,  t. 
i(e  now  set 

+ , (*)  = /C  (*,  t)  F it.  A)  z^it)  dt 


zix)=  S zixX 

n=  0 

From  (16)  it  is  easy  to  prove  by  induction  that 


n = 0,  1,  ... 


" "I  vw*y 


so  that  the  series  defining  z converges  uniformly  if  A is  bounded  away 
from  zero;  z ix)  satisfies  the  integral  equation,  is  twice  continuously 
differentiable,  satisfies  the  appropriate  differential  equation,  and 


(x) = 1 + 


This  establishes  the  existence  of  a solution  of  (3)  which  has  the 
property 


(17)  y,(x)=  y,(*)  [■■”©]  1 < a < X < oo,  Im  A > 0. 

It  remains  to  show  that  this  solution  is  a multiple  of  Hj^\Kx\ 

Since  y , is  a solution  of  (3), 

(18)  x~^  y j (*)  = c , (A)  H^^Xx)  + CjfA)  fij^KXx). 

Let  us  fix  A and  make  x -»  «>.  It  is  well  known  (Watson,  1922,  sec.  7.2) 
that  for  0 < arg  X.  < n 

^ - x.i>  ^ X— • 


^ L V 2 vJ 


as  X oo^ 
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I 


and  it  follows  from  (17)  and  (14)  that 

y , (*)  ~ y,  (*)  ~ exp[A/3  (*)]  ~ exp[t  A*  - t A V2]  as  * -♦  oo. 

Making  ac  -*  eo  in  (18)  we  see  that 

c , ( A)  - (1/2  »rA)^  e c ^ (A)  = 0. 

If  arg  A=  JT,  we  must  use  a slightly  different  asymptotic  formula  for 
but  the  conclusion  remains  the  same,  ^'e  thus  have 

(19)  ^j^'*(A*)  = (l/2»TA*)"^  y,(*)  j^l  +0 

where  1 < a < x < ec  and  Im  A > 0.  By  a similar  proof 

1.0 /-L 

V Ax 

when  1 < a < X < oo  and  Im  A < 0 . If  we  take  A positive  and  put  x = sec  /S, 
we  obtain  (2). 

441.  Transltioo  points 

Let  us  consider  again  Liouvillc’s  differential  equation 

(1)  y"  + [A*  p(x)  + r(x)]  y = 0 

with  large  positive  A.  As  in  sec.  4.1,  x is  a real  variable,  a < x < b,  p(x) 
is  real  and  twice  continuously  differentiable,  and  r (x)  is  continuous,  for 
a < X < b.  Instead  of  assuming  p (x)  to  be  positive,  we  now  assume  that 
p (x)  has  a zero  in  (a,  b),  A zero  of  p (x)  will  be  called  a transition  point 
of  the  differential  equation  (1).  For  the  sake  of  definiteness,  we  assume 
that  p (x)  has  a simple  zero  at  x = c,  and  no  other  zero  in  a < x < h,  and 
also  that  p ' (c)  > 0,  so  that  p (x)  < 0 when  a < x < c,  and  p (x)  > 0 when 
c <x  < b. 

We  have  seen  that  in  the  interval  c + f<x<6,  f>0,  where  p (x)  is 
positive,  solutions  of  (1)  are  asymptotically  of  the  form 

(2)  c^tpfx)]"*^  coslA  J [pCx)]'"^  dxl  + CjLpfx)]"*^  sin!  A J[p  (x)]^  dx\. 


(20) 


A/i«(Ax) 


(1/2»7Ax)~^  e^'^*  y^(x) 
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an«i  it  can  be  shown  similarly  that  in  the  interval  « < * < c — r , where  p (*) 
is  negative,  solutions  of  (1)  are  asymptotically  of  the  form 

(3)  explA/[-p  (*)]^  dx\  4 [-p  (*)]“’^expj-A  f[-p{x)]^dxl. 

The  validity  of  these  asymptotic  forms  depends  on  the  fact  that  p (^)  in 
4.1  (4)  is  a bounded  function,  and  hence  it  is  clear  that  neither  of  the 
asymptotic  forms  can  he  valid  at  * = c.  To  the  right  of  x » c,  (2)  shows 
that  every  solution  of  (1)  has  an  oscillatory  character;  and  to  the  left 
of  X = c,  (3)  shov's  that  every  solution  of  (2)  has  a monotonic  character. 

In  the  immediate  vicinity  of  x = c,  the  transition  takes  place  from  one 
type  of  behavior  to  the  other. 

In  this  situation  two  problems  arise.  The  first  of  these  is  the  problem 
of  finding  the  connection  between  the  constants  c,  and  on  the  one 
hand,  and  the  constants  c ^ and  c ^ on  the  other  hand,  if  (2)  and  (3)  re- 
present asymptotically  the  same  solution  of  (1)  in  different  intervals;  and 
the  second  problem  is  the  determination  of  the  asymptotic  form  of  the 
solutions  of  (1)  in  the  interval  (c  — r,  c + f). 

There  are  in  essence  two  ingenious  methods  for  solving  the  first 
problem.  The  first  of  these  methods  was  used  by  Jeffreys  in  1923,  and  re- 
discovered by  Kramers  a few  years  later.  It  is  based  on  the  remark  that 
sufficiently  near  to  x = c,p  (x)  may  be  approximated  by  the  linear  function 
(x  — c)  p ' (c),  and  r(x)  may  be  neglected.  The  resulting  differential 
equation  can  then  be  solved  in  terms  of  Bessel  functions  of  order  ± 1/3. 

A comparison  of  the  asymptotic  forms  of  Bessel  functions  with  (2)  and 
(3)  leads  to  the  desired  connection  formulas  between  c,  and  on  the 
one  hand,  and  c^  and  on  the  other  hand.  The  second  method  was 
developed  by  Zwaan  in  1929,  and  it  avoids  the  transition  point  altogether. 
ix  p (x)  and  r(x)  are  analytic  functions  of  x,  the  differential  equation  is 
integrated  along  a path  in  the  complex  plane  which  consists  of  the  real 
intervals  (a,  c — e)  and  (c  + <,  b),  and  of  a semi-circle  in  the  complex 
plane  joining  c — ( and  c 4-  e.  Along  this  path,  p of  4.1(4)  is  bounded  and 
Liouville’s  method  (or  a variant  of  it)  can  be  applied.  This  method  leads 
to  the  same  connection  formulas  as  the  first  one.  Both  methods  can  be 
extended  to  cases  when  p (x)  has  a zero  of  an  arbitrary  order.  They  are 
known  as  the  WKB  method,  or  sometimes  also  the  KKBJ  method.  For  a 
fuller  descrip^tion  see,  for  instance,  Morse  and  Feshbach  (1953,  p. 10920.). 
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The  second  problem,  the  determination  of  asymptotic  forms  of  the 
solutions  of  (1)  in  the  interval  (c  — e,  c + eX  is  much  more  difficult. 
There  seems  to  be  no  simple  elementary  function  which  describes  the 
transition  from  monotonic  to  oscillatory  behavior,  and  it  seems  plausible 
that  the  asymptotic  forms  willinvolve  some  higher  transcendental  function. 
Now,  the  simplest  differential  equation  of  the  form  (1)  and  having  a 
transition  point  is  the  equation 


(4) 


d^y 

~dP 


+ A*  *y  = 0; 


the  solutions  of  this  equation  are  fairly  well  known,  and  it  seems  tempt- 
ing to  seek  asymptotic  forms  of  the  solutions  of  (1)  in  terms  of  solutions 
of  (4).  Liouville’s  transformation  4.1(2)  transforms  4.1(1)  into  a differ- 
ential equation  with  approximately  constant  coefficients,  and  similarly, 
we  can  find  a transformation 

(5)  T/=0(*)y 

which  transforms  (1)  into  an  equation  approximately  of  the  form  (4).  Ihe 
transformation  (5)  carries  (1)  into 


d^T] 


j_  „2 — ^ + r I 
<p'  0/  ^ L 


1 2 ,-1 


d^  x{j' 


»7  = 0. 


In  order  to  reduce  this  differential  equation  approximately  to  the  form 
(4),  we  first  determine  \{j  so  that 


(6) 


r 


±L 


and  then  0 so  that 


= 0, 


0 = 


(7) 


(2  ~ 


<;60'^=p. 


With  and  i/r  so  determined,  the  differential  equation  becomes 


(8) 


A*  = p(^)t;. 
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where 


(9) 


1 


3 9S"*  I- 

4"  0^  “7^  ' 


Under  the  assumptions  on  p and  r made  at  the  beginning  of  this 
section,  there  is  a unique  three  times  continuously  differentiable  real 
function  0 which  satisfies  (7).  For  this  function,  tp'  is  bounded  away 
from  zero,  p(^)  is  i bounded  function,  and  we  shall  expect  that  the 
asymptotic  form  of  the  solutions  of  (8)  is 

(10)  c, 

where  H ^ (x)  and  H (x)  are  two  linearly  independent  solutions  of  (4). 

This  generalization  of  Liouville’s  method  was  originally  developed 
for  the  purpose  of  obtaining  asymptotic  forms  of  the  solutions  of  (1)  in 
the  interval  {c  — e,  c + e\  but  it  is  clear  that  it  can  be  extended  to  the 
entire  interval  (a,  b).  The  extension  enables  us  to  dispense  with  three 
different  asymptotic  forms  in  (a,  c — (c  — e,  c + e),  (c  + e,  b)  respec- 
tively, and  yields  a single  uniform  asymptotic  representation  of  the 
solutions  of  (1)  in  o < X < 6.  The  method  was  originated  by  I. anger,  who 
developed  it  in  a number  of  memoirs  of  which  a few  are  listed  at  the  end 
of  this  chapter.  Among  those  who  developed  further  Langer’s  method  we 
mention  in  particular  Cherry.  A survey  of  the  literature  regarding  this 
method  is  available.  (See  reference  at  the  end  of  this  chapter.) 

Before  describing  this  method  in  greater  detail  we  shall  list  briefly 
some  properties  of  the  solutions  of  (4). 

4.6.  Airy  functions 

'Ihe  differential  equation 
d^w 

(1)  “T-; zm  = 0 

dz 

can  be  reduced  to  the  differential  equation  satisfied  by  Bessel  functions 
of  order  1/3  (Watson,  1522,  sec.  6.4).  We  shall  use  the  notation 

(2) 

o 


ct)  = e 


27ri/3 
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Two  linearly  independent  solutions  of  (1)  are  the  80-calI»,d  Airy 
functions  of  the  first  and  second  kind, 

and  by  direct  computation 

(4)  Ai(z)  Bi ' (2)  - At' (z)  Biiz)  = rr 


We  also  have 

(5)  Ai(-z)- 


B 


(0-c7,/3(^)]. 


For  real  *,  the  integral  representations 


(6)  A i 


f \ 

cos  ( — + xt  1 dt 

r f ^ ^ 1 

i{,x)=  n I I exp  I ~ y + xt  J 


dt 


hold,  and  these  integrals  may  be  converted  into  contour  integrals  which 
remain  valid  when  x becomes  com.plex, 

Ai  (2)  and  fit  (2)  are  both  entire  functions  of  z:  they  are  real  for  real  z. 
Ai  (2)  has  a string  of  zeros  on  the  negative  real  axis,  and  it  has  no  zeros 
elsewhere  (Watson,  1922,  sec.  15.7).  For  any  integer  m,  w ^{z)  = Ai{co"z) 
is  also  a solution  of  (1).  By  direct  computation. 


(7) 


1 

m;_  (2)  10  ' (z)  - w‘^{z)w  ifz)^ — 


n{m—k) 

sin exp 


/m  + k 1 \ 
ini + — ) 

L \ 3 2yj 


3 
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and  hence  it  is  seen  that  and  are  linearly  independent.  On  the 

other  hand,  any  three  of  the  w ^ are  linearly  dependent,  in  particular, 

(8)  (i)  + <y  +,  (z)  + ^^(z)  = 0. 

Also 

(9)  Bi  (z)  = * [co^  Ai z)  — co  d* (ojz)]. 

Ihe  power  series  expansions  of  Ai(z)  and  Biiz)  follow  from  (3).  In 
particular 

(10)  3^’  /li(0)=  3'^«  5,  (0)  = —-^ 

- 3 di'  (0)  = 3"'/«  Bi ' (0^ 

The  asymptotic  behavior  of  di(z)  has  been  investigated  in  sec.  2.6 
for  -n/3  < arg  z < Tr/3.  The  analysis  can  be  extended  by  rotating  the 
path  of  integration  in  the  t-plane  of  2.6(3)> 

(11)  di(z)  = — z"*^  e"^  [1  + 0(^"’)1 

as  z -»  oo,  — jT  < arg  z < rr. 

Results  valid  in  sectors  containingthe  negative  real  axis  may  be  obtained 
by  using  (8). 

(12)  di(z)  = i z-X  le-^[l  + 0(C‘’)]  + «4l  + 0(C”)]l 

2 

as  z -►  e»,  n/3  < arg  z < 5 jr/3 

di(z)  -—  n-''^  z-^  le'^  [l  + 0(C"')l-»e^[l  + 0(C"’)]l 
2 

as  z -*  on,  -*5  7r/3  < arg  z < — ir/S, 
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From  (9), 

(13)  fii(2)=  + 

as  2 -*  oo,  — ;r/3  < arg  z < n/3. 

All  these  asymptotic  representations  hold  uniformly  in  arg  2 if  2 is 
restricted  to  a closed  sector  inside  the  open  sector  stated  above. 

It  follows  from  the  asymptotic  formulas  that  the  functions 

(14)  (1  + |2|  ^)  e ^ Ai  (2)  —JT  < arg  2 < w 

(15)  (1  + |2|*^)  e“^  Bi(2)  - 7t/3  < arg  2 < 7r/3 


are  bounded  functions  of  2.  The  reciprocals  of  these  functions  are  also 
bounded  provided  that  the  zeros  of  the  first  function  are  suitably  ex- 
cluded, for  instance,  if  — n + f < arg  2 < 77  — r,  t > 0,  in  (14):  the  second 
function  has  no  zeros. 

For  further  information  on  Airy  fimctions,  and  for  numerical  tables  of 
these  functions,  see  Miller  (1946). 

We  conclude  this  section  by  proving  the  following  inequality.  If  w (z) 
is  that  solution  of  (1)  which  satisfies  u)  (t)  = 0,  u; ' (t)  = 1,  then 


(16)  0 < 


2{xt)^  w{x) 


< 1 


for 


* > 0,  t>0,  f = 


,3/2 


>0, 


3 3 

The  proof  depends  on  a comparison  theorem  of  Sturm’s  type. 

ei-r_g'r-i 


v(*)  = 


2(,xt)^ 


satisfies  the  differential  equation 
/ 5 \ 
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and  the  same  initial  conditions  at  t as  tv.  From  (1)  and  (17)  it  follows 
that  both  V (x)  and  w (x)  are  0 if  x > 0,  t > 0,  x t.  Consider  the  function 
f(x)  = w(x)/v(x),  X ^ t..  Clearly,  /(*)  -»  1 as  * -»  f,  so  we  put  /(t)  = ].. 
Also,  V and  w have  the  same  sign  when  x t. 

d 5 txw 

(t  t \ St  //  ^ f\ 

w V — WV  ) - W V wv  = — ■ < 0 

dx  16  X 

and  w'  V — wv'is  a decreasing  function  ofa:.  Since  this  function  vanishes 
at  X = t,  it  follows  that  f (x)  < 0 if  x > t,  and  f {x)  > 0 when  0 < a:  < t. 
'Ihus,  f (x)  has  a maximum  at  x = t,  and  0 < f{x)  < 1 for  0 < x,  f < oo, 

4.7.  Asymptotic  solutions  valid  in  the  transition  region 

We  assume  that  (a,  6)  is  a bounded  interval,  a < c < b,  p (x)  is  real  and 
and  twice  continuously  differentiable  for  o <x  < b,  p\x)  < 0 for  a < x < c, 
p (x)  > 0 for  c < X < 6,  p (c)  = 0,  p ' (c)  ^ 0,  r (x,  X)  is  a bounded  function 
of  X and  A when  a < x < b and  A varies  over  a sectorial  domain  S:  |A|  > A,, 
0Q  < arg  A < and  for  each  fixed  A in  5,  r(x,A)  is  a continuous  function 
of  X for  a < X < 6.  We  then  consider  the  differential  equation 

(1)  y " + [A^  p (x)  + r(x.  A)]  y = 0, 

and  introduce  0(x)  as  the  (unique)  continuously  differentiable  real  solu- 
tion of  4 >5  (7),  i.e. 

(2)  — [</>(*)]^^  = /(*)=  dt  X > c 

3 ^ 


— [-<^U)]^"  = f{x)  = / ' [-p(t)]’^"  dt 
3 * 


X < c. 


Here  ^(x)  > 0 when  x > c,  0(x)  < 0 when  x < c,  and  all  fractional  powers 
have  their  positive  values. 

The  functions 


(3)  y , (*)  = [<^ ' (*)]"^  Ai  [-  A"^"  e 0 (x)] 

YAx)  = [.^'(x)]-''"  BiC-A"^  <f>{x)] 


m = 0,  ±1 
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satisfy  die  diffsrsotial  eqnstion 

(4)  y-(*)+  \*p(*)+y  1^*1  J 

I ® V 


where 


is  the  Schwarzian  derivative  This  differential  equation  follows  from 

4.6(1)  and  4.5(7).  Since  the  differential  equations  (1)  and  (4)  differ  in 
terms  which  are  comparatively  small,  we  consider  the  y's  as  leading 
terms  of  formal  solntibns. 

Under  the  above  assumptions,  the  differential  equation  (1)  possesses 
solutions  which  are  represented  asymptotically,  in  appropriate  sectors 
of  the  complex  It-plane,  by  the  Y m « — 1,  0,  1,  2.  The  proof  is  similar 
to  the  proofs  given  in  sections  3.3,  3.4,  and  4.3;  and  it  will  be  conducted 
in  several  steps. 

Let  Yifc)  be  any  solution  of  (4)»  let  K(?e,  t)  be  that  solution  of  (4) 
which  satisfies  the  initial  conditions 

(5)  K (t,  l)  - 0,  — (t,  t)  - 1, 

dx 


and  let  o £Xg  ^ b',  then  the  solution  of  the  integral  equation 

(6)  y(x)-y(x)  + J*  /((x,  t)F(t,  A)y(t)dL 

*0 

where 

(7)  F(x,  A)  « ~ X } “ r(*,  AX 
satisfies  (1). 

Proof:  Equation  (6)  is  a Volterra  inte^al  equation  with  a continuous 
kernel.  For  any  fixed  value  of  A in  5,  the  existence  and  uniqueness  of 
the  solution  follows  from  the  general  theory  of  Volterra  integral  equations 
(orcan  be  established  as  in  sec. 4.3).  This  solution  is  twice  costinoously 
differentiable,  and  substitution  in  (1)  shows  that  y satisfies  the  differ- 
ential equation  (1). 
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For  o < X,  t < b and  all  A in  St 

(8)  l/CU,  f)F(f,  A)| 

^ exp  1 1 Rel[-A^^^ri,U)]=^^-[-A^3  0(t)]V2{| 

- |A|  (|Ar'^«+  |0(x)|'^^'viA!-'^®  + |<^(OI'"-)[<A'(*)<;6'(0]''" 

where  A is  independent  of  x,  t,  A.  All  fractional  powers  have  their  princi- 
pal values. 

The  proof  of  (8)  is  different  according  to  the  position  of  x and  t 
relative  to  c and  relative  to  each  other,  and  according  to  arg  A.  We  shall 
give  the  details  for  a < x,  f < c,  — 7t/2  < arg  A < tt/2.  Since 

Yl  Y^-Y^  Y'^  =-A*^  {Ai'Bi-Ai  5t')  = «■"’  A^^ 
by  (3)  and  4.6(4),  we  see  that 

(9)  K{x,  [y^(x)  Y^{t)~Y^{t)  y^C*)] 

satisfies  (4)  and  (5).  Now,  <f>  < 0,  and  |arg(-A^®  9!>)|  ^ n/Z,  so  that  the 
expressions  4.6(14)  and  4.6(15)  are  bounded,  and  there  is  a S inde- 
pendent of  X,  t,  and  A such  that 

flexp}-|Re[-A^=’ 

' (1+  |A|'^  ’ 

with  a similar  estimate  for  y^.  Applying  these  estimates  in  (9),  we  prove 
(8)  in  this  case. 

In  order  to  prove  (8)  in  any  other  case,  express  y^  and  y^,  and  hence 
K,  in  terms  of  two  linearly  independent  solutions  of  (4),  of  which  one  is 
bounded  as  A -*  oo  (this  was  the  solution  y^,  in  the  above  case),  and  then 
apply  estimates  derived  from  4.6(14),  (15). 

We  shall  now  establish  the  existence  of  solutions  y^,  m = — 1,  0,  1,  2, 
corresponding  to,  and  asymptotically  represented  by,  the  y , It  will  be 
necesseuy  to  impose,  in  each  case,  restrictions  on  arg  A.  These  re- 
strictions correspond  to  the  introduction  of  branch-cuts  in  the  A^-plane. 

For  y^j  and  y we  assume  Re  A > 0,  and  define  these  functions  as 
solutions  of  the  Volterra  integral  equations 
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(11)  y^{x)^Y^{x)  + J‘K{x,t)F{t,  k)y^U)dt. 

F o>-  any  fixed  A in  S,  the  existence  of  the  solutions  of  these  integral 
equations  follows  from  the  general  theorj’  of  V'olterra  integral  equations 
(or  uf  successive  approximations).  We  shall  prove  thaty^  y ^ ~ Y ^ 

as  A -►  oo.  Re  A > 0.  * 

If  Re  A > 0,  Im  A ^ 0,  and  a < * < fc,  or  if  A ^ and  o < x < c,  we  have 
I'd  ^2  5^  this  case  we  may  out  y^  = I'o  >"2  ==  ^2  ^2’ 

obtain  ir  ei^val  equations  for  and  r..  Since  under  these  circumstances 
4.6(14)  4.6(15)  are  bounded,  and  bounded  away  from  zero,  we  have 


ypU) 


_cf>'  (t)  _ 


iAr»^«  + \<^{x)\''* 
|A|“’^"+  |.;6(t)| 


X exp  Re  I [-A*^  - [-X^^  0 (*)]"""  , 

and  a similar  estimate  for  |y  j(t)/y2(x)|.  Combining  these  estimates  with 
(8),  we  obtain  inequalities, 

y (t)  r 

(12)  - 9 — Kix,  t)  Fit,  A)  < 5 

" |A|  <;6'(0  i0(«)l'^ 

o < C < X £ b and  ImAj^O,  or  a < t < x < c and  A > 0. 


YAt) 

-f—  K{x,t)F{t,  X)  < 
Y,(x) 


|A1  <^'(0  l«/>(t)l' 


a < X < b,  |t-c|<|x-c|,  and  Im  A 7^  0, 
or  a < X < t < c and  A > 0 

for  the  ke-rels  of  the  integral  equations  for  and  z For  each  fixed 
A,  the  imegral  equation  for  z m = 0,  2 has  a bounded  solution.  Let 
Z ^ (A)  be  the  maximum  of  \z  (x)|  for  a < x < b.  From  the  integral  equa- 
tion and  (12),  (13), 

c r b dt 

Z (A)<1  + — Z (A)  / r-. 

- - A - / lp(t)|*^ 


I * 
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For  sufficiently  large  |A|,  |Z_  (A)|  < 2,  and  from  (10),  ill),  (12),  (13), 

(14)  y^  (*)  = Y ^ (x)  [1  + 0(A~')i  uniformly  in  x,  as  A -*  «>,  m = 0,  2 

provided  that  Re  A>0,  a < x < b,  and  also  that  in  case  of  a real  A,  x is 
restricted  to  a < x < c. 

In  the  case  of  a positive  real  A,  and  Y ^ have  zeros  when  x > c,  and 
(14)  fails  to  hold  near  such  zeros.  An  appropriately  modified  result  can 
be  derived  from  (10)  and  (11)  by  using  (9)  and  estimates  of  Y ^ (t),  Y ^ (^)* 
From  (9) 

y^{x)=Y^{x)[l  + nk-^^  S^Y^it)  Fit,  X.)y^it)dt] 


Y^{x)J^^  Y^it)Fit,  k)y^it)dt. 


Here  the  integrals  can  be  broken  up. 


(14)  may  be  used  in  the  first  integral,  and  thus  it  can  be  proved 

(15)  yo(^)=  Y^ix)[l  + 0i\~')]  + Y^ix)0(X''), 

uniformly  in  x,  as  A -»  oo,  c < * < ^»  A > 0. 

Similar  is  the  proof  of 

(16)  y^ix)=Y^ix)  [1  +C)(A’')]+  yo(x)0(A'’), 

uniformly  in  x,  as  A -*  oo,  c < x < 6,  A > 0. 

Equations  (14)  to  (16)  describe  completely  the  asymptotic  properties  of 
Xo  and  y^. 

'Ihe  work  on  y^,  is  similem.  We  define  y^  and  y„,  by  the  integral 
equations 

(17)  y ^ (x)  = Y^  (x)  - Kix,  t)  Fit,  X)  y^(t)  dt  m=  1,-1, 


and  analyze  the  solutions  of  these  equations  in  a manner  similar  to  the 
analysis  of  (10),  (11),  assuming  Im  A > 0 in  the  case  of  y j,  and  Im  A < 0 
in  the  case  of  y_,*  Y , and  y_,  have  zeros  when  A is  imaginary  and  x < c, 
and  in  this  case  the  asymptotic  forms  need  some  modification.  The  final 
results  are 


4.7 


DIFFERENTIAL  EQUATIONS  WITH  A LARGE  PARAMETER 


103 


(18)  y,  (*)  •=  y,  (a.)  [1  + 0(A”')],  uniformly  in  *,  as  A-*  oo 

ImA^O,  ReA^O  or  c ^ x ^ b,  —i\>0 

(19)  y,(*)=  y,(a:)[l  + 0(A‘’)]+  y_,(x)0(A'’), 

uniformly  in  x,  as  A -*■  oo,  a < x < c,  — t A > 0 

(20)  y.,  (x)  = y„,  (*)  [1  + 0(A*')],  unifomily  in  as  \ oc> 

< 6,  Im  A < 0,  Re  A 0,  or  c < x < fe,  t A > 0 

(21)  y_,(x)-  y_,(*)  [1  + 0(A'’)]+  y,(a:)0(A‘’), 

uniformly  in  as  A -►  oo,  o < * ^ c,  i A > 0. 

Equations  (14)  to  (16)  and  (18)  to  (21)  establish  the  result  staU  J 
the  beginning  of  this  section.  By  a more  careful  discussion  of  the  ino-gral 
equations  it  is  possible  to  show  that  also.y^  ~ y|^ . 

If  r (*,  A)  can  be  expanded  in  powers  of  A~’,  then  it  is  poss! 
develop  formal  solutions  of  (1).  The  approximations  discussed  it.  ’nis 
section  appear  as  the  leading  terms  of  the  formal  solutions.  As  in  tuu 
case  discussed  in  the  earlier  parts  of  this  chapter,  there  are  two 
of  formal  solutions.  The  first  kind  corresponds  to  4.2(3),  and  it  is  of  the 
form 

(22)  y(x)  2 a^(x)  A~"+  y'(x)  2 /3„(x)  A""'’ , 

n=  0 0 

where  y is  a solution  of  (4),  and  the  a (x)  and  B (x)  are  functions  of  x 
which  are  independent  of  A.  Recurrent  differentia]  equations  for  these 
functions  may  be  obtained  by  substituting  (22)  in  (1),  and  equating 
coefficients  of  like  powers  of  A.  This  approach  was  used  by  Langer  (1949). 

A second  kind  of  formal  solution  corresponds  to  4.2(9).  'ibis  solution 
is  of  the  form  (3),  except  that  instead  of  cf>{x)  we  have  a function  <f>ix.  A) 
which  depends  on  A,  and  possesses  a formal  expansion 

(23)  2 .?!.  (x)A‘", 

n = 0 

where  is  the  differentiable  real  solution  of  4.5(7),  and  recurrent 

differential  equations  for  0^,  ...  are  obtained  by  substituting 
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(34)  Aii-k^  4>) 

in  (1),  and  equating  coefficients  of  powers  of  A.  Such  solutions  were 
used  by  Cherry  (1950). 

The  differential  equation 

y " + q (,x.  A)  y = 0,  qix^  A)  = S q (x)  A*"" 

n = 0 " 

is  irore  general  that  (1),  because  it  contains  the  term  q , (x)A.  This  term 
causes  certain  additional  complications  which  were  also  discussed  by 
Longer  (1949). 

'I  he  extension  of  these  results  to  the  case  of  a complex  variable  was 
undertaken  both  by  Langer  (1932)  and  Cherry  (1950);  the  extension  to 
unbounded  x was  investigated  by  Cherry  (1950). 


4.8.  Uniform  asymptotic  representations  of  Bessel  functions 

We  conclude  with  an  application  of  the  results  of  the  preceding  sec- 
tion to  the  differential  equation  4.4(3),  thereby  obtaining  an  asymptotic 
representation  of  <7^(Ax)  which  holds,  uniformly  for  all  positive  x,  as 
A -►  oo.  Re  A > 0.  The  result  to  be  obtained  includes  the  results  of  sec. 
4.4  as  far  as  they  relate  to  J^\x\  and  in  addition  fills  the  gap,  b <x  < a 
(6  < 1 < a)  left  by  the  work  of  sec.  4,4. 

Let  - us,  then,  apply  the  methods  of  the  preceding  section  to  the 
differential  equation  4.4(3),  where 

1 1 

(1)  p(a;)=l 2 , • 

X 4x 

The  transition  point  is  at  x s 1.  The  function  <f>  is  determined  by  the 
differential  equation 

(2)  = , 

X 


and  4.7(2)  becomes 

(3)  - [-<^(*)]®^  = J ' (t'*  - 1)''*  dt 
3 r. 


= -a-x")«  + iog 


l + (l-x*)’< 


4 


X 


-/3(x) 


0 <x  < 1 
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(4)  - 1)’^-  cos-'  *-’=/(*) 

3 ^ 


i 


r 


1 < X < , 


see  also  4.4(4)  and  (14). 

Since  0(x)  is  an  analytic  function,  and  (x)  ^ 0,  clearly 

1 . . 1 .i'"  3 / d>"Y  1 

/’(x)  = -l^,xl-r(x.A)=-^--(^J 


is  a continuous  function  of  x for  0 < x < oo.  To  investigate  the  behavior 
of  F (x)  as  X ->  0 or  X ->  oo,  we  use  the  chain  rule 


10,  x!  10,  m| 


+ lu,  x} 


for  the  Scfawarzian  derivative,  with  u — when  x < 1,  and  u 
X > 1.  By  a simple  computation 


F(x) 


5/3'* 

18/3* 


5/3'  * 
18/3* 


4 + x* 
4(l-x*)* 


Fix) 


Sf  * 

18/'* 


4 + X* 
4(1 -X*)* 


= f when 

0 <x  < 1 

1 < X < oo. 


From  (3)  and  (4), 

/S  O (log  x),  /3  ' =0(x-')  as  x-*0 

/■=  0(x),  /■'  «C?(1)  as  X -♦  oo, 

and  it  follows  that 

(5)  Fix)  = Q [(x  log  x)-*]  as  x -»  0 
F (x)  = O (x”*)  as  X -♦  oo. 

Furthermore,  as  in  sec.  4.7, 


(6) 


lAl0'(t)  10(01  = 


C 


except  near  zeros  of  y^(x). 
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We  thus  see  that  we  may  take  x = 0 as  the  fixed  limit  in  the  integral 
equation  for  (x).  For  the  kernel  we  have  estimates  from  (5)  and  (6) 
which  show  that  the  kernel  is  integrahle  in  (0,  <»).  'I he  method  of  sec. 
4.7  may  he  applied  to  our  problem  in  spite  of  the  singular  point  at  x = 0 
of  the  differential  equation.  Equation  4.7(14)  with  m = 0 holds,  and  since 
it  follows  from  (5)  and  (6)  that 


Ypit) 

YpM 


K{x,  t)F{t) 


dt 


V Jo  £(iogt)y  viogx/ 


as  a;  -»  0, 

we  see  that  for  small  x,  the  term  0(A~*)  in  4.7(14)  may  be  strengthened 
to 

of i \ 

\ X log  X J 

We  then  conclude  that  the  differential  equation  4.4(3)  possesses  a 
soli  . Ion  Yq  (*)  for  which 

(8)  yo(x)  = <f,)  [1  + 0(A-')], 


holds  for  0 < X < oo,  A -♦  oo.  Re  A > 0,  except  that  in  the  case  of  a positive 
real  A,  and  x > 1,  the  error  term  needs  some  modification.  For  0 < x < b 
< 1,  the  error  term  may  be  strengthened  to  (7). 

In  order  to  identify  y^  in  terms  of  Bessel  functions,  let  us  take 
0<x<h<l  in  (8),  so  that  0(x)  < <^(6)  < 0,  the  Airy  function  may  be 
replaced  by  its  asymptotic  representation  4.6(11),  the  0-term  may  be 
strengthened  to  (7),  and  we  obtain 


where  a is  the  function  defined  in  4.4(5).  A comparison  with  4.4(6)  and 
4.4(13)  now  shows  that 


2 it'/* 

r(A+  1) 


J^Xx)  = 


^-A  ,-1/1 


0 < X < 


4^ 
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and  hence 

2_«/* 

(9)  — e-^  <^)[l  + 

1 VA  1/ 

with  the  same  remarks  about  the  error  term  as  in  (8).  By  applying  Stirling’s 
formula  to  r(A+  1),  this  result  may  be  put  in  the  simpler,  if  weaker  form 

(10)  J^iXx)^ 

uniformly  in  *,  0 < * < oo,  as  A -*  eo.  Re  A > 0,  except  that  the  error  term 
needs  some  modification  near  zeros  of  Ai{—X^^  <f>).  Note  that  in  (10) 
the  error  term  contains  the  error  of  Stirling’s  formula  and  cannot  be 
strengthened  for  small  x. 

In  the  process  of  deriving  (7)  we  have  seen  that  our  present  result 
includes  4.4(1).  Let  us  show  that  (7)  also  includes  the  sum  of  the  two 
equations  4.4(2).  lo  do  this,  we  assume  x > a > 1,  <^(x)  > ^(a)  > 0,  and 
apply  4.6(12)  to  show  that 

Jp^(Ax)~  cos[A/(x)-  n/4] 

1 < a < X < oo 
or 

c7^(Ax)  ~ n Ax^  a(x)  cos  [A/(*)  - 77/4]  1 < a < x < «>, 

and  this  is  in  agreement  with  4.4(2). 

The  main  result  of  this  section,  (6),  has  b6en  extended  to  complex 
values  of  x,  and  approximations  of  higher  order  have  heen  obtained  by 
Cherry  (1948). 
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